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PREFACE. 

MY  object  in  the  following  treatise  is  to  apply  to  the  Kinetic 
/Theory  of  Gases  a  method  of  analysis  different  from  that 
generally  employed.  It  has  been  treated  always  on  a  certain 
fundamental  assumption,  namely  ,  that  the  moWn^s  of  a.  gaa 
are,  aj.  regards  their  relative  motion,  independent  of  one 
another.  As  a  consequence,  we  may  say  as  the  expression, 
of  that  independence,  the  law  of  distribution  of  momenta 
assumes  the  exponential  form  e~A^,  and,  so  far  as  concerns 
translation  velocities, 


m  being  mass,  and  u,  v,  w  component  velocities.  ^/From   this 
independence  and  from  this  form  of  Q  are  deduced  Boltzmann's     ,, 
theorems,  namely  the  H  theorem,  and  that  of  tka-eqiiality  ST 
of  mean  kinetic  energy  for  each  degree  of  freedom. 

I  propose  to  give  to  Q  the  more  general  form  of  a  quadratic 
function,  namely 

Q  =  2m  (v?  +  vz  +  w2)  +  226  (uu*  +  w'  +  ww). 
Here  b  is  a  negative  function  of  the  distance  r  at  the  instant 
considered  between  the  two  molecules  whose  velocities  are 
u,  u',  etc.,  whichouiiction  is  inappreciable  except  for  very  small 
values  of  r.  *\  shall  endeavour  to  prove  in  Chapters  IV.,  V. 
that  without  the  b  coefficients  the  motion  cannot  be  stationary. 
It  has  been  proved  abundantly  that,  assuming  the  independence, 
the  motion  is  stationary  with  the  usual  form  of  Q.  I  question 
the  axiom,  not  the  demonstration. 

The  consequence  of  attributing  to  Q  the  proposed  new  form 
is,  that  molecules  near  to  each  other  have  on  average  a  motion 
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in  the  same  direction.  They  tend  to  form  streams.  That 
result,  if  it  can  be  established,  is  worth  investigation. 

For  ordinary  gases  under  ordinary  conditions  the  b  co- 
'efficients  are  probably  very  small,  and  their  effect  negligible 
in  such  investigations  as  those  of  Tait  and  Boltzmann  con- 
cerning diffusion,  viscosity,  etc.  But  I  think  that  the  law 
*~AQ  11?  ittft  fl-lte-rad  form  will  express  the  state  of  the  system 
without  restriction  as  to  density,  except  as  follows.  A  physical 
limit  there  must  be,  when  the  gas  liquefies  under  pressure, 
if  not  before.  For  it  will  not  be  contended  that  the  dis- 
tribution of  momenta  among  the  molecules  of  the  liquefied 
_:J/  gas  is  represented  by  the  same  exponential  form  as  in  the 
gaseous  condition.  An  analytical  limit  there  is,  when  Q  in 
its  altered  form  ceases  to  be  necessarily  positive,  that  is  when 
the  determinant  of  the  coefficients  ceases  to  be  positive.  It 
'can  be  shown  that  this  determinant  does  diminish  as  density 
increases,  or  temperature  diminishes.  But  I  have  not  cal- 
culated its  value.  It  is  therefore  no  more  than  a  conjecture, 
xL  though  perhaps  a  plausible  conjecture,  that  the  vanishing  of 
the  determinant  may  coincide  with  the  physical  change  in 
the  substance. 

It  appears  to  me  that  the  law  of  equality  of  mean  kinetic 
energy  for  each  degree  of  freedom  cannot  be  reconciled  with 
my  proposed  form  of  Q\  that  in  fact  the  law  holds  only  for 
the  limiting  case  of  a  very  rare  gas. 

It  is  no  light  thing  to  question  a  conclusion  maintained 
bv^Boltzmann,  if  indeed  he  does  maintain  this  conclusion  for 
''all  substances,  or  for  all  gases  irrespective  of  density.  I  can 
but  state  the  objections  to  this  theorem,  and  to  a  certain  aspect 
of  the  H  theorem,  as  they  appear  to  me.  The  reader  will 
judge  what  weight  is  to  be  attributed  to  them. 

S.   H.   BURBURY. 

10  May,  1899. 
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CHAPTER  I. 


OUTLINE   OF  THE  THEORY. 

1.  A  GAS  according  to  the  Kinetic  Theory  consists  of  a 
great  number  of  molecules  in  rapid  motion.     And  the  object  of 
the   theory  is   to  explain  on  this   hypothesis   certain  of  the 
physical  properties  of  gases. 

Any  quantity  of  gas  which  can  be  isolated  for  the  purpose 
of  experiments  is  to  be  regarded  as  containing  a  number 
of  molecules  practically  infinite.  It  is  not  possible  to  control 
or  to  observe  the  motions  of  individual  molecules.  But  the 
theory  assumes  that  such  motion  is  subject  to  the  usual 
dynamical  laws.  Also  that  if  the  gas,  as  an  aggregate  of 
molecules,  be  at  rest,  no  dissipation  of  energy  takes  place. 

2.  A  molecule  may  consist  of  one  or  more  than  one  atom 
according  to   the   chemical   constitution  of  the  substance  to 
which  it  belongs.     It  may  be  that  hereafter  we  shall  be  able  to 
explain  on  dynamical  principles  the  chemical  relations  of  atoms 
as  constituent  parts  of  a  molecule,  and  of  molecules  inter  se. 
And   some   progress   has   been   made   in   this   direction.      At 
present  the  theory  is  concerned  not  with  the  chemical  pro- 
perties, but  with  those  properties  of  gases  which  may  change 
without  any  change  taking  place  in  the  chemical  composition 
of  the  gas:  for  instance,  density,  pressure,  and  temperature. 
And  as  depending  on  these  latter,  it  is  concerned  with  the 
phenomena  of  viscosity,  diffusion,  and  conduction  of  heat  or 
electricity. 
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3.  A  molecule  moving  as  such,  one  and  indivisible,  has 
three  degrees  of  freedom  in  respect  of  its  motion  of  translation 
in  space.     And  it  may  conceivably  have  three  other  degrees  of 
freedom,  namely,  motion  of  rotation  about  its  principal  axes. 
In  respect  of  its  internal  constitution  it  may  have  many  degrees 
of  freedom,  but  we  know  nothing  of  the  internal  constitution 
of  molecules  except  what  chemists  tell  us. 

4.  Any  finite  quantity  of  gas  may,  if  protected  from  dis- 
turbance   from   without,   retain    for    an   indefinite    time    un- 
impaired its  physical  properties,  for  instance,  its  pressure  and 
temperature,  and  also  its  chemical  constitution.     If  therefore 
these  properties  are  not  changing,  and  are  to  be  explained 
on  dynamical  principles  as  resulting  from  the  motion  of  the 
molecules,  such  motion  must  be  stationary  motion. 


OF   HOMOGENEITY. 

5.  The  gas,  or  system  of  molecules,  may  be  in  a  certain 
sense,  i.e.  if  we  regard  all  the  molecules  at  any  instant  within 
a  finite  space,  homogeneous.  For  instance,  a  cubic  centi- 
metre of  air  has  sensibly  the  same  properties  in  all  respects 
whether  taken  from  one  part  of  a  horizontal  table  or  from 
another  part.  It  is  true  that  it  will  generally  not  have  in  all 
respects  the  same  properties  if  taken  from  near  the  floor,  as  if 
taken  from  near  the  ceiling.  It  will,  namely,  be  heavier,  i.e. 
will  contain  more  molecules,  in  the  former  than  in  the  latter 
case.  That  is  because  the  air  in  the  room  is  in  a  field  of 
uniform  force,  gravity.  But  this  force  will  not  necessarily 
affect  any  other  property  of  the  gas,  at  all  events  under 
ordinary  conditions  to  which  it  is  subject. 

If  we  assume  the  infinite  divisibility  of  matter,  or  that  any 
space,  however  infinitely  small,  contains  an  infinite  number  of 
molecules,  the  condition  of  homogeneity  may  apply  to  the 
molecules  contained  in  any  such  infinitely  small  space.  We 
may  in  that  case  say  without  ambiguity  that  the  gas  at  a 
point  P  has  the  same  properties  as  at  another  point  Q. 
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6.  Modern  physicists  however  teach  us  that  the  molecules 
have  finite  dimensions,  or  that  there  is  a  finite  distance,  how- 
ever small,  within  which  two  molecules  cannot  approach  each 
other.     That  being  the  case,  the  homogeneity  can  be  asserted 
only  in  the  following  form.     If  we  consider  a  quantity  of  gas 
within   a  sphere  of  radius  r  described  about   P  as   centre, 
and  a  quantity  within  an  equal  sphere  described  about  Q  as 
centre,  either  sphere  containing  a  number  of  molecules  practi- 
cally  infinite,   then   r  being  within   certain   limits,   the   first 
quantity  has  the  same  physical  properties  as  the  second. 

OF  THE   DENSITY   OF   A  GAS. 

7.  The  same  doctrine  of  the  finite  size  of  molecules  causes 
considerable  difficulty  in  defining  certain   other  properties  of 
the  system,  for  instance  the  density. 

The  density  of  the  gas  is  usually  defined  to  be  the  sum  of 
the  masses  of  all  the  molecules  in  unit  of  volume.  But  if  the 
molecules  have  finite  dimensions,  that  may  make  the  density 
at  any  point  relative  to  that  at  other  points  depend  on  the  unit 
arbitrarily  chosen.  It  is  not  therefore  a  complete  definition. 
We  can  however  say,  and  it  is  generally  sufficient  for  our 
purpose, 

(a)  If  there  be  N  molecules  within  a  definite  volume  S, 
the  average  density  at  any  point  within  S  is  N/S. 

(b)  Again,  we  can  define  the  density  at   a   point  which 
answers  some  general  description,  or  as  we  may  call  it  the  relative 
density.     For  instance,  suppose  /  to  be  a  continuous  function 
of  the  position  of  a  point,  then  we  may  define  the  density 
with   reference  to  f,  as  for   instance   at  points   where  /=  0, 
provided  that  there  exists  an  infinite  number  of  such  points. 
For  consider  any  n  of  them,  and  about  each  a  small  sphere 
of  volume  o>.     Let  no)   be   the   unit   of  volume.     Let  N  be 
the  number  of  molecules  supposed  of  equal  mass  which,  or  the 
centres  of  inertia  of  which,  are  within  some  one  or  other  of 
these  small  spheres.     Then  the  density  where  /=  0  is  N/nto  in 
the  limit,  as  n  is  increased  indefinitely,  new  remaining  constant. 

1—2 
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To  define  density  at  a  point,  as  a  function  of  the  coordinates 
of  that  point,  is  with  finite  molecules  a  more  difficult  matter. 

Let  r  be  a  line  of  any  arbitrary  fixed  length.  As  a  first 
step  to  our  definition,  let  us  say  the  density  at  a  point  P,  or 
x,  y,  z  is  the  number  of  molecules,  supposed  of  equal  mass,  con- 
tained within  a  sphere  of  radius  r  described  about  P  as  centre, 
divided  by  the  volume  of  that  sphere.  If  n  be  that  number,  p 
the  density,  then  our  definition  is 

v      %n 
p  at  P  =  -j  — 


But  that  makes  p  generally  a  function  of  r  as  well  as  of 
x,  y,  z.  How  are  we  to  get  rid  of  r  ?  Boltzmann  and  other 
writers  after  or  before  him  would  complete  the  definition  thus, 

p  at  P  =  j  —  when  r  becomes  infinitely  small,  any  space  how- 

ever infinitely  small  being  assumed,  for  this  purpose  only,  to 
contain  an  infinitely  great  number  of  molecules.  This  defi- 
nition does  not  materially  differ  from  that  which  we  should 
give  if  our  gas,  instead  of  consisting  of  discrete  molecules,  were 
a  continuous  fluid,  any  portion  of  which  may  be  conceived  as 
infinitely  subdivided. 

But  it  is  inconsistent  with  the  teaching  of  these  same 
writers,  who  maintain  that  molecules  have  finite  dimensions, 
that  in  fact  we  can  ascertain  limits  between  which  those  dimen- 
sions must  lie.  It  may  be  admitted  that  no  ill  consequences 
have  yet  arisen  from  the  use  of  this  inconsistent  definition, 
and  if  it  be  either  necessary  or  useful,  the  inconsistency  need 
not  be  regarded  as  a  serious  objection.  Inconsistency  may 
indeed  be  an  advantage  in  mathematics  as  it  often  is  in  politics. 
It  may  be  admitted  also  that  this  definition  occurs  to  us  by 
instinct  when  we  are  asked  to  define  density  at  a  point.  That 
may  be  because  we  have  passed  through  a  former  state  of 
existence  in  which  matter  was  supposed  to  be  continuous,  and 
therefore  the  definition  strictly  accurate.  And  our  instinctive 
resort  to  it  is  due  to  unconscious  memory  of  that  former  state. 

8.  We  may  avoid  the  difficulty  by  making  a  further 
assumption  with  regard  to  the  properties  of  our  system.  Let 
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r0  be  the  radius  of  a  sphere  described  about  P  as  centre  which 
contains  a  very  great  number  of  molecules.  Let  rx  be  greater 
than  r0  in  a  finite  ratio.  Let  as  before  n  be  the  number  of 
molecules  within  a  sphere  of  radius  r.  Now  let  us  assume,  as 
a  property  of  our  system,  that  for  all  values  of  r  between  r0  and 
r1}  the  ratio  n:  r3  is  constant.  Then  we  may  define  p  at  P  as 

r— —  using  this  constant  value  of  n/r3. 

This  method  implies  either  that  the  state  of  our  gas  does 
not  change  from  point  to  point  in  space,  or  that  it  changes 
so  gradually  as  to  be  sensibly  constant  throughout  the  space 
between  the  spheres  of  radii  r0  and  TV 

9.  Assuming  that  the  state  of  our  gas  changes  so  rapidly 
in  space   that  we   cannot   make   the   above   assumption,   and 
assuming  that  we  scruple  to  base  our  definition  on  an  assertion 
which  we  elsewhere  contradict,  it  seems  to  me  that  the  ex- 
pression "density  at  a  point"  or  "number  of  molecules  per 
unit  of  volume  at  a  point "  has  of  itself  no  meaning  whatever 
when  the  molecules  are  of  finite  dimensions,  but  a  meaning  may 
be  given  it  in  an  arbitrary  way.     One  way  would  be  to  take  an 
arbitrarily  chosen  volume  for  the  unit.     I  would  suggest  the 
following  as  the  best  form  of  definition.     Let /be  a  continuous 
positive  function  of  r  which  is  equal  to  unity  for  all  values  of  r 

less  than   a  certain  distance  a,   and  for  which  -^   is  always 

r00 

negative  when  r  >  a,  and  such  that   I  r*fdr  is   a  convergent 

J  a 

series.  Then  define  the  density  at  P  as  2mf,  r  being  measured 
from  P,  m  being  the  mass  of  a  molecule  and  the  summation 
including  all  molecules  in  space.  For  with  molecules  of  finite 
dimensions,  density  at  P  must  depend  upon  the  distances  of 
many  molecules  from  P,  but  cannot  be  affected  by  those  whose 
distances  from  P  are  very  great. 

THE   STREAM   VELOCITY. 

10.  As  with  the  density,  so  precisely  the  same  difficulty 
presents  itself  if  we   try  to  define  the  stream  velocity,  or  in 
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Professor  Tait's  language,  mass  motion,  of  the  gas  in  given 
direction  at  any  point  P.  M.  Ladislas  Natanson  in  a  paper, 
to  which  I  shall  refer  later,  makes  use  of  this  function 
without  giving  any  definition  of  it  except  by  calling  it 
"  vitesse  moyenne  ou  apparente."  He  is  however,  as  I  under- 
stand, willing  to  assume  that  any  infinitely  small  space  contains 
an  infinite  number  of  molecules.  As  the  same  difficulty  arises, 
so  the  same  solution  mutatis  mutandis  may  be  resorted  to,  as  in 
the  case  of  density. 


OF   THE   PRESSURE   OF   A   GAS. 

11.  The  motion  which  we  attribute  to  our  molecules  cannot, 
assuming  the  gas  to  be  at  rest  relatively  to  the  surrounding 
space,  be  related  to  any  particular  direction  in  space.     That  is, 
at  or  near  any  point  P  there  are  on  the  average  of  time  as 
many  molecules  moving  with  given  velocity  in  any  one  direction 
as  in  any  other.    For  if  otherwise  the  pressure  of  the  gas  would, 
as   shown   by  the   next   paragraph,   be   unequal   in   different 
directions,  which  is  known  not  to  be  the  case  if  the  gas  be  at 
rest. 

12.  Assuming  the  property  mentioned  in  the  last  paragraph, 
and  assuming  the  distribution  of  momentum  among  the  mole- 
cules to  be  known,  we  can  calculate  the  quantity  of  momentum 
transferred  across  an  imaginary  plane  within  the  gas  per  unit 
of  area  and  time.     Let  m  be  the  mass  of  a  molecule,  p  the 
number  of  molecules  per  unit  of  volume,  to  be  defined  as  above 
explained,  and  pf(V)  dV  the  number  per  unit  of  volume  whose 
velocity  of  translation  lies  between  V  and  V  +  d  V.     Then  by 
our  assumption  the  number  of  these  last  whose  direction  of 
motion  makes  with  the  normal  to  the  plane  angles  between 
6  and  0  +  d6  is 


Their  velocity  normal  to  the  plane  is  Fcos  0.  Therefore  the 
quantity  of  momentum  normal  to  the  plane  transferred  across 
the  plane  in  one  direction  per  unit  of  area  and  time  by 
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molecules  having  velocity  V ...  V+dV  is 


p  V*  Psi 
Jo 


sin<9cos20d0  =  £?rc/>F2. 

And  the  quantity   of  momentum   so   transferred   by   all   the 
molecules  is 


suppose, 
F2  being  the  mean  value  of  F2  for  all  the  molecules. 

13.  If  further  the  molecules  exert  on  one  another  no  forces, 
twice  the  quantity  of  momentum  so  transferred  through  unit  of 
area  of  the  plane  is  the  normal  pressure  of  the  gas  per  unit  of 
area  of  the  plane.  See  Watson's  Kinetic  Theory  of  Gases, 
second  edition,  pp.  56 — 58. 

For  we  may  suppose  the  gas  confined  in  a  vertical  cylinder 
under  a  heavy  moveable  piston  from  which  molecules  striking 
rebound  as  elastic  bodies.  The  momentum  normal  to  the  base 
of  the  piston,  i.e.  vertically  upwards,  which  if  it  were  an 
imaginary  plane  would  be  transferred  through  it,  is  in  fact, 
the  mass  of  the  piston  being  very  great,  reversed,  so  that  twice 
that  quantity  is  transferred  to  the  piston,  and  tends  to  make 
the  piston  rise  in  the  cylinder.  That  is,  the  piston  receives 
per  unit  of  area  and  time  from  the  impacts  of  the  molecules 
momentum  equal  to  ^mp  F2  vertically  upwards.  It  also  acquires 
by  the  action  of  gravity  per  unit  of  area  and  time  momentum 
Mg  vertically  downwards,  where  g  is  the  force  of  gravity,  and 
M  bears  to  the  whole  mass  of  the  piston  the  same  ratio  which 
unit  of  area  bears  to  the  base  of  the  piston.  For  equilibrium 

Mg  =  %mp  F2  =  p,  the  pressure  per  unit  of  area. 

If  the  impacts  of  the  molecules  be  sufficiently  numerous 
per  unit  of  time  the  piston  will  remain  sensibly  at  rest. 
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OF   THE  TEMPERATURE   OF   A   GAS. 

14.  We  have  then  in  the  absence  of  intermolecular  forces 

p  varies  as  p  Vz. 

If  the  molecules  do  exert  mutual  forces  on  one  another,  the 
above  value  of  the  pressure  is  not  accurate.  But  there  are 
sufficient  reasons  for  believing  that  for  ordinary  gases  under 
ordinary  conditions  the  pressure  is  not  materially  affected  by 
the  intermolecular  forces.  For  such  gases  therefore,  accurately 
or  approximately, 

p  varies  as  p  Vz. 

But  for  such  gases,  accurately  or  approximately,  the  laws  of 
Boyle  and  Charles  hold  good,  whence  p  varies  as  p0,  6  denoting 
the  absolute  temperature.  It  is  thence  inferred  that,  accurately 
or  approximately,  8  varies  as  V2.  And  that  statement  we 
must  accept  provisionally.  Temperature  measured  from  the 
absolute  zero  is  represented  by  the  mean  square  of  the  velocity 
of  translation.  On  this  theory  temperature  is  a  quality  which 
cannot  be  attributed  to  a  molecule,  but  only  to  an  aggregate  of 
molecules. 

OF   INTERMOLECULAR   FORCES. 

15.  It  is  however  certain  that  for  sufficiently  dense  gases 
intermolecular  forces  do  exist,  and  we  must  now  make  a  further 
assumption  regarding  them.     We  shall  assume,  namely,  that  the 
intermolecular  force  R  between  two  molecules 

(1)  Acts  in  the  line  joining  their  centres  of  inertia,  or 
within  very  small  limits  of  error  may  be  regarded  as  so  acting. 

(2)  Is  a  continuous  function  of  the  distance  r  between 
those  centres  of  inertia,  which  becomes  evanescent  for  all  values 
of  r  exceeding  a  very  small  finite  distance  r0,  which  distance 
again  is  negligible  compared  with  the  dimensions  of  any  mass 
of  gas  upon  which  we  can  make  experiments. 

(3)  The  force  must  for  sufficiently  small  values  of  r  become 
repulsive,  and  should  become  infinite  as  r  is  indefinitely  dimin- 
ished— but  save  for  these  restrictions  it  may  be  that  the  force 
is  attractive  for  some  values  and  repulsive  for  other  values  of  r. 
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R  shall  be  taken  as  positive  when  repulsive. 

If  c  be  the  smallest  radius  beyond  which  the  force  exerted 
by  a  molecule  at  the  centre  on  another  molecule  is  negligible, 
c  is  defined  to  be  the  radius  of  action  of  the  central  molecule, 
and  a  sphere  of  radius  c  described  about  the  molecule  is  its 
sphere  of  action,  or  its  effective  volume. 

THE   FUNDAMENTAL   ASSUMPTION. 

16.  As  above  stated,  the  number  of  molecules  of  our  gas 
per  unit  of  volume  may  be  different  at  different  points  of  the 
system  considered.     Representing  by  p  that  number,  we  may 
have  to  treat  p  as  a  function  of  x,  y,  z.     At  present  it  is  con- 
venient to  treat  only  the  case  in  which  p  is  constant  throughout 
the   system,   and   the   distribution   of    velocities    among    the 
molecules  is  the  same  throughout  the  system.     That  being  the 
case  we  may  define  as  follows. 

At  any  instant  let  pfx  (u)  du  be  the  number  per  unit  of 
volume  of  molecules  whose  component  velocities  in  the  direc- 
tion taken  for  x  lie  between  u  and  u  +  du.  Similarly  let 
pfy  (v)  dv  be  the  number  per  unit  of  volume  whose  component 
velocities  in  direction  y  lie  between  v  and  v  +  dv,  and  pfz  (w)  dw 
the  number  whose  component  velocities  in  direction  z  lie 
between  w  and  w  +  dw. 

We  might  say,  and  it  would  be  precisely  equivalent  to  what 
we  have  already  said,  that  fx  (u)  du  is  the  chance  that  a  given 
molecule  shall  have  velocity  in  x  between  u  and  u  +  du,  and 
fy  (v)  dv,fz  (w)  dw  are  the  corresponding  chances  for  y  and  z. 

17.  These  definitions,   however  apparently  unambiguous, 
may    be   based    upon    either    of    two    diametrically   opposite 
assumptions. 

We  may  make  namely,  assumption  A. 

A.  The  chance  of  any  molecule  having  velocity  in  x  between 
u  and  u  +  du  is  independent,  not  only  of  its  position  in  space, 
but  also  of  the  velocities  v,  w,  which  it  has  in  directions  y  and  z, 
and  further  except  in  the  case  mentioned  below,  it  is  inde- 
pendent of  the  positions  and  velocities  at  the  instant  of  all  the 
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other  molecules  of  the  system.  The  excepted  case  is  when  the 
two  molecules  are  so  placed  that  they  are,  or  very  recently  have 
been,  within  one  another's  sphere  of  action.  The  force  of  this 
exception,  and  the  necessity  for  it,  will  appear  in  the  considera- 
tion of  the  H  theorem,  Chapter  ill. 

On  that  assumption  it  follows  that  the  chance  of  any  given 
molecule  having  at  any  given  instant  velocities  in  x,  y,  z 
respectively  between  u  and  u  +  du,  v  and  v  +  dv,  w  and  w  +  dw, 

is  fx  (u)  fy  (v)fz  (w)  du  dv  dw. 

Further,  if  we  consider  any  pair  of  molecules,  and  denote 
by  the  suffix  1  quantities  relating  to  the  first,  by  the  suffix  2 
quantities  relating  to  the  second,  then  the  chance  that  they 
shall  at  any  instant  respectively  have  velocities 

in  x  between  v^  and  u±  +  duly 
u2  and  u2  +  du2 ; 

in  y  between  Vi  and  v1  +  dvlt 
v2  and  v2  +  dv2 ; 

and  in  z  between  wl  and  w^  4-  dwlt 
w2  and  w2  +  dw2) 

is         flx  00/iy  Mfu  (MI)/*  (u,)  . . .  du1dvldwldu2 . . .  dw2. 


This  or  its  equivalent  is  the  assumption  on  which  the  kinetic 
theory  of  gases  has  hitherto  generally  been  treated*.     I  shall 
refer  to  the  state  of  things  assumed  &s  Condition  A.    Without 
ydoubt  it  may  be  assumed   legitimately  with  respect  to   the 
*  limiting  case  of  an  infinitely  rare  gas,  that  is  one  in  which 

*  Dr  Watson  (Kinetic  Theory  of  Gases,  2nd  Edition,  p.  3,  equation  B) 
assumes  the  independence  of  the  chances  for  two  spheres  approaching  collision. 
And  that  is  sufficient  for  his  purpose.  Professor  Tait  ("  Foundations  of  the 
Kinetic  Theory  of  Gases,"  Transactions  R.  S.  Edinburgh,  1888,  Art.  21)  makes 
an  equivalent  assumption.  Dr  Boltzmann  (Vorlesungen  iiber  Gas  Theorie, 
Part  I,  p.  23)  assumes  that  the  motion  is,  and  ever  continues  to  be, 
"  molecular-ungeordnet."  I  shall  consider  later  (Chapter  in.)  what  may  be  the 
effect  of  that  assumption.  In  the  meantime  the  weapon  that  Dr  Boltzmann 
draws  from  that  armoury  is  precisely  the  same  as  the  one  used  by  Watson  and 
Tait,  namely,  that  for  all  pairs  of  spheres  or  molecules  approaching  collision  the 
chances  are  independent,  as  stated  in  my  condition  A.  These  writers  do  not 
define  the  exception,  the  occasion  for  definition  not  having  arisen,  but  they  do 
not  assert  the  independence  in  the  excepted  case. 
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the  dimensions  of  any  molecule,  or  if  it  be  a  centre  of  force, 
its  radius  of  action,  are  negligible  compared  with  the  mean 
distance  between  any  molecule  and  its  nearest  neighbour. 
It  is  probable  also  that  all  the  known  gases  under  ordinary 
conditions  of  pressure  and  temperature,  approximate  closely  to^/ 
the  limiting  condition  of  infinite  rarity.  And  therefore  the 
calculations  made  on  the  assumption  of  condition  A  or  its 
equivalent  by  Tait,  Boltzmann,  and  others,  of  the  rate  of 
diffusion,  conduction  of  heat  or  electricity,  &c.,  for  such  gases 
remain  unaffected. 

18.     Or  we  may  make  assumption  B.  _ 

B.  The  chance  of  a  given  molecule  having  at  any  instant 
assigned  velocities  is  not  independent  of  the  positions  and 
velocities  of  all  the  other  molecules  at  the  instant.  On  this 
assumption  B,  instead  of  deducing  the  chance  of  the  members 
of  a  group  of  n  molecules  having  respectively  at  any  instant  the 
velocities 

Ui . . .  u^  +  dur^  &c. 

from  the  assumed  chances  for  individual  molecules,  we  must 
reverse  the  process.  Let  the  chance  that  at  the  given  instant 
the  coordinates  of  the  n  molecules  shall  lie  respectively  between 

x1  and 

yl  and  y^  + 

2l  and  Z-L  + 
for  the  first  molecule, 

#2  and 

&c. 
for  the  second,  and  so  on, 

and  their  component  velocities  between 
u±  and  M!  +  du^ , 
Vj  and  Vi  -f  dvl , 
wl  and 

for  the  first  molecule, 

U*  and 

&c., 
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for  the  second  and  so  on,  be  denoted  by 

F(a!Jy1 . . .  ZnU^  ...Wn) dx^  . . .  dwn. 

Then  the  chance  that  the  x  velocity  of  the  first  molecule 
shall  lie  between 

^i  and  MJ  -f  du-L , 

whatever  be  the  positions  and  velocities  of  the  other  n  —  1 
molecules,  is 


r/Y+oo 

1 1 1       (dx^dx^ . . .  dzndv1dwldu2 . . .  dwn) F(scl'i/l  ...  u± ...  wn). 

JJJ    -00 

19.    By  way  of  illustration  let  us  suppose,  under  condition  A, 


where  C  is  the  usual  constant,  chosen  so  that 

r00 
C        e~hu*  du  =  1. 

J  -00 

Then  on  assumption  A  the  chance  that  the  members  of  the 
group  of  n  molecules  supposed  to  be  of  equal  masses,  shall 
respectively  have  their  coordinates  and  velocities  between  the 
above  limits  is  of  the  form 

here  Q  =  X  (w2  +  v2  +  w2),  and  h  is  constant ; 

and  in  dealing  with  assumption  B,  let  us  assume,  the  molecules 
being  all  of  the  same  mass  and  the  same  structure, 


where  Q  =  u^  +  v-f  +  ivf  +  uf  +  . . .  + 


&c. 

Here  the  b  coefficients  shall  be  functions  of  the  distance  r,  at 
the  instant  between  the  pair  of  molecules  to  which  the  suffixes 
relate,  which  functions  become  evanescent  when  r  exceeds  a 
certain  very  small  distance.  In  that  case  omitting  the  constant 


M  I       e~ 


.  .  .  dwn 

--  • 

e    2Z?»        ,  as  shown  in  the  Appendix  (g). 
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Here  D  is  the  determinant  of  n*  constituents 

2       612     613     ., 
612     2      633 
&«     &«,    2 


and  Dn  is  its  first  coaxial  minor. 

The  form  of  f(u)  under  assumption  B  is  the  same,  except 
for  the  factor  of  h,  as  under  assumption  A^/6ut  the  physical 
state  of  the  system  may  be  very  different.  /For  instance,  under 
assumption  B,  two  molecules  at  a  distance  r  apart,  so  small 
that  the  corresponding  b  is  not  negligible,  are,  if  b  be  negative, 
on  average  moving  in  the  same  direction.  According  to  Boltz- 
mann  (Vorlesungen,  p.  21)  the  motion  in  this  case  is  not 
"  molecular  ungeordnet." 

20.  Of  these  two  assumptions  A  and  B,  B  includes  A  as  a 
particular  case.  It  is  therefore  better,  if  and  in  so  far  as  it 
may  turn  out  to  be  possible,  to  establish  a  proposition  on 
assumption  B,  rather  than  on  assumption  A.  My  object  in 
the  present  treatise  is  to  show  that  it  is  generally  possible.  It 
will  appear  also  that  A  requires  us  to  assume  our  gas  to  be 
infinitely  rare,  whereas  B  imposes  upon  us  no  stringent  con- 
dition with  regard  to  the  density.  The  only  condition  in  fact 
which  appears  to  be  required  at  present  is  this :  the  coefficients 
b  must  be  such  that  for  all  possible  values  of  the  u's,  v's,  and 
w's,  Q  shall  be  positive.  Mathematically  that  is  expressed  by 
saying  that  the  determinant  D  and  all  its  coaxial  minors  must 
be  positive. 

D  will  appear  to  be  generally  a  function  of  density  and 
temperature,  and  if  ever  Z)  =  0  the  mathematical  treatment 
changes,  whether  or  not  the  physical  system  changes  its  state, 
and  if  so  whatever  the  nature  of  that  change  may  be. 

It  will  appear  also  that  if  u,  vy  w  be  the  component  velocities 
of  a  molecule  of  mass  ra,  u,  v',  w'  those  of  any  other  molecule, 

\m  (uu  +  w'  +  ww') 

is  related  to  the  Virial  defined  in  Chapter  II.  of  the  intermolecular 
forces  acting  on  m.     See  Arts.  58,  85. 


CHAPTER  II. 

CLAUSIUS'  THEOREM. 

21.  IT  will  be  necessary  to  consider  at  length  later  the  law 
of  distribution  of  velocities  or  momenta  among  the  molecules 
for  any  given  value  of  the  mean  kinetic  energy — that  is  the 
form  of  the  function  which  we  denoted  above,  Art.  12,  by/(F). 
The  following  proposition,  due  originally  to  Clausius,  is  inde- 
pendent of  that  law. 

If  m  be  the  mass  of  any  one  of  a  system  of  molecules  in 
stationary  motion,  #,  y,  z  the  coordinates  of  its  centre  of 
inertia  at  any  instant,  we  have,  2  denoting  summation  for  all 
the  molecules  of  the  system,  and  2w#,  %my,  %mz  being  zero 

d  /          dco 


d 


with  similar  equations  for  y  and  z. 
It  follows  that 


^       %V*Vf  ^    -~»r    /  ^WV,  J 

X,    F,  Z  being   the   component   forces,  whether  external  or 
intermolecular,  acting  on  a  molecule. 
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22.  We  will  DOW  however  reserve  X,  Y,  Z  for  the  com- 
ponents of  external  force  only,  and  let  R  be  the  intermolecular 
force  between  two  molecules  acting  in  the  line  joining  their 
centres,  and  let  RXJ  Ry,  Rz  be  its  components.  So  that  Clausius' 
equation  becomes 


.........  (1). 

Half  the   expression  on  the  right-hand  side   of  the   last 
equation  is  called  the  Virial.     And  separately 


is  the  Virial  of  the  external  forces,  and  —  ^^(Ryp  +  Ryy  +  Rzz) 
is  the  Virial  of  the  intermolecular  forces. 

23.  Let  one  molecule  be  at  P  or  x,  y,  z,  another  at  P'  or 
a?',  y',  z'.  Let  r2  =  (as  -  xj  +  (y-  yj  +(z-  /)2.  Then  for  the 
force  on  P  due  to  P',  R  being  taken  as  positive  when  repulsive, 


Similarly  for  the  force  on  P'  due  to  P, 

Rx  =  --  R,    RV  =  --  R,    RZ  =  —    ~~  R, 
r  r  r 

and  therefore  for  two  molecules  at  P  and  P 


And  so  the  term  in  Clausius'  equation  due  to  the  inter- 
molecular forces  is 


............  (2), 

each  pair  of  molecules  being  counted  once. 

24.  We  have  next  to  consider  the  external  forces,  whose 
components  are  X,  Y,  Z.  I  shall  follow  Van  der  Waals' 
method. 

Let  the  whole  system  be  enclosed  by  a  surface  S,  through 
which  the  molecules  cannot  pass.  Let  p  be  the  normal  pressure 
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per  unit  of  area  on  S,  directed  inwards.  At  any  point  P  on  S 
let  X,  fj,,  v  be  the  direction  cosines  of  the  normal  s  to  S  at  P, 
I,  m,  n  those  of  a  the  line  joining  P  with  the  origin.  Then  for 
a  molecule  at  P 

-  (Xx  +  Yy  +  Zz)  =  p  (l\a  +  mpa  +  nva) 
—pa  cos  (s,  a) 
da 

***7B- 

if  ds  be  an  element  of  the  normal,  measured  outwards, 

da* 


Let  us  now  assume  p  to  be  constant  over  the  whole  surface 
8.     Then 


-*/// 


throughout  the  space  enclosed  by  $,  by  Green's  Theorem, 


if  to  be  the  volume  enclosed  by  8. 
We  see  then  that  the  Virial 


or  by  (1)  f  jw  =  2  Jm  (^2  +  w>  +  w2)  +  ^^Rr  .........  (3); 

or,  if  T  be  the  mean  kinetic  energy  of  the  system  of  molecules, 


and  --  represents  the  mean  kinetic  energy  of  a  molecule. 

If  we    take   the   volume   &>    for   a  generalised  coordinate, 
supposing  it  to  change  in  magnitude  but  not  in  shape,  the 
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other  coordinates  of  the  system  may  be  xa>  ,  ya>  ,  &c.,  where 

T     dT 

x,  y,  &c.  are  numerical.     Then  f  -  =  -j-  ,   and   if  y   be   the 

A  (o      day 

potential  of  the  intermolecular  forces  jSS-Rr/o)  =  —  -~  ,  and  our 


equation  assumes  the  Lagrangian  form  p  —  -5  --  -^  ,  p  being 
the  applied  force  which  maintains  &>  constant. 

25.  If  the  surface  8  be  an  imaginary  closed  surface,  the 
system  of  molecules  in  the  external  space  being  homogeneous 
with  the  system  within  S,  the  pressure  per  unit  of  area  of  S 
is  given  by  the  above  equation, 


The  whole  Virial  or  ^pco  —  ^^Rr  is  thus  equal  to  the 
whole  mean  kinetic  energy  of  the  molecules  within  S.  It 
does  not  follow  that  the  two  parts  of  the  Virial,  namely  f  pay 
and  —  22^Rr,  are  separately  proportional  to  SJm  (u2  +  v2  +  w2}. 


26.  If  S  be  the  external  elastic  boundary  of  the  system 
through  which  molecules  cannot  pass,  we  might  calculate  p  by 
the  same  method  as  we  calculated  it  in  the  case  of  there  being 
no  forces,  namely  p  =  ^mp'V'2,  where  p'  is  the  mean  density,  and 
\mV*  the  mean  kinetic  energy,  of  molecules  near  S.  This 
result  is  not  inconsistent  with  Art.  24,  because  the  mean 
potential  energy  near  S  is,  in  the  case  now  supposed,  different 
from  the  mean  potential  energy  throughout  the  enclosed  space, 
F*  and  so~F'2  ±  F2  and  p'  $  p. 


The  value  of  22#r  for  Elastic  Spheres. 

27.  We  have  assumed  R  to  be  a  finite  force.  In  the 
limiting  case  of  so-called  elastic  bodies,  the  force  which  acts  in 
collision  between  two  such  bodies  is,  according  to  the  usual 
convention,  an  infinite  force  which  reverses  the  velocities  of  the 
colliding  bodies  normal  to  their  tangent  plane  at  collision  in  an 
infinitely  short  time.  As  a  consequence  of  the  time  being 
infinitely  short,  no  change  takes  place  in  the  coordinates  during 
B.  2 
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collision,  and  no  body  can  ever  be  in  collision  with  more  than 
one  other  body  at  the  same  instant. 

It  is  of  course  admitted  that  no  such  body  exists  in  nature. 
t  the  behaviour  of  gas  molecules  in  regard  to  their  mutual 
ion  is  probably,  at  least  in  the  case  of  rare  gases,  approxi- 
mately the  same  as  it  would  be  if  they  were  elastic  spheres  of 
this  form.  And  jhe_assumption  made  in  the  case  of  the 
spheres  that  the  kinetic  energy  of  translation  of  two  spheres  in 
collision  is  in  the  aggregate  unchanged  by  the  collision,  must 
if  molecules  are  stable  be  true  on  average  for  any  form  of 
molecule. 


28.  It  is  then  possible  to  calculate  ^Rr  for  a  system  of 
elastic  spheres,  regarding  the  force  as  finite,  and  sufficiently 
great  to  produce  the  reversal  of  the  normal  velocities  in  a 
finite  time  T,  and  treating  the  coordinates  as  constant,  to  proceed 
to  the  limit  when  the  force  becomes  infinite,  and  r  zero. 

Let  the  colliding  spheres  be  each  of  mass  m,  and  of 
diameter  c.  Let  p  be  the  number  of  spheres  per  unit  of 
volume. 

Let  T  be  now  the  mean  kinetic  energy  of  a  sphere.  And 
let  u,  v,  w,  u',  v',  w',  be  the  component  velocities  of  any  two 
spheres,  V  their  relative  velocity.  Then 

F2  =  (u  -  uj  +  (v  -  vj  +  (w-  wj, 
and  mF2  =  2m  (v?  +  v*  +  w*)  on  average 


Let  6  be  the  angle  between  the  relative  velocity  and  the 
line  of  centres  at  collision.  Then  for  each  sphere  the  normal 
component  of  relative  velocity  is  £  V  cos  6,  and  it  is  reversed  in 
time  T.  The  assumed  force  is  then  R  =  wFcos  O/T,  and  it  acts 
during  the  time  T.  The  number  of  collisions  per  unit  of 
volume  and  time  is,  given  V,  7rc2Vp.  And  the  number  for 
which  6  lies  between  0  and  6  +  d6  is 
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Also  r  =  c.     Hence  for  all  collisions  of  a  pair  of  spheres  with 
given  V 


=  27rrac3F2p  f*sin  0cos2  6dO 
Jo 


And  for  all  collisions  of  a  pair  of  spheres  with  whatever 
value  of  V, 


or 


But  each  collision  is  now  counted  twice,  namely,  once  from 
the  point  of  view  of  each  of  the  spheres  concerned.  Counting 
each  collision  once,  or  from  the  point  of  view  of  one  sphere  only, 

tSRr  =  frrt?pT  .................  .  ......  (4), 

whence  we  obtain  for  each  sphere 


and  for  all  spheres 


and  therefore  £  22-Rr  =  if  ^P  .  pT. 

I  shall  write  K  for  f  Trc3/).     If  a  be  the  radius  of  a  sphere 
which  on  average  contains  one  molecule,  evidently 

•-*!. 

or  K  is  four  times  the  aggregate  volume  of  all  the  spheres  in 
unit  of  volume. 

It  follows  from  Art.  24  that 


or  !>  =  §(l+J*)pr  .....................  (5). 


On  a  system  of  Elastic  Spheres  in   Vertical  Column. 

29.     The  following  investigation  depends  upon  the  usual 
assumption  made  with   respect   to   collisions   between   elastic 

2—2 
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bodies,  that  a  certain  momentum  is  instantaneously  transferred 
through  a  finite  space.  The  responsibility  for  that  assumption 
rests  on  the  inventor,  whoever  he  was,  of  the  theory  of  rigid 
elastic  bodies.  Consider  a  vertical  column  of  gas,  or  material 
system,  whose  molecules  are  equal  elastic  spheres,  under  a 
constant  vertical  force/. 

If  p  be  the  pressure  per  unit  of  surface,  as  the  height  of  a 
point  in  the  column  above  a  fixed  plane,  m  the  mass  of  a  sphere, 
p  the  density,  p0  the  density  at  the  base,  we  have 


also  by  (5)  p  =  %(I+± 

Whence,  if  we  make 

T  =  constant 

wefind 


or  p  = 

Again,  consider  N  spheres  crossing  the  plane  x  =  0  with  u 
for  vertical  component  of  velocity.  Of  these  some,  say  N  —  N', 
will  reach  the  plane  x=dx  without  collision.  N'  will  undergo 
collision  before  reaching  dx.  But  for  these  N'  there  will  in 
stationary  motion  be  substituted,  as  the  result  of  collisions,  N' 
other  spheres  with  the  same  vertical  component  u. 

Now,  if  the  impact  were  direct,  i.e.,  the  line  of  centres  at 
collision  vertical,  the  effect  of  the 
collision  would  be  to  substitute  for  a 
molecule  with  vertical  velocity  u  at 
height  x  a  molecule  with  the  same 
vertical  velocity  u  at  height  x  +  c, 
where  c  is  the  diameter,  or,  as  we 
may  express  it,  the  substituted  sphere 
would  gain  a  vertical  height  c,  with- 
out losing  in  respect  of  that  distance 
any  kinetic  energy  to  the  force/.  We  might  without  affecting 
the  general  motion  suppose  the  two  spheres  to  change  places 
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after  collision,  so  that  one  would  lose,  and  the  other  gain,  the 
height  c  without  change  of  kinetic  energy.  If,  therefore,  all  the 
N'  collisions  were  direct,  the  average  height  of  the  N  spheres 

(J  /Y> 

or  their  substitutes  at  the  end  of  the  time  —  would  be,  not  dx,  but 

jp 
dx  +  -^  c,  while  their  loss  of  kinetic  energy  would  be  Nmfdx. 

But  all  impacts  will  not  be  direct ;  we  must  consider  then 
the  result  of  indirect  impacts.  For  this  purpose  consider  two 
classes  of  collisions,  in  one  of  which  the  sphere  A  has  vertical 
component  u  before  collision,  and  in  the  other  A'  has  vertical 
component  u  after  collision.  The  effect  of  a  pair  of  collisions  one 
from  each  class,  is  to  substitute  A'  for  A  as  the  sphere  with  verti- 
cal component  u.  Now  let  I  (see  figure)  denote  the  vector  line 
of  centres  at  collision,  and  cos  (ul)  the  cosine  of  the  angle  which 
I  makes  with  the  vertical.  Then  in  the  first  of  the  pair  of 
collisions  the  centre  of  A  is  below  the  point  of  contact  by 
£c  cos  (ul).  In  the  second,  the  centre  of  A'  is  above  the  point 
of  contact  by  ^c  cos  (ul).  There  is  no  reason  why  the  point  of 
contact  should  be  higher  or  lower  in  one  case  than  in  the  other. 
It  will  be  on  average  at  the  same  height.  Therefore  on  average 
of  all  pairs  of  collisions  substituting  A'  for  A  with  vertical 
velocity  u,  A'  is  above  A  by 


c  cos  (ul)  =  r,  suppose. 

Let  q  be  the  relative  velocity  of  the  two  colliding  spheres. 
Then  considering  q  as  radius  vector  of  a  cone  of  axis  u,  and  I  as 
radius  vector  of  a  cone  of  axis  q,  we  get 


r  =  c  cos  (uq)  cos  (ql) 


=  f  c  cos  (uq), 

because      cos(ql)  =  f 2  cos2  0  sin  6  ddl  f  2  cos  6  sin  0  d<9  =  f . 

Let  v  be  the  absolute  velocity  of  the  sphere  whose  vertical 
component  of  velocity  is  u,  so  that 


/    \ 
cos  (uv)  =  - . 

v    '      v 
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Then  r  =  §c  cos  (uq) 


u 


=  f  c  -  cos  (vq). 

Let  i/r  be  the  velocity  of  the  other  colliding  sphere,  E  the 

angle  between  v  and  ty.     Then 

„    u  v  —  ty  cos  E 
r=|c- 

v  q 

We  have  to  multiply  this  by  the 
'v  number    of    collisions    which   N 

spheres  having   velocity  v  undergo   with    spheres   of  velocity 
i/r . . .  ty  +  d^r,  making  with  v  angles  E ...  E  +  dE  in  time  dt,  or 

dx 

— ,  and  then  integrate  for  all  values  of  ^  and  E. 
u 

Let  pf(^r)d^r  be  the  number  of  spheres  in  unit  volume 
with  velocity  ^  . . .  ^  +  dty.     The  result  is 


/, 


u  v  —  v  cos 

sin 
o  " 


=  f  ?rc3p  .  Ndx 


Therefore  at  time  dt  the  average  height  of  the  N  spheres  or 
their  successors  above  the  plane  cc—  0  is  (1  +  K)  dx. 

But  the  energy  which  allowing  for  substitutions  they  lose 
in  the  ascent  is  Nmfdx.  The  loss  takes  place  only  during 
free  path.  It  follows  that  the  loss  of  energy  due  to  the  ascent 
dx  is,  allowing  for  substitutions,  mfdx(l—/c)  per  sphere. 

Now  suppose  that  at  x  =  0  the  number  per  unit  of  volume 
of  spheres  having  ^mu?  ...d  (u2)  for  energy  of  vertical  velocity 
is 

Ke-hmu<1  d  (u?)  ........................  (a), 

where  K  is  constant.  Then,  by  what  has  been  proved  above, 
the  number  which  at  height  'dx  have  ^muz  ...  d  (u2)  for  energy 
of  vertical  velocity  is 
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and  the  number  which  at  height  dx  have 

\mu*  ...d  (i?)  -  mfdx  (1  -  K) 
for  energy  of  vertical  velocity  is 


that  is,  K€-hmu'du?   ........................  (6), 

because,  neglecting  /t2,  &c., 


and  therefore  2hmffcdx  =  —  d/c. 

The  two  groups  (a)  and  (6)  are  equally  numerous,  and 
therefore  either  can  by  ascending  or  descending,  allowing  for 
substitutions,  exactly  replace  the  other.  Now  this  is  the 
reasoning  by  which  in  the  ordinary  case,  when  K  =  0,  we  prove 
p  =  p^-*™-?*.  It  now  proves  p  =  p0e~2hm'fx6-'t. 

I  shall  refer  to  this  problem  later  for  illustration. 

30.  The  pressure  per  unit  of  surface  is  increased  in  the 
ratio  1  :  (1  +  ^K)  as  the  molecules,  from  being  material  points, 
become  spheres  with  finite  diameter  c. 

But  the  pressure  per  unit  of  surface  is  the  quantity  of 
momentum  which  is  carried  through  unit  of  surface  in  unit  of 
time,  Art.  12.  Now,  so  far  as  this  momentum  is  carried 
through  the  surface  by  molecules  during  their  free  path,  it  is 
not  altered  in  the  least  by  K  acquiring  finite  value.  The 
increase  of  the  transfer  of  momentum  consists  in  the  case  of 
elastic  bodies  in  the  process  above  explained,  namely,  the 
instantaneous  transfer  of  momentum  through  a  finite  distance 
which  occurs  on  collision. 

The  result  of  Art.  29  may  also  be  explained  thus.  The 
mean  effect  of  all  collisions  on  any  sphere  is  equivalent  to  a 
force  K,/  acting  from  the  greater  to  the  less  density. 


CHAPTER  III. 

THE   DISTRIBUTION   OR1   MOMENTA    ^MONO   THE   MOLECULES, 
ASSUMING   CONDITION   A. 

31.  THE  simplest  form  of  molecule  that  we  can  conceive  is 
an  elastic  sphere,  such  as  described  in  Arts.  28,  29. 

An  elastic  sphere,  so  conceived,  may  be  considered  as  having 

three  degrees  of  freedom  only,  namely,  motion  of  translation 

in  space.     For^the  spheres  being  supposed  perfectly  smooth, 

motion  of  rotation  about  an  axis,  if  non-existent,  will  not  be 

yfc  produced,  and  if  existing  will  not  be  altered,  by  collisions. 

Let  u,  v}  w  be  the  component  velocities  of  a  molecule,  m  its 
mass.  Then  the  molecules  being  elastic  spheres  of  mass  m,  let 
the  number  per  unit  of  volume  of  molecules  whose  velocities 
in  direction  x  lie  between  u  and  u  +  du  be  proportional  to 
€~hmu<idu.  This  is  called  Maxwell's  distribution/It  has  been 
proved  by  several  writers  that  assuming  this  distribution 
to  exist,  and  assuming  condition  A  (which  is  necessary J^Jihe 
distribution  will  not  be  disturbed  by  collisions  between  the 
spheres. 

My  object  in  the  present  chapter  is  not  to  give  another 
proof  of  this  proposition,  but  to  show  in  what  manner  condition 
A  affects  the  truth,  as  well  as  the  proof,  of  the  proposition.  I 
shall  therefore  follow  the  general  method  elaborated  by  Boltz- 
mann  and  Watson,  as  being  more  convenient  for  my  purpose 
than  other  known  proofs,  e.g.  Professor  Tait's. 

32.  Consider  two  sets  of  spheres,  one  set  having  mass  M, 
the  other  mass  m.     Let  the  number  per  unit  of  volume  of  M 
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spheres  be  pM,  of  m  spheres  pm.  Let  us  form  a  diagram  of 
velocities,  the  velocity  of  any  sphere  of  either  set  being  repre- 
sented in  magnitude  and  direction  by  a  line  drawn  on  the 
diagram  from  the  origin  0  to  any  point  P  or  x,  y,  z.  At  P 
suppose  an  element  of  volume  dxdydz.  The  number  per  unit 
of  volume  of  spheres  M  whose  velocities  are  represented  by 
lines  drawn  from  0  to  points  within  that  element  shall  be 


y,  z)  dxdydz, 
or  shortly  pMFdxdydz. 

Call  these  the  class  F  (dxdydz). 

Similarly  the  number  per  unit  of  volume  of  spheres  m  whose 
velocities  are  represented  by  lines  drawn  from  0  to  points 
within  the  element  of  volume  dgdrjdt  at  the  point  £,  77,  f  shall 

be  p-f(E, 

or  shortly  pmfd%dr)d£, 

and  we  will  call  these  the  class 


We  may  also  express  the  above  definitions  as  follows.     The 

/wiance  that  a  sphere  M  shall  belong  at  any  instant  to  the  class 

*  Fdxdydz  is  Fdxdydz.     The  chance  that  a  sphere  m  shall  belong 

at  any  instant  to  the  class  fdgdrjdZ  is 


33.     Now  assume  condition  A. 

Then  the  two  chances  are  independent.  And  therefore  the 
chance  that  a  pair  of  spheres,  M  and  m,  shall  belong,  M  to  the 
class  Fdxdydz,  and  m  to  the  class  fdgdrjdg,  is 

FfdZdrjd&xdydz. 

On  our  assumption  of  condition  A  this  is  true  even  if  the 
two  spheres  are  on  the  point  of  collision,  and  therefore  very  near 
each  other.  For  that  case  is  not  within  the  exception  defined 
in  Art.  17.  If  we  assume  B  instead  of  A,  it  may  be  true  only 
if  the  two  spheres  are  at  a  considerable  distance  from  each 
other.  But  in  this  chapter  I  assume  condition  A.  On  this 
assumption  the  number  per  unit  of  volume  of  pairs  of  spheres, 
M  and  m,  which  have  velocities  between  the  limits  aforesaid, 
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and  for  which  m  lies  within  the  element  of  volume  dco,  so  situated 
that  M  and  m  are  approaching  collision,  is  Ffd^drjd^dxdydzdco. 
If  the  element  dco  be  so  situated  that  M  and  m  are  separating 
after  collision,  the  case  is  within  the  exception  of  Art.  17,  and 
therefore  condition  A  does  not  apply.  Let  us  now  define  dco. 
Let,  namely,  R2  =  (x  -f)2  +  (>-  ^)2+  (z-  ?)2,  so  that  R  is  the 
relative  velocity  of  M  and  m. 

Let  c  be  the  sum  of  their  radii.  About  the  centre  of  M 
suppose  a  circular  area  described  with  radius  c  in  a  plane 
perpendicular  to  R,  and  let 

dco  =  7TC2Rdt. 

Then  the  number  per  unit  of  volume  and  time  of  collisions 
between  spheres  M  of  the  class  Fdxdydz,  and  spheres  m  of  the 
is 


And  the  whole  number  per  unit  of  time  of  collisions  which 
members  of  the  class  Fdxdydz  undergo  with  spheres  m  of  any 
class  is 

pMpmFdxdydz 

Each  of  these  diminishes  the  number  of  the  class  Fdxdydz 
by  one.  So  the  rate  of  diminution  per  unit  of  time  of  the  class 
Fdxdydz  due  to  collisions  between  its  members  and  spheres  m 


pMpmFdxdydz  (  ( j  fd^drjd^Trc^R. (6). 


is 


The  Effect  of  Collisions. 

34.  We  have  now  to  find  the  number  of  spheres  M  which 
by  collisions  with  spheres  m  pass  into  the  class  Fdxdydz  per 
unit  of  Volume  and  time.  Let  P,  p  be  the  centres  of  M  and  m. 

On  a  collision  the  only  thing  that  changes  is  the  direction 
of  jR.  And  that  change  depends  upon  (1)  the  angle  ty  which 
the  line  of  centres  Pp  at  the  instant  of  collision  makes  with  R, 
(2)  the  angle  <j>  which  the  plane  through  a  line  drawn  through 
the  centre  of  M  parallel  to  R,  say  Pn,  and  the  point  of  contact 
of  M  and  m  makes  with  a  fixed  plane  through  Pn. 
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The  direction  of  R  after  collision  is  in  the  same  plane 
through  Pn,  and  makes  the  angle  2^|r  with  Pn.  That  is  the 
easiest  way  to  picture  to  ourselves  the  effect  of  a  collision. 

In  order  that  a  collision  may  take  place,  it  is  necessary  that 
a  line  ps  drawn  through  p  parallel  to  R  should  cut  the  circular 
area  above  mentioned.  We  may  divide  that  circular  area  into 
elements  whose  type  is  2c2  sin  ty  cos  ^frd^frd^.  Assuming  con- 
dition A,  the  line  ps  is  as  likely  to  cut  the  circular  area  in  any 
one  point  as  in  any  other.  Therefore  the  chance  that  the  new 
direction  of  R  shall  make  the  angle  2i|r...2^+  dty  with  the 
original  direction  of  R  is  proportional  to  2  sini/r  cos  <fyd"fy,  that 
is  to  sin  2-v/rcfyr.  It  follows  that,  assuming  condition  A,  for  any 
given  direction  of  R  before  collision,  all  directions  after  collision 
are  equally  probable,  and  on  average  the  whole  kinetic  energy 
is  divided  equally  between  the  spheres. 

35.     The  direction  of  R  after  collision  is  then  a  function  of 

x,  y,  *>  f  ,  V,  6  ^,  </>• 

Let    the   new   values   of  these   variables   be  denoted  by  the 
corresponding  accented  letters.     Evidently  ty,  <£  are  unaltered, 

or  <f)f  =  <j),  -\IT'  =  ^IT. 

All  those  pairs  of  spheres  for  which  before  collision  the 
variables  x  ...  are  between 


x  ...  x  +  dx  ...  f  .  . 
and  i/r,  0  between  limits 

T/r  .  .  .  T/T  +  di/r,    (f)  ...  </>  +  d<f> 

will  after  collision  have  these  variables  between  limits  denoted 
by  the  corresponding  accented  letters. 

Now  let  pMF'dxdydz'  be  the  number  per  unit  of  volume  of 
M  spheres  which  at  any  instant  belong  to  the  class  F'dxdy'dz. 
Similarly  after  collision  m  passes  into  the  class  f'dg'drj'dg, 
whose  numbers  aiQ/'dgdri'dg. 

36.  If  after  collision  the  velocities  of  M  and  m  were 
reversed,  without  change  of  their  positions,  they  would  by  a 
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reverse  collision  pass  into  the  classes  Fdxdydz  and  fd^drjd^ 
respectively  but  with  reversed  velocities. 

It  follows  that  the  number  of  pairs  of  spheres,  M  and  m, 
which  before  collision  belong  to  the  classes  F'dx'dy'dz'  and 
f'dgdifdZ'  respectively,  and  by  collision  with  given  values  of 
ty  and  </>  pass  into  the  classes  Fdxdydz  and  fdgdtjdZ  respectively, 
is,  again  assuming  condition  A, 

pMpmF'dxdy'dz'f'd^'drj'd^  2  sin  ^r  cos  "fydtydfytfR. 

37.  Now  x',  y'y  z,  £',  77',  f  are  functions  of  x,  y,  z,  f ,  77,  f,  ^, 
</>,  and  of  those  quantities  only.  Also  the  system  of  two  spheres 
passes  with  unchanged  kinetic  energy  from  one  state  to  the 
other. 

Whence  it  follows  by  the  general  proposition  proved  later, 
Art.  50,  and  can  be  proved  independently,  see  Watson's  Kinetic, 
Theory  of  Gases,  second  edition,  p.  9,  that 

dxdy'dz'd%dri'd%r  —  dxdydzd^drfd^ (7). 

Therefore  the  number  of  pairs  of  spheres,  M  and  m,  which 
per  unit  of  volume  and  time  pass  out  of  the  classes  F'dx'dy'dz' 
and  f'd£d*fd£f  into  the  classes  Fdxdydz  and  fd^drjd^  respec- 
tively is  with  given  i|r  and  <f> 

pMpmF'f'dxdydzdgdijdZ  2  sin  i|r  cos  tyd-fydfytfR. 

And  if  now  f,  77,  f,  -\Jr,  <f)  vary,  we  find  for  the  number  of 
M  spheres  by  which,  owing  to  collisions  between  M  and  m,  the 
class  Fdxdydz  is  increased  per  unit  of  time  the  expression 


pMpmdxdydz 

=  pMpmdxdydzjjj 

But  the  number  by  which  the  same  class  is  diminished  per 
unit  of  time  by  collisions  with  the  ra's  we  found  to  be 


pMpmFdxdydz  \  \  \ 
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If  therefore  -y-  dxdydz  denote  the  change  per  unit  of  time 
(tit 

of  Fdxdydz,  due  to  collisions  of  M  with  m, 

—*  dxdydz  =  dxdydz  Iff  d^drjd^  (F'f  -  Ff)  v*R  . .  .(8). 
Similarly 

dxdydz  (F'f  -  Ff)  -ir<?R (8a). 


38.     If  F'f  =  Ff  in  all  cases  in  which  a  pair  of  spheres  M 
and  m  can  by  collision  pass  from  the  classes  F  and  /  to  the 

7  TJ 

classes  F'  and/'  respectively,  or  vice  versa,  -^-  =  0.     Similarly 

df  *^ 

dg^O.     And  this  being  true  for  all  values  of  F  and  /  the 

motion  is  stationary.     Assuming  then  condition  A,  F'f  —  Ff  is 
a  sufficient  condition  for  stationary  motion. 

The  proof  that  it  is  a  necessary  condition  is  due  to  Boltz- 
mann.  and  is  as  follows. 

Let  H  =  { [  f"88  dxdydzF  (log  F-  1) 


Then 


J  TTT  "tf 

And  substituting  for  -=-  and  -^  their  values  above  found, 
f  f  f"  dxdydz  ffr~d&rt  (((F'f  -  Ff)  (log  F  +  log/) 

./././  -oo  J J J  -oo  ././ 

=  HI  dxdydz  HI  d^drjd^  (F'f  -  Ff)  log  (Ff)  irc^R. .  .(9). 
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But  in  the  course  of  this  integration  x,  y,  z  and  f  ,  97,  f  assume 
all  possible  positions  on  the  diagram  of  velocities,  and  therefore 
assume  the  positions  corresponding  to  any  F'  and  /',  and  there- 
fore F  and  /  among  other  values  assume  the  values  F'  and  /'. 
Whence 

~j-=jj(dxdydzjjjd£dr)d];(Ff-  F'f)  log  F'f'ir<*R...(\Q). 
and  adding  together  (9)  and  (10) 


/ 
; 


And  this  is  necessarily  negative,  if  not  zero,  and  then  only 
;ero  when  F'f  =  Ff  for  every  case  in  which  two  spheres  can 
pass  by  collision  from  the  classes  F  and  /  to  the  classes  F'  and 

/'  respectively,  or  vice  versa.     But^in  ^stationary  motion  -y- 

\       d/t 

must  be  zero.  Therefore  in  stationary  motion,  given  condition 
A,  F'f  =  Ff,  is  a  necessary  condition.  A  distribution  of 
velocities  in  which  F'f'  =  Ff  for  all  cases  in  which  a  pair  of 
spheres  can  pass  by  collision  from  the  state  F  and  /  to  the 
state  F'  and/'  or  vice  versa  shall  be  called  the  normal  distri- 
bution.  Any  other  distribution  shall  be  called  an  abnormal 
distribution.  The  H  theorem  proves  that  when  the  distribution 

•      i  ^dH• 

is  abnormal  —  j-  is  negative. 

at 

This  is  known  as  Boltzmann's  minimum  theorem,  or  the 
H  theorem,  and  H  is  Boltzmann's  minimum  function,  or,  as  he 
sometimes  calls  it,  the  Entropy  function.  The  theorem  is  as 
here  given,  founded  on,  and  depends  for  its  existence  upon, 
condition  A,  without  which  we  cannot  use  the  product  Ff  to 
express  the  frequency,  or  chance,  of  the  simultaneous  occurrence 
in  collisions  of  the  two  states  denoted  by  /and  F. 

v  On  the  H  theorem.     An  objection  considered. 

39.  If  when  the  entire  system  of  elastic  spheres  has,  with 
H  diminishing,  reached  a  certain  state,  all  the  velocities  were 
simultaneously  reversed  without  change  of  the  coordinates,  the 
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system  would  exactly  retrace  its  course  with  H  increasing.  It 
is  therefore  possible,  when  a  system  is  set  in  motion,  for  H  to 
increase  in  that  motion.  This  was  at  one  time  considered  an 
ejection  to  the  theorem,  because  the  theorem  was  supposed  to 
prove  that  under  no  possible  circumstances  could  H  increase. 
It  however  does  not  profess  to  prove  that  H  always  diminishes, 
but  that,  assuming  a  certain  condition  H  must  diminish  on 
average.  For  the  assumption  underlies  the  whole  proof  as 
hitherto  given. 

Now  to  meet  this  objection,  consider  two  spheres,  M 
having  velocities  x,  y,  z,  and  m  having  f,  77,  f.  Construct 
a  spherical  surface  of  radius  c  about  the  centre  of  M,  and 
consider  the  element  of  volume 

2c2  sin  ty  cosi|r  dty  d<f>  Rdt  —  dco. 


If  m  be  within  da,  the  two  spheres  are  either  (1)  approaching 
collision,  or  (2)  separating  after  collision,  according  as  dco  is  on 
one  or  the  other  side  of  M.  Condition  A  asserts  that  the 
number  per  unit  of  volume  of  pairs  of  spheres,  M  and  m, 
belonging  respectively  to  the  classes  F  and  f,  which  are  ap- 
proaching a  collision  of  that  kind — i.e.  a  collision  with  the  given 
T|T  and  <£ — is  Ffdx...d£dco,  see  Art.  17  and  note.  But  con- 
dition A  does  not  assert  that  the  number  per  unit  of  volume  of 
such  pairs  which  are  separating  after  a  collision  of  that  kind  is 
Ffdx...d%dco,  because  that  case  is  within  the  exception  of 
Art.  17.  In  fact  if  x,  y,  z  and  f ,  77,  f  denote  the  precollision 
velocities  of  the  separating  pairs,  x',  ?/,  &c.  their  velocities  after 
collision,  the  separating  pairs  are  the  identical  pairs  which  dt 
seconds  ago  had  respectively  velocities 

x...x+dx . 
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and  were  approaching  a  collision  of  the  kind  in  question.  The 
number  per  unit  of  volume  of  the  separating  pairs  is  therefore 
by  our  assumption  of  condition  A 

Ff,  dx...dS,  da, 
or  which  is  the  same  thing, 


and  not  F'f'dx  ...d^'dw.  (We  are  assuming  at  present  that 
Ff^f  F'f'-)  Now  in  the  reversed  course  the  pairs  which  in 
the  original  course  were  separating  have  become  the  approach- 
ing pairs,  and  vice  versa.  Therefore  in  the  reversed  course  the 
number  per  unit  of  volume  of  pairs,  M  and  m,  which,  belonging 
respectively  to  the  classes  Ff  and  /',  are  approaching  a  collision 
of  the  kind  in  question,  is  Ffdx'...d£' da,  and  not  F'f'dx'... 
da),  as  according  to  condition  A  it  should  be.  We  see  then 
that  condition  A  is  not  satisfied  in  the  reversed  motion,  and 
why  not. 

The  fact  then  that  in  the  reversed  motion  H  increases  with 
the  time  is  no  objection  to  the  H  theorem  as  a  mathematical 
proposition,  because_the  reversed  motion  does  not  satisfy  the 
condition  on  which  the  H  theorem  is  based. 

40.  We  now  see  the  necessity  for  the  exception  referred  to 
in  Art.  17.  For  if  condition  A  continues  to  exist  throughout 
the  direct  motion  without  that  exception,  it  must  also  continue 
to  exist  throughout  the  reversed  motion.  Therefore  by  the  H 

theorem  —7-  is  negative  or  zero  in  both  motions.     But  if  it  be 
uit 

negative,  and  differ  from  zero,  in  one  of  those  motions,  it  must 

7  TT 

be  positive  in  the  other.     Therefore  -7-  must,  if  we  assume 
condition  A  without  exception,  be  zero  in  both  motions. 

In  other  words,  if  we  assume  condition  A  with  the  ex- 
ception, the  mathematical  consequence  of  our  assumption  is 
that  H  is  either  minimum  or  diminishing,  as  the  theorem  asserts. 

Df  we  assume  condition  A  without  the  exception,  the  mathe- 
latical  consequence  of  our  assumption  is  that  H  is  minimum 
already.     This  condition  is  satisfied  by  Maxwell's  distribution. 
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41.  We  must  now  consider  Boltzmann's  own  assumption 
that  the  motion  is  "  molecular-ungeordnet."  It  being  assumed 
that  in  the  direct  course  the  motion  is  molecular-ungeordnet,  is 
it  molecular-ungeordnet  in  the  reverse  course,  or  not  ?  I  think 
Boltzmann's  answer  to  this  question  would  be  in  the  negative, 

the  ground  that  we,  or  Maxwell's  corps  of  demons,  have 
"  ordered "  (geordnet)  the  motion  by  the  very  act  of  reversal. 
If  this  be  so,  then  "  molecular-ungeordnet "  has,  as  applied  to 
this  theorem,  precisely  the  same  properties  as  my  condition  A. 
See  Vorlesuncien,  p.  42. 

Let  us  endeavour  to  construct  synthetically  a  system  which 
shall  without  doubt  be  rnojecular-ungeordnet.     The  molecules 
being   distinguished  by  numbers,  I  ask  (say)  Dr  Watson   to 
assign  velocities  to   them   according   to   any  law  he   pleases. 
Then   I,  in   complete   ignorance  of  those  assigned  velocities, 
scatter  the  molecules  at  haphazard  through  space,  and  they 
shall  start  from  the  positions  which  I  so  give  them  with  the 
velocities  so  assigned  by  Dr  Watson.     That  is,  primd  facie,  a   . 
molecular-ungeordnet  system — in  fact  it  is  as  near  an  approach"^ 
to  chaos  as  is  possible  in  an  imperfect  world. 

Clearly  also  in  this  system  condition  A  is  satisfied,  and  the 
system  will,  or  in  all  probability  will,  move  from  its  initial/ 
position  with  H  diminishing.  Equally  clear  is  it  that  if  in  the 
initial  position  all  the  velocities  were  reversed,  condition  A 
would  be  satisfied,  and  the  system  would,  or  in  all  probability 
would,  move  off  in  the  opposite  direction  with  #  diminishing*. 
It  follows  then  that  if  the  selected  distribution  of  velocities  be 
abnormal,  the  state  in  which  we  have  placed  the  system  is  one 
from  which  the  system  can  and  probably  will  move  with  H 
diminishing,  but  into  which  it  must  move  with  H  increasing. 

*  This  statement  is  not  inconsistent  with  that  of  Arts.  39  and  40  that  — 

changes  sign  on  reversal  of  the  velocities.  For  in  those  articles  the  pairs  of 
molecules  having  velocities  x',  £',  &c.  are  separating  after  an  actual  collision, 
before  which  they  had  velocities  x,  £,  &c.  Corresponding  to  the  classes  x,  £, 
their  number  is  Ff.  In  the  synthetic  system  no  collision  has  taken  place,  and 
the  number  is  F'f. 

3 
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dH 
dt 


JTT 

If  the  system  were  to  move  through  that  state,  — =—  would  be 


discontinuous. 

It  may  indeed,  if  the  distribution  of  velocities  .happens  to  be 
the  normal  one,  be-m  Ma.Ywell!a.  distribution.  It  cannot  be  one 
of  a  series  of  states  through  which  the  system  passes  in  its 
approach  to  Maxwell's  distribution,  as  the  H  theorem  requires. 
The  system  then  which  I  have  supposed  to  be  constructed  is 
not,  as  appears  to  me,  in  Boltzmann's  sense  molecular-unge- 
t.  It  should  seem  that  the  molecular  ungeordnet  state 
must  cease  to  be  such  on  reversal  of  the  velocities. 

What  t.V>A  H  t-hfigrem  proves  then  is  this,  that  the  distri- 
bution of  velocities  expressed  by  the  equation  F'f  =  Ff  is  the 
only  distribution  which  can  be  permanent  consistently  with  the 
existence,  and  the  continued  existence,  of  condition  A  or  its 
equivalent.  It  is  assumed  without  proof,  and  as  I  hope  to 
show  in  the  succeeding  chapters  is  not  generally  true,  that 
condition  A  can  and  does  continue  to  exist. 

The  rate  at  which  H  diminishes  has  been  calculated  by 
Tait,  Watson,  and  others  for  certain  very  interesting  cases.  See 
Watson's  work  above  referred  to,  Art.  15. 

42.  Boltzmann  (see  Berlin  Sitzungsberichte,  February 
1897)  considers  that  we  have  in  thi_s  behaviour  of  the  system  a 
true  irreversible  process.  That  term  is  generally  applied  only 
to  processes  involving  dissipation  of  energy,  that  is  conversion 
of  some  other  form  of  energy  into  heat.  If  however  heat,  or 
temperature  is  to  be  explained  as  the  kinetic  energy  of  mole- 
cules or  atoms  moving  according  to  usual  dynamical  laws, 
every  process,  even  if  it  does  involve  conversion  of  other  forms 
of  energy  into  heat,  is  theoretically  reversible  ;  would  in  fact  be 
reversed  if  at  any  instant  all  the  velocities  were  reversed,  as  we 
assumed  them  to  be  in  discussing  the  H  theorem.  If  that  be 
so,  by  asserting  that  any  process  is  irreversible  we  assert  only 
that  by  no  means  within  our  power  can  we  reverse  it,  because, 
that  is,  we  cannot  control  individual  molecules.  It  is  in  this 
sense  that,  as  I  understand,  the  diminution  of  H  in  accordance 
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with  the  theorem  would  be  called  by  Boltzmann  an  irreversible 
process. 

Boltzmann  also  shows  (Vorlesungen,  pp.  58 — 60)  that  —  H 
differs  by  a  constant^  only  from  the  logarithm  of  the  chance  of 
the  system  being  in  its  actual  state,  and  therefore  as  H 
diminishes  in  absolute  value,  the  system  passes  from  a  less  to 
a  more  probable  state.  See  note  at  end  of  chapter. 

43.     A  solution  of  the  equation 
for  the  two  sets  of  elastic  spheres  is 

a*+j/Hzs)Y 

(12). 


That  is  found  by  treating  the  kinetic  energy  as  the  only 
thing  which  remains  unchanged  by  the  collision.  It  follows 
from  this  solution  that  for  the  same  sphere  x*  =  y2  =  z2,  and 
that  for  different  spheres  Ma?  =  raf2  &c.,  or  the  mean  kinetic 
energy  is  the  same  for  each  sphere  whatever  its  mass. 

44.  The  complete  solution  of  the  equation  Ff=F'f,  as 
given  by  Boltzmann,  Vorlesungen,  p.  131,  is 


where  A  is  a  function  of  the  coordinates,  and 

Q  =  m{(a-uy  +  (/3-v)*  +  (y-w)*\   ............  (13), 

u,  v,  w  being  constant  velocities  ;  F  has  the  corresponding  value. 

In  stationary  motion  -£  =  0,  that  is 


,  = 

dx  dt  da      dtdu  } 

Suppose  now  the  system  to  be  in  a  field  of  external  force 
whose  potential  is  ^,  e.g.  the  column  of  Art.  29. 

3—2 
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If  no  other  forces  act  on  the  system  than  those  derived  from 
%,  we  have 

/da     du\          dy  /1/(X 

ra    -JI--JI    =  -  i     .....................  (14)» 

\dt      dtj          dx 

and  therefore  the  equation  becomes 

*«  +  &c.  =  0    .........  (15), 

the  solution  of  which  is          A  =e~2*x    .....................  (16). 

For  instance,  if  the  molecules  be  elastic  spheres  of  infinitely 
small  diameter,  or  material  points,  between  which  no  collisions 
occur,  there  are  no  forces  except  those  derived  from  %,  and  we 
deduce  by  the  above  method,  A  =  e~2/t*. 

The  expression  now  asserts  that  the  chance  of  a  molecule 
having  velocities  a  ...  a  +  da  &c.,  and  being  in  a  position  where 
potential  of  the  external  force  is  %...%  +  dx,  is 


It  follows  that  the  distribution  of  the  velocities  is  the  same  for 
all  values  of  %,  that  is  in  all  parts  of  the  system,  and  that  the 
number  of  molecules  per  unit  of  volume  is  proportional  to 
e-2&x.  This  agrees  with  the  result  of  Art.  29,  because  in  the 
case  now  supposed  K  =  0. 

45.  If  the  spheres  have  finite  diameter  c,  so  that  collisions 
do  occur,  there  are  other  forces  acting  besides  those  derived 
from  2£,  namely,  the  infinite  forces  assumed  to  act  for  infinitely 
short  times  during  collisions.  And  it  is  no  longer  true  that 

d  ,          \          dy 

md^a-u^-dX' 

and  we  cannot  obtain  the  result 

f=  6-2hx  e-hm{(*-uy*+&c.}} 

except  by  ignoring  the  infinite  forces  altogether.  It  is  true 
indeed  that  they  act  in  equal  and  opposite  pairs,  but  no 
solution  can  be  satisfactory  which  takes  no  account  of  them. 
See  further  as  to  this  problem  post,  Art,  103, 
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46/    If  the    intermolecular   forces  are  finite,  and   have  a 
potential  ty,  we  have 

d  /  dy     d^lr 

m  ^i  (a  ~  u)  =  ~  ~T*  -  ~f-  • 
at  dx      dx 

And  our  formula  leads  to  the  conclusion  that 

.  ...................  (17), 


expressing  now  the  fact  that  as  before  the  distribution  of 
velocities  is  on  an  average  the  same  for  any  molecule,  in 
whatever  position  it  happens  to  be,  not  only  with  respect  to 
the  external  field,  but  also  with  respect  to  other  molecules  —  but 
that  the  number  of  molecules  per  unit  of  volume  in  a  position 
in  which  the  whole  potential  is  ^  +  -^  is  on  average  e"2^"1"1^. 
T  ggo-Rr*  r^as^n  why  we  should  hesitate  to  accept  this  result,,  or 
how  we  can  avoid  accepting  itr  except  that  Dr  Watson  (in  the 
work  above  referred  to,  pp.  70,  71)  hesitates  to  accept  it, 
though  I  do  not  understand  him  to  deny  its  validity.  He 
puts  the  case  that,  with  finite  intermolecular  forces,  a  molecule 
P  may  have  a  sphere  of  action,  within  which  may  be  many 
other  molecules.  Any  individual  molecule  outside  of  the 
sphere  of  action  exerts  on  P  a  force  which  is  negligible.  But 
the  infinite  number  of  molecules  outside  of  the  sphere  of  action 
may  exert  on  P  a  sensible  force.  And  he  shows  that  this  may 
be  included  in  the  fixed  centre  forces.  But  he  refrains  from 
drawing  any  conclusion  for  the  molecules  within  the  sphere  of 
action. 


BpUzmann  and  WuAsj^fis  generalisation. 

47.  Our  molecules  regarded  as  elastic  spheres  have  each 
only  three  degrees  of  freedom,  that  namely  of  motion  of  trans- 
lation. 

Let  usLjiow_at  tribute  to  our  molecules  the  most  general 
form  possible. 

Let  a  molecule  m  be  a  material  system,  whose  position  and 
•A^     state  at  any  instant  are  defined  by  n  generalised  coordinates 
ql  ...  qn,  and  the  corresponding  momenta^  ...pn. 
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Let  pm  be  the  number  of  molecules  of  the  kind  m  per  unit 
of  volume. 


Let  Pmf(q\  ••'  qnpi  •••  Pn)  dq±  . . .  dpn , 

or  shortly  pmfdqi  ...dpn, 

denote  the  number  per  unit  of  volume  of  molecules  m  for  which 
the  variables  lie  between  the  limits 

q1  ...  q!  +  dq} 

qn  ...  qn+dqlb 

Pl...Pl  +  dp    ^ ^' 


Pn  '- 

or,  as  we  shall  express  it,  are  in  the  state  c. 

We  may  express  this  in  a  different  way  that  the  chance  of 
an  m  molecule  being  in  the  state  c  isfdq1 ...  dpn. 

It  is  assumed  that  the  molecule  is,  so  far  as  its  own 
internal  forces  are  concerned,  stable.  That  is  that  if  it  be  now 
in  the  state  c,  and  if  no  forces  except  its  own  internal  forces 
act  upon  it,  it  will  after  the  lapse  of  some  time,  greater  or  less, 
find  itself  in  the  state  (c)  again,  except  so  far  as  the  coordinates 
defining  its  position  in  space,  which  are  included  in  q1 ...  qn , 
are  changed  by  its  motion  of  translation.  It  is  also  assumed 
that  /  will  not  be  altered  by  reversing  the  sign  of  pl ... pn 
simultaneously.  And  therefore  f  cannot  contain  odd  powers  or 
products  of  the  p's. 

In  the  same  region  let  there  be  another  set  of  mole- 
cules M,  each  having_r  degrees  of  freedom,  whose  coordinates 
andjnomentaare  denoted  respectively  by  Q1 . . .  Qr  and  Pl . . .  Pr. 
And  about  these  make  the  same  assumption  as  above  stated. 

Let  pM  be  the  number  of  M  molecules  per  unit  of  volume. 
Let  pMF(Q,...  <2,A  ...  Pr)dft ...  dPr, 

or  shortly  p^FdQj.  ...dPr> 
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denote  the  number  per  unit  of  volume  of  molecules  M  for 
which  the  same  variables  lie  between  the  limits 


Qr...Qr+dQr 


(C), 


Pr-..Pr 

or,  as  we  may  express  it,  are  in  the  state  C. 

48.  Let  <j>(ql...  qnQi---  Qr)  be  a  function  of  the  coordinates 
of  two  molecules  m  and  M,  such  that  when  <f>  is  positive  there 
is  no  mutual  action  between  the  two  molecules. 

Let  us  now  make  an  assumption  equivalent,  with  respect  to 
the  molecules  now  under  consideration,  to  .condition  A.  Let 

us  assume  namely  that  so  long  as  </>  is  positive  and  -^-  negative, 

the  twoj?hances  F  and/j£g_independent.  And  therefore  the 
number  per  unit  of  volume  of  pairs  of  molecules,  one  belonging 
to  the  m's  and  the  other  to  the  M 's,  for  which  at  any  instant 
the  variables  lie,  for  m  within  the  limits  c,  and  for  M  within 
the  limits  C,  is 

PmpjfFfdqi ...  dpndQi ...  dPr. 

When  <f)  passes  through  zero  the  two  molecules  pass  into 
their  sphere  of  mutual  action,  and  an  encounter  commences 
between  them.  The  chances  F  and  f  generally  cease  to  be 
independent. 

It  may  be  the  case  that  <£,  having  become  negative,  would 
if  the  two  molecules  are  not  influenced  from  without,  remain 
negative  for  all  time.  The  two  molecules,  that  is,  may,  each 
remaining  stable,  form  a  permanent  union  revolving  round 
ne  another  in  a  closed  orbit.  I  shall  not  consider  that  case, 
but  assume  that  after  a  certain  time  r,  <f>  will  again  pass 
through  zero  and  become  positive.  The  encounter  ceases  when 
<f>  again  becomes  zero.  All  those  pairs  of  molecules  for  which 

at  a  given  instant  (f>  lies  between  zero  and    -£-  dt,  -^  being 


40  KINETIC   THEORY   OF   GASES.  [CHAP. 

positive,  will  within  the  time  dt  after  that  instant  commence, 
and  as  I  am  now  assuming,  will  in  due  time  afterwards  com- 
plete, an  encounter.  The  number  per  unit  of  volume  and  time 

of  such  encounters  must  therefore  contain  the  factor  -^- . 

'  49.  But  the  two  molecules  have  between  them  only  n  +  r 
degrees  of  freedom.  We  must  therefore  use  </>  for  one  of  the 
generalised  coordinates  denning  the  position  of  the  double 
system. 

Let  us  substitute  <f>  for  qn.  Then  all  those  pairs  of  mole- 
cules, m  and  M,  for  which  at  any  instant  the  variables  other 
than  qn  lie  within  the  limits  c  and  C  respectively,  and  for 

which  at  the  same  instant  <f>  lies  between  zero  and  -J  dt,  will 

dt 

in  the  interval  of  time  dt  after  that  instant  encounter  one 
another.  The  number  of  such  encounters  which  commence  in 
time  dt  is  therefore 

pmPMFfaqi  -  -  •  dqn^dp,  ...dpndQl...  dPr  -$  dt. 

As  a  consequence  of  the  encounter,  after  the  time  T  the 
variables  will  be  found  between  the  limits 


q n-i  ...  q n-i  +  dq n_3 
PI  PI~\~  dp^ 

Pn    Pn+dpn 

for  m,  and 


(c') 


(C') 


Pr'...Pr'+dPr' 

7    I  I 

for   M,  and   d>  between   zero   and  -±-  dt,  and  the  encounter 

at 

ceases  as  (/>  again  passes  through  zero.     The  pair  of  molecules 
pass  out  of  their  sphere  of  mutual  action. 
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According  to  our  notation  the  number  per  unit  of  volume 
of  pairs  of  molecules,  ra  and  M,  for  which  the  variables  lie 

between  these  limits  c',  C',  with  -£  positive,  is 

(Lt 

...  <-..'  -  dn'dQ,'  ...  dPr'          dt. 


Conversely  if  at  the  beginning  of  the  time  r  the  variables 

7    I  / 

were  in  the  limits  c',  C',  <£',  -^  dt  with  reversed  velocities, 

they  would  at  the  end  of  the  same  time  T  be  found  within  the 
limits  c,  C  with  reversed  velocities.  We  may  call  these  reversed 
encounters.  The  number  of  reversed  encounters  less  the  number 
of  direct  encounters  per  unit  of  volume  and  time  is 

,'  .  .  .  dpn'dQ,'  .  .  .  dPr'  ^  dt 

.  -  dq^dfr  .  .  .  dpndQl  .  .  .  dPr  -    dt. 


50.     It_can_now  be  .proved  (the  system  being  conservative) 
that 

dq,'  ...  dq'^dp,'  ... 


For  the  accented  variables  <?/,  &c.  are  all  functions  of  the  un- 
accented variables  and  of  the  time  r.  And  as  T  is  the  same  for 
every  system  within  the  initial  limits,  we  have 

dS  dS 


where  S_is  Hamilton's  Principal  function  for  the  motion  in 
question.     Similarly 

di^_  dS  _ 

dq\  dq<i 
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Now  writing  d<f>  for  -     dt  and  d(f>'  for        dt, 


i  ...  dq'^dtfdpl  ...dpn'dQS  ...  dPrf 

dQi'  ...  dQr'  x  %'  ...  dpn'dP,'...  dPrf 
dQi  ...  dQr'dq1  ...  dqn_ld<f)dQl  ...  c£Qr 

dPl'dP^    dp^ 

A«tvT       7  j  *..      _//-v       > 

-  c%  d^2      rfyr 

where  the  last  factor  denotes  the  functional  determinant 

dpi    dp^         dp^ 
dql     dqz       '   d<f>  ' 

dpi 

dql  "' 


and  since 

fj  Q  fJ       '  /-/2  Cf 

-,&c., 


the  functional  determinant  may  be  written 

d*S 


In  the  same  way 
dql  .  .  .  dqn-ld<j)dpl  .  .  .  dpndQ^  .  .  .  dPr 

i...  dQr  x  dqS  ... 


^  dp, 

^  —  j~  f  j~  >  '"  > 
dq^  dqz 

and  the  functional  determinant  has,  but  for  its  sign,  the  same 
value  as  before.  Hence  the  continued  products  of  the  differ- 
entials are  numerically  equal,  that  is 


dpi  ...  dp^dQi  ...  dPr' 

l  .  .  .  dpndQl  .  .  .  dPr  ......  (18), 
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and  therefore  the  excess  of  the  reverse  over  the  direct  encounters 
per  unit  of  volume  and  time  is 

pMpm  (F'f  -Ff)dq,...  dq^dfdfr  ...dpnXdQi...  dPr, 
or 


i  ...  dPrffJ(F'f  -  Ff)dqi...  dq^d+dfr  ...  dpn  ; 
and  is  zero  if  in  all  cases     Ff  =  Ff. 

A  corresponding  equation  holds  for  ~~  .     See  further  for  this     ., 

at  — 

proposition  Watson's  Kinetic  Theory  of  Gases,  second  edition, 
Arts.  8,  10,  12,  13. 

51.  As  this  property  is  true  for  all  values  of  F  and  /with 
the  corresponding  F'  and/',  wejeejhat  if  F'f  '  =  Ffin  all  cases 
where  a  pair  of  molecules  can  by  encounter  pass  out  of  the 
classes  F,  f  into  the  classes  F',  f  respectively,  or  vice  versa, 
we  have  a  sufficient  condition  for  stationary  motion. 

The  H  theorem  of  Art.  38  can  easily  be  adapted  to  prove 
that  it  is  also  a  necessary  condition.  Further  we  may  take 
for  solution  of  the  equation  F'f  =  Ff, 

F=Ce-*T,    f=C'<r*Tf, 

where  T  is  the  kinetic  energy  of  the  molecule  M,  and  T'  that  of 
m  ;  and  if  T^J^rar1  V>p  rpHnppH  to  the  sum  of  squares  of  the 
velocities,  we  get  a  result  analogous  to  that  obtained  for  elastic 
spheres  in  Art.  43,  namely  that,  as  a  consequence  of  condition 
A,  the  mean  kinetic  energy  is  the  same  for  each  degree  of 
freedom. 

52.  It    is    necessary   here    to    point    out    what    further 
assumptions,  besides  condition  A,  we  have   made   concerning 
our   molecules   in   the   above    demonstration.      We   expressly 
assumed  in  Art.  47,  that  each  molecule  left  to  itself  is  xtnhle, 
that  is,  that  so  far  as  the  mutual  forces  between  its  parts  are 
concerned,  the  coordinates  and  momenta,  whatever  values  they 
have  at  this  instant,  will  or  may,  at  some  future  time  pass 
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simultaneously  through  the  same  values  or  infinitely  near  to 
them,  the  coordinates  x,  y,  z  defining  the  position  in  space  of 
the  centre  of  inertia  being  excepted*/ We  also  expressly  assumed 
that  if  two  molecules  pass  into  an  encounter,  they  pass  out  of 
it  again,  that  is  do  not  form  a  permanent  union.  We  have 
further  assumed  tacitly  that  each  molecule  is  stable,  not  only 
under  the  influence  of  its  own  internal  forces,  but  also  under 
the  influence  of  the  new  forces  which  may  be  developed  during 
its  encounter  with  another  molecule.  Suppose  for  instance  for 
an  ra  molecule  /=  e~hs,  where  s  is  a  quadratic  function  of  the 
momenta  p±  ...  pn  with  coefficient  functions  of  ql...qn.  A 
necessary  condition  of  stability  is  then  that  s  shall  be  positive 
for  all  possible  values  of  the  p's.  And  that  requires  that  the 
determinant  of  the  coefficients  in  s,  which  we  will  call  d,  and  all 
its  coaxial  minors,  shall  be  positive.  See  Appendix  (q\ 


If 


&c., 


then 


d  = 


Similarly  in  order  that  an  M  molecule  shall  be  stable  if  left 
to  itself,  we  have  if 


F=  e~hs,  and  S  =  A 


&c., 


D  = 


A 


and  D  must  be  positive.  And  so  long  as  the  two  molecules 
do  not  influence  each  other,  the  condition  is  that  Dd  shall  be 
positive,  which  is  necessarily  true  if  it  is  true  for  d  and  D 
separately.  But  when  the  two  molecules  are  within  the  sphere 
of  their  mutual  action,  the  corresponding  criterion  of  stability 
is  that  s  +  S  shall  be  positive  whatever  values  the  velocities  may 
have.  If  A  be  the  determinant  of  the  coefficients  in  S  +  s, 
A,  and  all  its  coaxial  minors  must  be  positive. 
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But  if  we  form  the  determinant  A  for  the  coefficients  in 
s  4-  S,  it  will  generally  during  an  encounter  or  when  <f><  0 
contain  constituents  of  the  form  pP.  The  condition  may  fail 
for  A,  although  it  held  when  </>  >  0  for  d  and  D  separately. 
So  that  the  two  molecules,  though  stable  separately,  may  be 
unstable  together.  The  encounter  may  be  the  death  of  the 
molecule. 

NOTE  ox  THE  H  THEOREM,  Art.  42. 

/      The  H  Theorem  proves  strictly  that,  assuming  the  independence  of 
the  chances  (Art.  17),  -r-  is  generally  negative.     It  is  claimed  for  it,  as  a 

physical  theorem,  that  the  diminution  of  H  is  analogous  to  the  increase  of 
Entropy  in  irreversible  processes.  To  make  that  analogy  complete,  the 

T  TT 

chance  of  -=--  being  positive  ought  to  diminish  indefinitely  as  time  in- 
creases. I  will  state  the  difficulties  that  this  view  presents  to  me,  without 
saying  that  they  are  insuperable. 

A  finite  number  of  molecules  in  a  finite  space  are  set  in  motion  at  an 
initial  epoch  at  random,  in  the  manner  described  in  Art.  41,  and  then  left 
for  infinite  time,  undisturbed  from  without,  under  conservative  forces. 

Initially  —=-  is  almost  certainly  negative.     After  time  t  the  coordinates 

and  momenta  of  every  molecule  are  determinate  functions  of  t  and  of  all 
the  initial  coordinates  and  momenta.  If  (7,  C'  be  two  successive  states  of 
the  system,  each  in  its  turn  a  consequence  of  the  same  initial  state,  in 
what  sense  is  C  less  probable  than  C'  ?  It  may  be  that  the  initial  state, 
formed  as  we  have  formed  it,  is  less  likely  to  be  C  than  C'.  But  the 
chance  of  state  C  at  time  t  is  not  proportional  to  the  chance  of  the  initial 
state  being  C.  It  is  then  improbable  that  the  initial  state  should  be  one 

in  which  -=-  is  positive.     It  is  not  proved  to  be  improbable  that  it  should 

be  one,  a  necessary  consequence  of  which  is  that  —7-  will  be  positive  at 
time  t. 

Again,  the  continued  independence  of  the  chances,  if  conceded, 
attributes  to  every  subsequent  state  of  the  system  the  same  properties 
as  the  initial  state  has,  which  would  be  right  if  the  system  were  con- 
tinually receiving  external  disturbances.  Can  it  be  conceded  for  the 
isolated  system  ?  Let  S  be  the  Principal  function,  x,  x1  type  coordinates 

of  two  molecules  at  time  t.     If  x-=x>  nearly,  are  the  momenta  -=-,  -p-r, 

dx    dj/1 

correlated   or  independent?    I  think  it  can  be  proved  that  they  are 

T  |7F 

correlated.     If  so,  the  formula  8  of  Art.  37,  expressing  —,   ,  is  inexact. 

And  the  theorem  fails  to  prove  —=-  negative,  except  when  the  factors 
F'f  -  Ff  are  great,  that  is  when  H  deviates  widely  from  its  normal  value. 


/as* 


CHAPTER  IV. 

CORRELATION   OF   VELOCITIES.      FINITE   FORCES. 

53.  I  PROPOSE  to  prove  in  this  and  the  next  chapter  that 
in  a  system  consisting  of  molecules  of  finite  dimensions  in 
statkmary  motion,  it  is  not  true  for  molecules  very  near  to  one 
other,  that  the  chances  of  their  having  velocities  between 
assigned  limits  are  independent,  as  condition  A  assumes :  but 
on  the  contrary  if  the  forces  be  repulsive,  they  tend  to 
move  on  average  in  the  same  direction ;  and  therefore  that  in 
dealing  with  their  encounters,  or  (in  the  case  of  rigid  bodies) 
their  collisions,  we  cannot  legitimately  base  our  conclusions  on 
condition  A,  except  in  the  limiting  case  when  the  density — 
i.e.  the  ratio  which  the  aggregate  volume  of  all  the  molecules 
in  unit  of  volume  bears  to  unit  of  volume, — is  infinitely  small, 
though  conclusions  based  on  condition  A  are  approximately 
accurate  for  gases  under  ordinary  conditions. 

In  the  present  chapter  I  shall  deal  only  with  the  case  of 
finite  intermolecular  forces. 

Let  m,  &c.  be  the  masses  of  a  number  of  molecules  in 
stationary  motion.  Let  x,  y,  z  denote  the  coordinates  of  m, 
x,  y',  z'  those  of  any  other  molecule  as  m'.  Let  r  be  the 
distance  from  m  to  m'.  Let  /  be  a  function  of  r  which  is 
everywhere  positive,  finite  and  continuous,  which  is  of  negative 
degree  except  for  very  small  values  of  r,  and  such  that 


r 

Jo 


4<7rr*fdr  is  finite. 


Then  for  any  m,  2/  is  the  sum  of  the  values  of/,  referred 
to  the  position  of  m  as  centre,  for  all  the  other  molecules ; 
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and  if  p  be  the  number  of  molecules  per  unit   of  volume, 

.00  7  /• 

2/  =  I    ^Trr'pfdr.     Also  2r  -4~  has  corresponding  meaning. 
J  o  **r 

Extension  of  dausius*  Theorem. 

64.     We  can  now  put  the  Clausian  equation  in  a  modified 
form  as  follows, 


for  the  function  within  brackets  has  a  determinate  value  at 
every  instant,  and  cannot  on  average  increase  or  diminish  with 
the  time,  the  motion  being  stationary. 

Again  for  the  same  reason, 


Here  the  first  S  denotes  that  every  molecule  successively  is 
to  play  the  part  of  m,  and  the  second  2  includes  all  the 
molecules  except  m  each  with  its  own  x',  y',  2',  and  each  with 
its  own  f  referred  to  the  position  of  m  for  the  time  being. 

•  Subtracting  I.  from  II.  we  obtain 


+  2m  ^  £/(*'  -  x)  +  2m  ^  2/(y'  -  y)  +  2m  ^ 

+  s.Js(^-.)|+s.Jsy-,)| 

rf-2:  ^  x  rf/" 


^-'->i=o (m). 


df    dfdr 
A%*la'  dt=Trdt- 
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And  if  X,  fj,,  v  be  the  direction  cosines  of  mm,  or  r,  drawn  from 
vti  to  m',  x'  —  x  =  \r  &c.,  and 

dr        dx'         dy'        dz' 
dt         dt          dt         dt 


and  the  last  line  of  III.  becomes 

df  /    dx         dy        dz\  (    dx         dyf        dz 
-          ^H-7*-£  +  »^rJ  [X-^  +/*-£-+  y-j- 

dt      ^  dt        dt)  \     dt          dt         dt 

doc         dy 


That  is,  since  X/A  =  \v  =  ^  =  0  on  average,  and 
X2=/,2=z/2  =  l 

on  average,  the  last  line  of  III.  becomes, 

df  /    dx         dy       dz\  /    doc          dy'        dz 
-f-  \*>-r:+»>-i7  +  v^7)(\-7r  +P-J7+V-J 

dr\dt         dt        dt)  \    dt          dt         dt 


v 

Again,          -jSr. 

-hi*^M*^ 

=  2/  by  the  conditions  for  /,  if  -^  =  0.* 

dzx 
Again,  m-^  =  X)  &c. 

if  JT,  F,  ^  be  the  components  of  all  the  forces,  whether  external 
or  intermolecular,  acting  on  m. 

*  If    the    molecules    be    scattered   uniformly  ^  =  0  on  average.     If  not 

uniformly,  it  will  be  generally  negative,  because  we  take  a  molecule  for  centre, 
and  there  are  more  molecules  where  p  is  greater.  But  (see  Art.  106)  it  will  be 
small.  If  negative  it  helps  my  present  argument,  but  I  shall  treat  it  as  zero. 
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Therefore       m  -=-£  2/V  —  x  =  X^f(x  —  x), 


m  — 
Making  these  substitutions,  and  arranging  the  terms,  we  find 

T  •{  Sm  -j-  £/(#'  —  x)  +  &c.  >  =  0 
dt  (       dt 


(dx  d^_,dydy^     dz_dz\ 
U*  dt+  dt  dt+dt  dt) 


That  is 

dt  dt      dt  dt 


'  -  y)  +  2 
dy       dz\/dx          dy'       dz1 


=  0  .........  (IV). 

55.     Now  if  there  are  no  intermolecular  forces, 


is  for  an  infinite  system  negligible,  and  equation  IV.  is  satisfied, 
whatever  /may  be,  by  making 

dx  dx'     dy  dy      dz  dz'  _ 
ti~dt+~~'¥~~'~   ' 
A  dx        dy 

and    (xdt+^dt+l' 

on  average  for  each  value  of  r,  and  therefore  for  each  value 
of  /.  This  expresses  condition  A,  which  thus  appears  as  a 
consequence  of  there  being  no  intermolecular  forces. 

a  4 
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56.  But  if  there  are  inter  molecular  forces,  we  shall  find 
that  Qondition  A  cannot  exist. 

In  the  limiting  case  when  the  encounters  are  binary,  or  no 
molecule  is  within  the  sphere  of  action  of  more  than  one  other 
molecule  at  the  same  time,  we  have 

XZf(x'  -x}+  Y^f(y'  -  y)  4-  Z2/(/  -  z)  =  -  Erf. 

In  the  general  case  consider  a  spherical  shell  of  radii 
r  ...r+  dr  described  about  m.  Let  E  be  the  intermolecular 
force  between  m  and  a  molecule  m'  within  the  shell,  0  the  angle 
made  with  x  by  the  line  mm  drawn  from  m  to  m'.  Let  fr  be 
the  value  of  /  for  the  radius  r,  and  let  2r  denote  summation 
for  all  molecules  m  within  the  shell.  Then  if  X,  &c.  relate 
only  to  the  forces  exerted  on  m  by  molecules  within  the  shell, 

Z2/(V  -  x)  =  -frEr  (2r  cos  0)2, 

the  mean  value  of  which,  if  the  molecules  be  distributed  through 
the  shell  at  haphazard,  is  —  J2r/?r. 

Again  X  includes,  not  only  the  forces  acting  on  m  due  to 
the  molecules  in  the  shell  r  ...  r  4-  dr,  but  also  all  other  inter- 
molecular  forces  acting  on  m,  e.g.  forces  due  to  molecules 
within  the  shell  r'  .  .  .  r  +  dr,  which  forces  we  will  denote  by 
R'.  The  part  of  X2f(x'  —  x}  due  to  these  forces  is 

—frEr^r  cos  02^  cos  6', 
where  Q'  corresponds  to  6  above.     If  we  assume 

2rcos02r<cos0'  =  0  .....................  (1), 

on  average,  we  shall  have 


But  whether  X^f(x'  —  x)  be  equal  to  —  ^rErf  or  not,  it  is 
nevertheless  negative,  and  not  zero.  What  has  been  proved 
for  X  is  equally  true  for  Y  and  for  Z.  So  that 

X2f(af  -  x)  +  F2/(2/'  -  y)  +  Z2f(z'  -  z) 

is   in  any  case  negative,  and   in  the  limiting  case  of  binary 
encounters,  or  wherever  we  can  assume  (1),  is  equal  to  —  ZErf. 

Since,  as  we  have  seen,  in  equation  IV.  the  term 

2Z2/(a?'  -  a)  +  &c. 

is  negative  and  not  zero,  it  follows  that  the  sum  of  the  remaining 
two  terms,  that  is, 
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?  ^f         i. 

*  J\dt  dt  **  dt  *dt  dt) 


dy        dz\  /    dx         dy 

- 


is  positive  and  not  zero.  But  this  is  inconsistent  with  condition 
A.  It  follows  that  we  must  discard  condition  A  when  inter- 
molecular  forces  exist. 

57.     First  solution  of  equation  IV.     Let  for  any  value  of  r 
(dx^  dx'     dy^  dy'     dz  ~  dz'\ 

m(di-rdt+dt*r  it  +  dt  2''  di)  =  47r^  on  avera«e' 

where  <f>  is  a  function  of  r  such  that  4^?^,  and  therefore  also 
,  vanishes  at  either  limit,  when  r  is  zero  or  infinite.    Also 


Ci  ¥    fill 

assume  all  terms  of  the  form  -^-  ~-  ,  &c.  to  be  zero  on  average. 

Then  as  in  Art.  54, 

df  /.  dx  dy  dz\  /  dx  dy'  dzf\ 
^U.-Tr4-/i-£  +  z/-TrX-T7+At^-  +  ^jTl 
dr\  dt  ^  dt  dtj\  dt  dt  dt  J 

,  because  \*  =  ~jjp  =  ^  =  ^ 

df  * 
*    dr 
dr 


and  using  this  in  IV.  we  find 

2ZS/  (x  -x)  +  &c.  -  J2  r*irrV  ^fdr  =  0, 

which  is  satisfied  by  assuming  for  every  r 


also  I   4t7rr2p(f)dr  =  —  i  I    4t7rrsp  -^  dr, 

Jo  Jo  df 

,00 

whence  |    4t7rr*p<f>dr  =  ^.Rr, 

-  o 

in  all  cases  in  which  we  can  assume  (1). 

Second  solution  of  equation  IV.     Since  /  is  arbitrary  we 

may  equate  separately  to  zero  the  term  in  IV.  involving  r  -£  , 

4—2 
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and  the  sum  of  the  two  terms  involving  f.     That  is,  we  make 
for  each  value  of  r, 

dx        dy 


in  which  X  includes  the  intermolecular  forces  acting  on  m  due 
to  molecules  at  whatever  distance,  but  x'  —  x,  &c.  relate  to  a 
particular  distance  r. 

The  result  (IV  b)  and  the  first  solution  show  that  the 
velocities  of  m'  molecules  distant  r  from  m  are  correlated 
with  that  of  m,  or 

(dx  ^   dx      dy  ^   dy'     dz  v«  dz'\ 

m  (dt  2"  *  +  ft  ^  dt  +  Tt  ~*Tt  )  >  °'  on  avera«a 

And  the  correlation  is  a  function  of  r. 

In  this  solution  47rr2/)<£,  and  in  the  former  one  —  ^4?rr3/)  -^  , 

is  a  determinate  part  of  the  Virial  of  the  intermolecular  forces 
acting  on  m,  though  not  necessarily  equal  to  Rr,  except  in 
cases  where  (1)  can  be  assumed. 

When  r  exceeds  the  "  radius  of  action  "  of  the  molecule  m, 
and  so  R  =  0,  X  nevertheless  includes  the  x  force  due  to 
molecules  distant  r'  from  m,  for  which  R  is  not  zero.  But 
ultimately  as  r  increases,  £r  cos  6  2/  cos  &  —  0  for  all  values  of 
r'  for  which  R  is  sensible.  So  that  as  r  increases, 

m  ~Ji  ^r  -JT  +  &c.  becomes  in  any  case  zero. 
at       at 

But  the  correlation  may  possibly  retain  a  sensible  value  at 
distances  from  m  greater  than  its  radius  of  action. 

Equation  (IV  a)  expresses  that  although  the  velocity  of  m 
is  correlated  with  that  of  m',  yet  it  is  not  on  average  correlated 
with  the  velocity  of  m'  resolved  in  direction  mm',  that  is  not 

correlated  with  X  ^'  +  /*  ^  +  „  —  .   For  (IV  a)  with  (IV  b)  the 
at          at         at 

simplest,  though  not  perhaps  the  best,  solution  is  the  following. 
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If  mP  represent  in  magnitude  and  direction  the  velocity  of 
m,  and  if  it  be  taken  for  polar  axis  of  a  sphere  of  radius  r 
described  about  m  as  centre,  then  the  average  motion  of 
molecules  on  the  surface  of  that  sphere  is  on  the  meridian 
towards  P,  and  proportional  to  the  sine  of  the  colatitude. 

58.  If  the  force  R  be  attractive  instead  of  being  repulsive, 
-IT  -7-  4  &c.  will  be  negative  instead  of  positive.  The  velocities 

of  m  and  m'  may  be  said  in  this  case  to  be  contrar elated.  If 
there  be  two  sets  of  molecules,  A  and  B,  and  every  A  repels 
every  other  A ,  every  B  repels  every  other  B,  but  A  and  B 
attract  each  other,  the  system  will  tend  to  form  at  every  point 
two  streams,  the  A's  moving  in  one  direction,  and  the  B's  in  the 
opposite,  which  we  may  compare, 

'  si  parva  licet  componere  magnis,' 

with  the  mutually  attracting  and  repelling  electricities  supposed 
to  constitute  an  electric  current. 

It  should  be  noted  that  for  any  molecule  m, 

'dx  dx' 


the  summation  including  all  values  of  r,  is  if  (1)  be  assumed 
equal  to  the  Virial,  taken  as  positive  when  the  forces  are  re- 
pulsive, of  the  intermolecular  forces  acting  on  m.  (See  Art.  85.) 


CM  */*   CM  ST 

I  have  treated  -^      -  as  zero  in  the  absence  of  intermolecular 
dt  dt 

forces.  Strictly,  n  the  number  of  molecules  in  the  system  being 
finite,  and  the  centre  of  inertia  at  rest,  it  must  be  negative,  but 
it  may  be  neglected  when  n  is  great. 

59.     I  have  thus  proved  that  if  there  be  finite  intermolecular 
forces,  condition  A  cannot  prevail,  and 

dx  dx      dy  dy      dz  dz 
~^~     ~+~ 
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has  on  average  a  finite  value,  a  function  of  r,  and  is  positive  if 
the  forces  are  repulsive.  The  investigation  does  not  apply  to 
elastic  spheres,  in  which  the  intermolecular  forces  are  infinite 
during  collision.  I  shall  therefore  devote  the  two  following 
chapters  to  proving  the  same  thing  for  equal  elastic  spheres  by 
an  independent  method.  It  will  be,  I  regret  to  say,  a  much 
more  laborious  method.  If  however  the  reader  is  now  satisfied 
that  I  have  established  my  case  for  finite  intermolecular  forces, 
he  will  have  little  difficulty  in  provisionally  taking  for  granted 
that  the  same  thing  can  by  some  means  or  other  be  established 
for  the  limiting  case  of  elastic  spheres,  in  which  the  forces 
become  infinite  for  very  small  values  of  r.  And  omitting 
Chapter  v.  proceed  at  once  to  Chapter  vi. 


CHAPTER  V. 

CORRELATION   OF   VELOCITIES.      ELASTIC   SPHERES. 

60.  I  PROPOSE  in  this  chapter  to  prove  for  elastic  spheres  a 
proposition  corresponding  to  that  proved  in  the  last  chapter  for 
molecules  between  which  finite  forces  act.  The  elastic  spheres 
shall  be  each  of  unit  mass  and  diameter  c.  The  same  pro- 
position for  spheres  of  unequal  masses  or  diameters  is  much 
more  laborious,  and  when  proved  does  not  bring  us  much 
nearer  to  any  probable  form  of  molecule. 

Equal  Elastic  Spheres. 

Assume  the  chance  that  a  group  of  n  spheres  shall  simul- 
taneously have  component  velocities 

A  -  •  •  A  +  <^A  [  f°r  ^ne  firs^ 
7i  . . .  71  +  <£> 


\ 

0n  ...J3n  +  dj3n  [   for  the  rcth 
7»  •••  Vn  +  djn  ) 
to  be  Ce 


If  we  assume  condition  A,  we  virtually  assume  that 
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Now  it  will  be  found  that  when  the  diameters  are  finite 
the  motion  is  not  stationary  when 


but  becomes  stationary  when 

Q  =  2  (a2  +  /32  +  72)  +  226  (««'  +  0ff  4-  77'). 
Here  a,  a',  &c.,  are  component  velocities  at  the  instant  con- 
sidered of  any  pair  of  spheres,  and  b  is  a  function  of  the  distance 
between  them  at  that  instant  which,  except  for  very  small 
distances,  is  evanescent.  I  admit  of  course  that  given  con- 
dition A  always  existing,  the  distribution  of  velocities  among 

^   the  spheres  denoted  by  e~^2  («a+0"+y*)  is  not  disturbed  by  collisions. 

ylx  That  has  been  proved  many  times.  But  the  point  is  that 
condition  A  itself  cannot  continue  to  exist  when  collisions 
take  place.  In  order  that  it  may  continue  to  exist,  it  is 
necessary  that  the  system  be  continually  receiving  disturb- 
ances, or  else  that  the  diameters  of  the  spheres  shall  be 
infinitely  small,  so  that  no  collisions  occur. 

Definitions. 

61.  Let  P  be  any  point,  r  the  distance  of  any  other  point 
from  P.  Let  /  be  a  function  of  r  satisfying  the  following  con- 
ditions, namely, 

(1)  /is  equal  to  unity  at  all  points  within  a  small  sphere 
of  radius  a  described  about  P  as  centre. 

(2)  /  is  finite,  positive  and  continuous  everywhere. 

(3)  /  is  of  negative  degree  at  all  points  outside  of  the  '«' 

Too 

sphere  and  such  that  the  series  denoted  by  2/,  or  I    4<7rrzfdr, 

J  a 

is  ultimately  convergent,  but  contains  no  high  powers  of  -  . 

We  may  take  for  'a'  the  radius  of  a  sphere  which  on 
average  contains  one  molecule. 

Let   now  f=  yf>.rl=-?f>  f—  ^j??j  ^ne    summation   in- 

cluding all  the  molecules  of  the  system,  and  /  having  for  each 
molecule  the  value  proper  to  its  distance  from  P.  Then  f,  rj,  f 
are  continuous  functions  of  xt  y,  z,  the  coordinates  of  P. 
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Let  M^jjjdxdydzzy^  +  g)  ,  in  which  (g+  gis  the 

value  of  that  function  at  a  given  instant  within  the  element  of 
volume  dxdydz,  ay  is  the  value  of  the  product  ay  for  the 
molecule  (if  any)  which  is  within  that  element,  or  the  sum  of 
those  products  for  all  molecules  within  the  element  if  there  be 
more  than  one,  and  the  integration  includes  all  space.  Then, 
for  any  element  of  volume,  ay  changes  from  instant  to  instant 
by  the  passage  of  molecules  into  or  out  of  the  element. 

There  are  two  analogous  functions,  namely, 


and  also  three  others  of  the  form 

2    &c- 


It  will  be  sufficient  in  this  chapter  to  deal  with  MXZ)  and  we 
may  for  the  present  drop  the  suffix,  and  write  M  for  Mxz. 


Summary  of  the  Argument. 

62.  It  is  necessary  first  to  treat  of  the  case  in  which  the 
molecules,  though  possessing  finite  mass,  are  of  infinitely  small 
dimensions,  so  that  collisions_between  them  will  not__j>ccur. 
That  being  the  case,  the  system  will  be  in  stationary  motion  if 
the  velocities  are  distributed  according  to  Maxwell's  law,  Q 
being  a  sum  of  squares  only.  This  we  may  admit  whatever 
views  we  may  hold  as  to  the  possibility  of  any  other  law.  I 
shall  assume  Maxwell's  distribution  to  exist.  It  follows  from 
Q  being  a  sum  of  squares  only,  that  the  fact  of  any  molecule,  or 
any  group  of  molecules  having  any  given  velocities,  affords,  if 
the  whole  number  of  molecules  be  infinite,  no  presumption 
whatever  with  regard  to  the  velocities  of  any  other  molecules. 
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We  then  find  for  this  system  the  mean  values  of  f  2,  (  -7-  J  ,  &c. 
In  this  system  of  material  points  —7-  =  0,  which  gives  rise  to 

(MI 

equations  A  and  B  below.  It  is  then  proved  that  if,  with  the 
same  masses,  the  spheres  have  finite  diameter  c,  so  that  collisions 

occur,  -=-  ,  which  was  zero,  becomes,  if  f  2,  &c.  retain  the  same 
cit 

mean  values  as  before,  positive,  and  the  motion  not  stationary. 
But  the  motion  must  become  stationary.  Therefore  in  the 
system  of  finite  spheres  f2,  &c.  cannot  retain  the  same  mean 
values  as  before.  It  follows  that  in  this  system  Q  in  e~hQ  cannot 
be  the  sum  of  squares  of  the  velocities  only. 

A  System  of  Material  Points. 

63.  The  chance  that  the  molecules  shall  at  any  instant 
have  for  their  component  velocities  in  xt  ax  ...  «!+€?«!,  &c.,  is 
A€~hQdcL!  ...  doint  where  A  is  a  constant  and  in  the  system  now 
considered 

Q  =  of  +  a/  +  &c. 

But  from  Art.  61, 


It  conduces  to  simplicity  to  assume  fn  =  1,  but,  as  will  be 
seen,  it  in  no  way  affects  our  conclusions.    That  being  assumed, 

«n  =  ?2/-  ai/i  -  "2/2  -  &c- 
Substituting  this  value  for  an  in  Q,  we  obtain 

+/22)  «22 


2/1/3a1a3  +  &c. 

If  we  now  integrate  e~hQ  for  0.1  .  .  .  «n_i  successively  with  limits 
±  oo  ,  we  get  the  result  next  stated. 

Let  D  be  the  determinant  of  the  coefficients  in  Q,  each  of 
the  coefficients  of  products  being  divided  by  2,  that  is 
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D  = 


/,/, 


Let  Dp4  be  the  minor  formed  by  omitting  the  pih  row  and 
column.      Similarly   let   Dpqr8   be   the   minor   formed   by 
omitting  rows  p  and  r,  and  columns  q  and  s.     Then  by  inte- 

grating for  «!  ...  ««_!,  Q  is  reduced  to  -=r-  f2,  whence  f2  =  ^r^. 

See  Appendix  ( 
reduced  to 


Also  by  integrating  for  0%  ...  an_!  only,  Q  is 


Again  we  find  (see  Appendix 
Ai  =  1  +/i2  +/22  +  ...  /Vi  = 

A221  =  ?/2  -/,», 

Am  =/i/», 


since  /n2  =  1, 


(19). 


We  have  then  for  the  chance  that  «!  and  f  shall  lie  within 
assigned  limits  irrespective  of  the  values  of  «2  .  .  .  a,^  the 
expression  Ae~hq'da.id%  in  which 


and 


We  see  now  that  the  assumption  fn  =  1  has  not  affected  the 

/2 

result,  because  writing  •£-  for  /2  does  not  alter  p,  q,  or  r. 

Jn 

64.     Let  now  f  be  the  value  of  that  function  at  P  the 
origin,  f  +  9£  its  value  at  a  neighbouring  point  P'.     Then  at 

d 
P'}  e-h#  becomes  €~h<?  +  d£  3fc  €~h(*-    ^d  if  f  be  very  small,  this 

is  equal  to 
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But  if  at  P,  f  is  not  negligible,  we   must  use  instead  of 
a,  /3,  7  the  following  functions,  namely, 

M  =  a  -  f ,  v  =  @  -  i),  w  =  y-£,  whence  a  =  u  +  j;  &c., 


and  Qf  =  q  (u  +  %)2-2pu  +  %  %  +  r%2, 

or  Of  =  Qu2  +  2a  —  v£u  +  (Q  - 


and  let  e^'  =  <£  (M  f  ). 

Then  if  </>  (^  f  )  be  the  value  of  this  function  at  P,  ^  (w  f  )  its 
value  at  P', 


Again,  let   <f>  (u)  du  be   the   chance   that  at  P  M  shall  lie 
between  u  and  u  +  dw,  whatever  f  may  be,  that  is 


and  the  same  chance  at  P'  is 

*.  («)  -  r 

J    — 

Also          |Mtt 
if  f          denotes  differentiation  with  u  constant.     Therefore 


because  </>  (w  f  )  vanishes  at  either  limit.     Therefore 


=  9f  f    ^^  (» 

J     -00 

And  this  is  proportional  to 
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= 


whence  fa  (u)  =  <£  (u)  (1  +  2A' 

65.     Now  with  the  above  definition  of  Jtf  we  may  write 
-IT-  =  0,  that  is 


the  summation  being  for  all  the  elements  of  volume,  each  with 

i  -.  (d£      dt\ 

its   own   ay  ana  its  own  (  -^  4-  -=-    . 

\dz      dx) 

Now  we  shall  find  that  on  average  -=-  ay  is  of  the  opposite 

at 

/ d£  .  d£\         ,   d  ( d£  .  dt\        „  .,  .         ... 

sign  to  (  -p-  +  -y2    ,  and  -=-    -p-  -f-  -r5-    is  of  the  same  sign  with 
\dz      dx)  dt\dz      dx) 

ay.    So  that  — =-  consists  of  two  series  neither  of  which  is  zero, 
at 

but  which  are  equal  and  opposite. 

66.     Let  us  consider  the  plane  of  xz,   and  f  as  varying 
in  direction  z,  f  as  varying  in  direction  x.     Let  us  first  deal 

with  ~~    only.      Consider   the   two  infinite  planes   z  =  0  and 
dz 

z  =  dz,  and  suppose  for  a  moment  ~  constant  between  them. 

dz 

The  number  of  molecules  at  or  near  the  plane  z  =  dz  for  which 
u,  or  a  —  g,  lies  between  u  and  u  -f  du  exceeds  the  correspond- 
ing number  at  or  near  the  plane  z  —  0  by  the  quantity 


. 
dz 

That  follows  from  Art.  64.  To  fix  the  ideas  let  u  be  positive. 
Let  us  separate  from  this  class  of  molecules  those  whose  z 
velocity  is  7,  and  suppose  7  negative.  Let  us  then  call  the 
separated  class  the  class  (uy).  Then  the  number  of  molecules 
of  the  class  (uy)  which  pass  through  the  plane  z  —  dz  per  unit 
of  area  and  time  exceeds  the  number  of  the  same  class  which 
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pass  through  the  plane  z  =  0  per  unit  of  area  and  time  by  the 
quantity  Ae~h'u*  e~hv  duZh'uy-^-dz.  And  the  number  of  mole- 

cules of  the  class  (wy)  which  enter  S  the  space  between  the 
planes  per  unit  of  time  exceeds  the  number  of  the  same  class 
which  pass  out  of  S  per  unit  of  time  by 


'u2  e~W  dudyZh'uy       , 

(jiZ 

which  is  positive  because  uy  is  negative. 

If  SNdudy  be  the  number  of  the  class  (uy)  within  S,  we  have 

^  dudy  =  -  Ae~h'u*  €-hv2  dudyZh'uy  ^  ......  (20). 

at  dz 

But  if  uy  be  the  mean  value  of  uy  for  all  the  molecules  within 
S,  u  and  y  now  taking  all  possible  values, 

_      /Y+0°  /  f  f+0° 

uy  =  1  1      Nuydady  /  1  1       Ndudy, 
and  therefore 


7       /*/* 

because  -j-  I  \Ndudy  =  0  on  average, 

=  -  A      e-^e-WWu^dudy       by  (20) 


I  have  assumed  in  this  proposition  that  the  space  S  within 
which  -j-    is   constant,  is   the   continuous  space  between  two 

infinite  parallel  planes.  That  will  not  occur  in  fact.  But  it  is 
not  essential  to  the  argument  that  S  should  have  that  shape, 
or  that  it  should  consist  of  one  closed  surface.  The  proposition 
holds  equally  if  S  consists  of  all  those  elements  of  volume  for 

which  at  a  given  instant  -5  has  the  same  value.     And  it  is  true 

dz 

on  average  for  each  element  of  volume. 
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d—     d 


67.     Again,  ^  wy  =  ^  (a  -  f  )  7 

d  —     d  „ 

-S**-5^ 

Now  a  molecule  moving  vvithin  $  with  positive  7  moves 
from  less  to  greater  f,  -^  being  positive;  if  with  negative 
7,  it  moves  from  greater  to  less  f.  For  all  molecules  moving 
within  S,  -=-  yg  is  positive  and  equal  to  yz  -£•  .  If  S  be  divided 

Cit  (LZ 

into  parallelepipeds  of  base  dxdy  and  altitude  PP 

§  ******  = 

f  P  and  fp>  being  the  values  of  f  at  the  upper  and  lower  surfaces 
of  the  parallelepiped.  By  the  motion  of  all  these  molecules 
7f  is  increased. 

But  at  the  upper  end  where  f  is  greater,  molecules  with 
positive  7  pass  out  of  the  space,  at  the  lower  end  they  pass  in. 
For  negative  7  the  converse  is  the  case.  By  these  entries  and 

exits  7f  is  therefore  diminished  at  the  rate  I  \<f  (%P  -  fp-)  dxdy, 

or  by  the  same  quantity  by  which  it  is  increased  by  the  motion 
of  the  molecules  within  S.  On  the  whole  therefore  if  f  alone 
vary 

' 


d  -      d  —    d  - 

And  5*»~  a       ' 


"a"* 

We  have  then  from  Art.  66, 

d  —         1  df 
dt^=   -Zhdz' 

7  -1  7  1* 

Similarly,  if  f  only  vary,  £  **  "*  ~  3  £ 
and  therefore  if  f  and  f  both  vary, 


64  KINETIC   THEORY   OF   GASES.  [CHAP. 

and  (A)  becomes 

f  f  k  J  J    \(&£     d$\d  l  (d%     d%\*\      A 

\dxdydz  \(^-  +  -T-}-j\a!Y  —  -XT  {-£-  +  -£-}+*=$   ••  -(B). 
JjJ  (\dz     dx)  dt  2h  \dz      dx)  } 

It  follows  that  «7  -y  (  -=£  +  ^/-  )  must  be  on  average  positive, 
dt  \dz     dx) 

1  /df     d 
and  equal  to  ST  (•?•*•- 


68.     We  must  however  consider  the  value  of  -j-  (-—  -f  -f  } 

dt  \dz     dxj 

more  carefully  thus  : 

Let  P,  P'  be  two  points  on  a  sphere  S  whose  coordinates 
referred  to  the  centre  0  are  respectively  x,  y,  z,  and  oc,  y,  —  z,  so 
that  PP'  is  an  edge  of  the  parallelepiped,  whose  base  is  dxdy 
at  x,  y.  Let  fp,  %p>  be  the  values  of  f  at  P,  P'  respectively. 

Then  if  S  be  the  volume  of  the  sphere,  and  ~  the  mean  value 

of  -f  within  the  sphere, 
dz 


the  integration  dxdydz  being  throughout  the  space  S.     Let  the 

radius  of  8  be  '  a  '  (Art.  61).    Then  we  may  take  -f  as  equal  to 
T.  dz 


Suppose  a  molecule  at  C  whose  velocities  are  a,  @,  7.     The 
change  with  the  time  of  S  -—  due  to  the  motion  of  that  molecule 


~fj>*  The  change  due  to 

the  motion  of  all  external  molecules  is  zero  on  average. 

Let  now  /  be  the  value  of  /  referred  to  P  as  centre.     Then 

%P  =  -^  .     And  let  /'  be  the  value  of  /  referred  to  P'  as 

2/"a 
centre,  so  that  fp-  =    ^-    .     Then  due  to  the  motion  of  the 

molecule  at  C 

d 
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because   on   average  ^-2/=0  &c.     The   minus   sign   is   here 

required  because  /  is  referred  to  P,  and  therefore  -jj-  has  the 

opposite   sign  to  that  which  it  would  have  if  referred  to  C. 

i  .     df        zdf 
That  is,  -f-  =  --  -r  • 
dz        r  dr 

Now  to  form  equation  (B)  we  shall  have  to  multiply  by  ay. 
And  that  will  make  the  terms  in  a,  ft  disappear,  because 
a?cty  =  aftay  =  0,  on  average  and  therefore  on  average 


And  therefore 


the  integration  being  over  so  much  of  the  plane  of  xy  as  is 

within  S,  that  is  -  ay  -=-  -f-  =  aV  =r/.-  -f- ,  -f-  being  the  value 
'  dtdz         '  2fada    da 

df 

of  -4-  when  r  =  a. 
dr 

By  symmetry 

_a  i^  =  a2  .11^ 
tt£  dx  Zf  a  da 

It  follows  from  (A)  and  (B)  that  in  stationary  motion  of 
the  medium  of  material  points  on  average 


,  Case  of  Spheres  with  Finite  Diameters. 

69.  We  now  pass  to  the  case  in  which  our  molecules, 
instead  of  being  material  points  incapable  of  colliding  with  each 
other,  become  equal  elastic  spheres  of  finite  diameter  c.  It  will 
be  shown  that  as  the  result  of  their  collisions  with  each  other, 

—r- ,  which  in  the  system  of  material  points  was  zero,  becomes, 
so  long  at  least  as  f,  (-—}  ,  &c.  retain  the  same  values  as  before, 
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positive,  and  the  motion  not  stationary.  It  is  worth  while  t< 
show  first  in  a  general  way  by  the  aid  of  a  diagram  how  thi 
effect  is  produced  by  the  collisions. 


Let  us  take  the  plane  of  the  paper  for  the  plane  of  xz  and 

suppose  f  alone  to  vary  in  z,  -f-  being  at   a   given  instant 

dz 
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positive  near  a  certain  point.  That  being  the  case,  there  is  at 
the  instant  in  question  near  that  point  a  mean  motion,  or 
stream,  of  the  molecules  in  direction  x  increasing  as  z  increases. 
We  may  suppose  this  instantaneous  stream  to  be  positive  for 
positive  values  of  z,  negative  for  negative  values  of  z.  It  is 
indicated  in  the  diagram  by  the  long  arrows. 

Let  us  consider  two  kinds  of  collisions,  (1)  and  (2),  the 
directions  of  motion  of  either  sphere  before  and  after  collision 
being  shown  by  the  small  arrows.  For  simplicity  these  directions 
before  collision  are  here  drawn  parallel  to  x. 

If  we  examine  collision  (1),  we  see  that  the  molecule  which 
before  collision  was  moving  towards  positive  x  is  deflected 
upwards,  i.e.  towards  the  positive  stream,  and  the  molecule 
which  before  collision  was  moving  towards  negative  x  is  deflected 
towards  the  negative  stream.  In  either  case  the  effect  is  to 
increase  the  stream  pro  tanto. 

If  we  examine  collision  (2)  we  see  that  the  reverse  is  the 
case.  The  effect  of  the  collision  is  to  diminish  the  stream.  But 

now>  z»  being  positive,  there  are  more  collisions  of  the  kind 
dz 

(1)  than  of  the  kind  (2)  per  unit  of  volume  and  time.  And 
therefore  collisions  tend  on  the  whole  to  increase  the  stream. 

For  the  same  reason,  it  will  be  found,  they  make  —j-  positive. 

It  is  true  that  in  a  rare  medium  the  effect  here  indicated  is 
very  small.  But  the  question  is  only  of  its  sign,  not  of  its 
magnitude. 

70.  Since  every  collision  changes  the  direction  of  motion 
of  the  colliding  spheres,  the  expression  in  Art.  66  derived  from 

dt     dt\  d(ay) 
dz     dx)     dt 

contains,  in  addition  to  the  result  calculated  in  that  article, 
a  new  term  derived  from  collisions  which  we  have  now  to 

calculate. 

5—2 
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Let  V  be  the  half  relative  velocity  of  two  molecules,  A,,  /-t,  v 
its  direction  cosines  before  collision,  X',  //,  v  after  collision. 
Let  XF '=  VX)  vV =  Vz.  The  number  of  collisions  per  unit  of 
volume  and  time  for  given  V  is  27rc2/>  V.  Hence,  for  the  change 
with  the  time  of  VXVZ  due  to  collisions  with  given  V,  we  have 


the  bar  denoting  mean  values.     We  have  now   to   find   the 
values  of  XV  and  \v  for  all  collisions,  given   V. 

At  the  instant  of  collision  let  the  centres  of  the  two  colliding 
molecules  be  A,  B\  so  that  AB  is  the  line  of  centres.  Let  BD 
be  the  relative  velocity.  Or,  if  DE-BE,  one  molecule  has 
velocity  DE,  and  the  other  BE,  in  addition  to  the  velocity  of 
their  common  centre  of  inertia,  whatever  that  may  be.  Let 
the  angle  DBA  =  6.  IfAO  =  BO,0  is  the  point  of  contact. 

Let  BX,  BZ  be  the  directions  of  the  axes  of  oc  and  z,  and  let 


the  angle  between  the  planes  DBA  and  DBZ  be  <£,  and  the 
angle  between  the  planes  DBA  and  DBX  be  </>'. 

Then,  if  x,  y,  z  be   the   coordinates   of  A  referred  to  0 
as  origin,  and  if  the  direction  DE  be  that  of  X,  fj,,  v, 


s\  s\  _  \ 

z  —  —  v  2  cos  6  +  £  V 1  —  v2  sin  6  cos  <f>, 
x  =  -  X I  cos  6  +  I  Vl  - X2sin  6 cos  <£' 
The  complete  definition  of  cf>  shall  be  this  : — 


(22). 
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When  the  plane  DBA,  turning  round  DB,  contains  the  axis 
of  2,  $  —  0  for  that  position  of  BA  which  makes  the  least  angle 
with  the  axis  of  positive  z.  And  <£'  has  a  corresponding  defi- 
nition for  the  axis  of  a?  (see  Art.  72,  post). 

71.     If  f0,  £o  be  the  values  of  f  and  f  at  0,  then  at  A 

'—  fo  =  I  (  -  v  cos  6  +  Vl  —  i/2 sin  6  cos  <£) -f 

(23). 


Then  by  the  same  reasoning  as  employed  in  Art.  67  the 
number  per  unit  of  volume  of  pairs  of  molecules  whose  com- 
ponents of  half  relative  velocity  are 


Vy...Vy  +  dVy,  Vz...Vz 
is  at  A 


+0  dVdVdV 


Here  owing  to  the  variation  of  f  in  z,  Vx  is  involved  asym- 
metrically with  respect  to  Vy  and  Vz.  Let  us  suppose  now  V 
constant,  but  the  angles  6  and  <£  to  vary.  The  general  mean 
value  of  VXVZ  is  of  course  zero.  But  the  mean  value  of  VXVZ, 
given  V,  in  the  asymmetrical  system  is  to  be  found  by  integra- 
ting the  above  expression  according  to  6  and  <£,  and  is  not  zero. 
Let  us  denote  the  mean  value  of  VXVZ,  or  \vV\  for  all  values 
of  0  and  <£,  V  being  constant,  by  {VxVz)t  the  mean  value  of  the 
same  function  when  V  also  varies  being  denoted  by  Vx  Vz.  We 
have  then,  treating  f  alone  as  varying, 

(VXVZ)  =  —  (^  2  sin  6  cos  6  dO  ind6  (1  +  4,h'Vx  (f  -  £>)) 

7T  Jo  JO 

in  which  £  —  f  0  has  the  value  given  in  (23),  and  the  general 
means  of  functions  of  Vx  and  Vz,  or  XF  and  vV,  are  taken 
after  integration  according  to  6  and  <f>.  That  is,  writing  VX/V 
for  \  &c. 
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(VxVz}-=  VXVZ  I*" 2  sin  (9 cos  ddO 
Jo 

^  (^  dO  f^  2sin0cos  6  4>h'l  (  -  -pfcos  0 
Jo  A  \ 

4-   ~V  '* 

^ — —  sin  G  cos  <f> 


3         F      <fe" 

because  the  general  mean  value  of  VXVZ,  or  Xz/F2,  is  zero. 

We  have  thus  proved  that,  on  the  average  of  all  colliding 
pairs  of  spheres,  (F^Fz)  is  before  collision  of  the  opposite 

7£  TCf 

sign  to  -^  or  ~- .     This   does  not   depend    on   the   result   of 

collision.    It  is  therefore  true  if,  instead  of  being  elastic  spheres, 
the  molecules  be  centres  of  force,  or  whatever  be  their  form. 

72.  We  have  next  to  calculate  the  mean  value  for  all 
collisions  of  XV,  or  the  value  of  \v  after  collision.  This  will 
depend  on  the  form  of  the  molecule,  and  we  shall  now  deal  only 
with  elastic  spheres.  It  will  be  found  that  XV  is  of  the  same 

.A  (d%  ,  d%\ 

sign  with    -^  +  -r-    • 
\dz     dx) 

We  use  the  same  notation  as  before.  When  the  plane  DBA 
turns  round  DB  until  it  contains  the  axis  of  z,  let  AZ)  AZ)  which 
are  on  opposite  sides  of  DB,  be  the  two  positions  in  which  that 
plane  is  cut  by  the  circle  which  A  describes.  Let  Az  be  nearer 
to  the  positive  axis  of  z  than  Az'.  Then,  for  AZ)  </>  =  0;  for 
AZ,  <£  =  TT.  Similarly,  when  the  plane  turning  round  DB  con- 
tains the  axis  of  x}  Ax,  Ax  are  the  two  positions  of  A  in  that 
plane,  and,  if  Ax  be  nearer  than  Ax  to  the  positive  axis  of  x,  (£'  =  0 
for  Ax,  (/>'  =  TT  for  Axf. 

Let  e  be  the  angle  between  the  plane  of  DB  and  z  and  the 
plane  of  DB  and  x,  so  that  <£'  =  </>  —  e,  <£  =  <£'-{-  e.  It  will  be 
seen  that,  if  \v  is  negative,  cos  e  is  positive ;  if  \v  is  positive, 
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cos  e  is  negative.     In  either  case 
cos  e  = 


we  have  then,  to  find  X'  and  v, 

v  =  -  v  cos  26  +  Vl  —  v*  sin  20  cos  </>, 


(25). 


Whence,  substituting  Vx'  for  X'F,  and  F/  for  z/F,  we  derive 
the  two  symmetrical  systems 


=  -  F2  cos  20  +  \/Fx*  +  TV  sin  20 


cos 


Vx'  =  -  F.COS20- 
F/=-FzCos20- 


V  V 

z 


VF.-  +  vv 

F.FZ 


sin  20  cos  <f> 


,2  +  F; 


—  sin  20  cos  <f>' 


Fa/=-Frcos20  +  \/F;s  +  Fe2  sin  20  cos  0 

The  term  involving  sin  e  will  disappear  in  the  subsequent  inte- 
gration, and  is  omitted.  The  above  systems  hold  for  all  values 
of  Vx  and  Vz. 

73.     We  now  find  that  for  given  F,  taking  into  considera- 
tion only  the  variation  of  f  ,  or  -p  t 

(Fa/F/)  =  -  l^dB  f*d<l>  2  sin  0  cos  0  (1  +  4^/Fa;f^|r0)Fa;/F;. 

7TJ  0  Jo 

F  '    F' 

But  with  the  above  values  of  X',  v',  or  -~  ,  -^-  , 


[  dO  (*dA  2  sin  0  cos  0F/F/  =  0, 
Jo       JQ 

as  is  easily  seen.     Therefore 

(T7F/)  =  -  f  ^0  I"  ^  2  sin  0  cos  0  ^'F.f^fo 

TTJo  Jo 

=  i  4>h'c    ^d6  f'd6  sin  0  cos  0F* 

JO 

I       F« 

K  IF 


cos  0 


sm 
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by  (23).     In  this  expression  we  have  to  substitute  from  (26) 
the  values  of  Vx'  and  F/,  and  integrate  according  to  6  and  </>. 

To   effect  the  integration  according  to  $,  we  reject  odd 

77" 

powers  of  (f>,  and  for  cos2  </>  write  <r .     We  have  then 

( V*Vz)  =  -  P 4/i'c  [:&  r<ty sin 0 cos  6VX 

TT  az         Jo       Jo 


r,     n .  xlv, 


r  F 

1  d%     t    [**      [" 
=  —  -j—  4<hfc  I    dd  I    dd>  sin  0  cos 

TT  a^         J0       Jo 


+  F2  ) 

^  y-  sin  0  cos  </>  v  F^  F/ 


iTT  ./  Y  2  4-  F  * 

-ycos6>+       y     y  sing  cos 

0  =  (-  F,  cos  2(9  +  VF7TT7  sin  20  cos  </>), 
g  =  (-  F.cos  20  -  V^F.  ^n  20  cos 

or  (  Vx'  F/)  =  -  ^  4Ji'c  ?-^p-  .  f^sin  (9  cos2  0  cos2  2(9  d<9 
dz  JQ 

jfc         TT"^2       T^T 
+  g  4/i'c  -^  i  J    sin  6  cos2  (9  sin2  20  d<9 

/7  fc  XT'  2  /  TT  2    .    y  2\          /•  J,r 

-^4^c    '  V"r  i      sin20cos0si 

^  F  2J0 

4Arc  sin2  6  cos  0  sin  2(9  cos  2(9  d<9 


-  ^  4^'c  ^^  (^sin  (9  cos2  (9  cos2  2(9  d<9 
o«  F     Jo 

Jt         y*y*      rfr 

+  S  4^'C  ^-If   1       sin  6  cos2  6>  sin2  2(9  d6 
dz  *  Jo 


-       4/t'c  F^FJ-      sin2  0  cos  6  sin  20  cos  20  dO 


4-  52  4^c  -V-i  sin2  0  cos  0  sin  20  cos  20  dO 

dz 
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fsin2  6  cos  0  sin  20  cos  20  dO 
dz  Jo 


'    V*?        ^2 

F  *F 

..       -..  ?' 

but 


'       F4     ~3.5' 


or 


Hence  the  expression  becomes 

d%     ,    Fa;2Fg2-ll 
* 


3.5.7 


7  r^  7£ 

By  symmetry  if  we  treat  ~  as  varying  instead  of  -j-,  we 
obtain  the  expression 


3.5^  F 

and  therefore,  given  F, 


And  therefore  comparing  this  with  (24) 

+        ......  (27). 


74.  To  find  the  rate  of  change  with  the  time  of  Fa;Fz 
given  F,  we  multiply  the  last  expression  by  27rc2/>F  (F  being 
the  half  relative  velocity),  that  is  by  the  number  of  collisions 
which  take  place  with  given  F  per  unit  of  volume  and  time. 

Let  us  denote  the  change  due  to  collisions  by  z-  . 
The  result  is,  given  F, 
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And  then  we  introduce  the  factor 


7T     7T 

and  integrate  for  all  values  of  F^,  Vy,  Vz  and  so  we  obtain 


75.     This  being  the  value  of  ^-  VXVZ)  what  is  the  value 


If  Z7a,  Uyt  Uz  be  the  components  of  the  common  velocity  of 
the  two  spheres,  a,  /3,  7  the  velocities  of  one,  a!,  f¥t  y  those  of 
the  other, 

oL='Ua+Vx>         a'=Ux-Vx, 

p=Uy+Vy>       /3'=Uy-Vy, 


Therefore         ay  +  ay  =2(UXUZ+VX  Fz). 
Hence  since  Ux,  Uy,  Uz  do  not  vary  by  collision, 

9      /  //X  ^      3      TT     TT 


.(28). 


It  follows  that  on  average  of  all  collisions 
9  /     \     9  TT~  ., 


and  therefore  on  average 

fdg     df\d'\l  (I        1   \  /dg     d£\> 

a  +  %)  dt^  -  IA  ^  (3  +  o)  is  +  9  ......  (GI)- 


76.     We  have  next  to  consider  the  change  in  the  term 

a7^l7r+  j    )  due  to  the  variation  of  a,  or  y  by  collisions. 
7  dt  \dz      dx) 

As   before   let  us  treat  f  only  as  varying.     Let  us  use  the 
construction   of  Art.  68.     Consider   two   spheres   in   collision 
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their  point  of  contact  being  at  G  the  centre  of  the  a  sphere. 
The  z  ordinates  of  their  centres  are 

^  (—  v  cos  9  4-  Vl  —  v2  sin  6  cos  </>), 


and  £  (~  v  cos  #  +  Vl  —  z/2  sin  0  cos  0). 

To  find  the  value  of  ~-  fp  due  to  the  collision,  let  us  refer/ 

to  the  point  P.     Its  value  at  the  point  of  contact  being  /, 
its  value  at  the  centres  of  the  two  spheres  respectively  is 


/+     (-  v  cos  0  +  Vl-i>2  sin  6  cos  <£) 


, 


and  /—-*(—  *>cos#  +  Vl  —  y2  sin  0  cos 

with  corresponding  terms  in  ?-  and  -p  which  will  disappear 
on  integration. 

Now  for  the  first  sphere  a  is  increased  by  the  collision  by 
Vyf  —  Vx,  for  the  second  it  is  diminished  by  the  same  quantity. 
Therefore  as  the  result  of  this  collision 

d£p  =  d  ^  =  ~^fc%(-  "  cos  e  +  ^^'  sin  e  cos 

For  the  minus  sign  here  introduced  see  Art.  68. 

Substitute  for  Vx  —  Vx  its  value  from  (26),  multiply  by 
2  sin  6  cos  6,  and  integrate  according  to  6  and  <f>.  That  gives 


f**^  f* 


sin  (9  cos 


where  F={-  v  cos  0  -f  Vl  —  y2  sin  0  cos  <£}, 

G=(-\2  cos2  0  -    ,  X?/      sin  20  cos  <b\  . 
\  \/l-j^  / 

That  is 


Now  multiply  by  27rc2/>  F,  the  number  of  collisions  given  V 
per  unit  of  volume  and  time.     That  gives 

8fP          16          =r^-  1  df 

» 
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and  ^- 

and    S//^<^ 

16     ,  Tr  T7  1   1  d/  0 
=  -  =-=  irtfpVxVt^j-  -/-  £, 
15  Zfada 

where  the  integrations  are  throughout  the  section  of  the  '  a  ' 
sphere  made  by  the  plane  of  xy,  S  is  the  volume  of  the  a  sphere, 

,  1  df  .    ,  ,  ,  1  df    , 

and  -  -f-  is  the  value  of  -  -£-  when  r  =  a. 
a  da  r  dr 

But 

as  in  Art.  68,  therefore 

d  16 


d  d£         16     _  T.  T.   1  1  d/ 
Similarly          -  -^  ?  =  -  —  Trc?pVxV2  =r>-  ^-  , 
9^cZa;          15  Zfada 

d    d       d\          32  1 

and 


This  expresses  the  change  of  -^  +  -^  due  to  collision  between 

two  spheres.     Let  a,  /3,  7  be  the  velocities  of  one  of  the  two 
colliding  spheres,  a',  ff,  y  those  of  the  other. 

Then    -=—  contains  the  two  terms 
at 


But  a  =^+7*          a^C^-Fa- 

7=^+n      y=^-^, 

as  in  Art.  75. 

Therefore     (ay  +  a  V)  VXVZ  =  2  Vx*  Vf  on  average, 
because  UXUZVXVZ  =  0  on  average. 
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Therefore  for  each  sphere  separately 


8  8lc?/ 

=  -  ^  7rc3p«V  -  4-  > 
15  '  a  da 

because         ay  =  iT,8^*  +  F,fF.f  +  JJ/F.1  +  ^'F/ 
=  4  F/  F,2  on  average. 

We  have  then  on  average,  remembering  Art.  68  ......  (C), 

4 


Adding  this  to  Ci,  we  find  for  the  term  in  -y-  due  to  collisions, 

dt 

3M     2          [([,    ,    ,    1   /dt     dfV 
Si  =  8  ^  JJj***  2S  (df  +  SJ 

We  interpret  this  as  follows.  As  the  result  of  collisions, 
molecules  moving  in  a  given  direction,  as  that  of  a,  are  on 
average  deflected  towards  that  side  where  the  average  motion 
of  the  molecules,  that  is  the  stream,  is  for  the  time  being  in  x. 
And  so  collisions  tend  to  increase  the  stream.  Compare  Art.  (57). 

77.     Now  referring  to  equation  (A),  we  see  that  --^-  ,  which 

in  the  medium  of  material  points  was  zero,  and  so  the  motion 
stationary,  has,  when  the  spheres  have  finite  diameter,  become 
positive,  namely 


dx 


of  which  the  first  two  terms  are  together  zero,  as  in  (A),  and  the 
third  is  positive.     Therefore  if  in  the  system  of  spheres  of  finite 
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diameter  (-^jr)  &c.  have  the  same  values  as  in  the  system  of 

dM  . 
material  points,  —  ^—  is  positive. 

Cut 

But  in  stationary  motion  --=-  must  be  zero.     Therefore  in 

dt 


stationary  motion  of  the  system  of  finite  spheres  (-p)    &c., 

\dz  J 

cannot  have  the  same  mean  values  which  they  have  in  the 
system  of  material  points. 

78.  It  might  be  suggested  that  the  explanation  of  the 
difficulty  consists  in  the  fact  that  when  the  molecules  have 
finite  diameters  the  quantity  of  momentum  transferred  across 
any  plane  per  unit  of  area  and  time  is  increased  in  a  certain 
ratio,  1  :  1  +  £/e,  where,  in  the  case  of  elastic  spheres,  K  =  f  7rc3/o. 
But  this  increases  both  terms  in  (A)  in  the  same  ratio,  and 

therefore,  so   far  as   this  property  is  concerned,  -=—  remains 

zero.  This  effect  arises  from  the  usual  convention  concerning 
elastic  bodies,  according  to  which  there  is  on  every  collision 
an  instantaneous  transfer  of  momentum  through  a  certain  space. 
It  has  nothing  to  do  with  the  change  of  direction  of  the  relative 
velocity,  which  also  occurs  on  collision,  and  which  we  now  have 
considered.  The  first  effect  would  take  place  in  precisely  the 
same  way,  if  after  every  collision  the  direction  of  the  relative 
velocity  were  restored  by  Maxwell's  corps  of  demons  to  what  it 
was  before  collision. 


CHAPTER  VI. 

CORRELATION   OF   VELOCITIES.      ELASTIC   SPHERES. 

79.     THE  solution  of  the  problem  of  the  last  chapter  is 
this.     In  the  system  of  spheres  of  finite  diameter  f  has  become 

,   /df     d£\2  ,       ,  (d%     dt,\*     fd?     d? 

f  +  f  ,   and  (-f-  +  -f-)    has  become  f^  +  ^     +  (-f-  +  -J- 

\dz      dxj  \dz      dxj       \dz       dx 

That  gives  by  arranging  the  terms  in  (D) 


We  may  assume  that  a7j^  +          ^s  no^  altered   by  the 


substitution  of  f  +  f7  for  f,  or  ay  -=-  (-*-  +  -~-  j  =  0. 


The  first  two  terms  of  this  equation  are  together  zero 
by  (A)  Art.  65.  The  remaining  two  terms  can  be  made  zero  by 
suitably  choosing  the  ratio  of  the  means 


. 

\dz      dx         \dz     dx 
We  have,  namely,  writing  K  for  f  -rrc3/?, 

K  [i*s  ,  3?  v  ,  (**L  .  i?  vi  _  (*e  ,  ^  V 

K\\dz+  dx)  +\dz  +  dxl}-\dz  +  dx)  ' 
whence 

rf  +S)'=TM3+gy=«(§+gy 

V^       rfa?/       \-K\dz       dx)         \dz       dxj 
if  K-}  Ac3  &c.  are  negligible. 
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80.     Instead    of    using    Mxz  =  III  dccdydzuy  (  -p  +  -^-j  ,  or 
the  analogous  terms  Mxy  or  Myz,  we  may  make  use  of 


or  the  analogous  terms  Myy  or  Mzz,  and  we  obtain  the  same 
result.  I  have  set  out  the  mathematical  work  in  the  Appendix 
by  reference  to  Arts.  66,  68.  We  may  therefore  now  use  the 
complete  expression  for  M,  namely 


and  we  then  reduce  -=-  to  zero  by  making 
ot 


dj'\*  _     K     /df  y     /^V  _     *     /5f  V  fr 
dx)  ~T^K\dx)  '    \d)  ~T^K\d)  c 


dy 

Now  so  long  as  the  distribution  of  velocities  is  represented 
by  e~AQ,  and  Q  is  the  sum  of  squares  of  the  velocities,  f  2  &c.  and 

(-r^)  &c.  must  have  either  the  same  or  less  values  than  those 
\dz] 

found  by  the  method  of  Art.  63,  for  material  points.  In 
stationary  motion  therefore  in  the  system  of  finite  spheres  Q 
cannot  be  the  sum  of  squares  only.  It  must  contain  products 
of  the  form  b  (OLOL  -f  /3/37  +  77').  And  the  question  now  is,  what 
must  be  the  form  of  the  coefficients  b,  or  as  I  shall  call  them 

coefficients  of  correlation,  in  order  that  the  ratio  I  ~  j    :  I  ~  j 
may  have  the  required  value. 

81.  By  way  of  illustration  let  us  suppose  n  molecules  within 
a  sphere  of  radius  jR,  their  velocities  being  distributed  according 
to  Maxwell's  law  €~h®,  where  Q  is  the  sum  of  squares  of  the 
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velocities,  and  let  us  impart  to  each  molecule,  in  addition  to  the 
velocity  which  it  has  in  that  motion,  the  small  velocity  X'  in  x. 

This  being  done,  the  chance  that  one  of  these  molecules 
shall  have  velocity  represented  by  lines  drawn  from  the  origin 
to  some  point  or  other  within  an  element  of  volume  dfl,  whose 
centre  is  distant  p  from  the  origin  in  direction  making  the  angle 
T|T  with  the  axis  of  x,  is 

Ae-W  (1  +  2hpX'  cos  ^)  dfl. 

The  chance  that  another  molecule  shall  have  velocity  repre- 
sented by  lines  drawn  from  the  origin  to  some  point  or  other 
with  an  element  of  volume  d£lr  distant  q  from  the  origin,  in 
direction  making  the  angle  <f>  with  p,  is 

Ae~h<P  d£l'  (1  +  2hqX'  (cos  ^r  cos  (f>  +  sin  ^  sin  <f>  cos  e)}, 

where  €  is  the  angle  between  the  plane  of  p  and  q  and  the  plane 
of  p  and  x.     The  chance  that  both  molecules  shall  have  the 
velocities  aforesaid  is 
^V-KuH-rf)  dQdn'  (i  +  zhpX'  cos  i/r)  (1  +  ZhqX'  (cos  -f  cos  <£ 

+  sin  i|r  sin  </>  sine)}. 

Multiply  this  by  J  sin  ^  d^r,  and  integrate  for  all  values  of 
e  and  T/T,  and  we  obtain 

(47,2  \ 

1+  —  X'*pqcos<f>), 

expressing  now  the  chance  that  the  velocities  p  and  q,  of  two 
molecules,  in  whatever  direction  either  is,  shall  make  the  angle 
</>  with  each  other.  Now  let  u,  v,  w  be  the  components  of  p ; 
»',  v',  w'  those  of  q.  Then 

f-  +  (f  =  u*  +  v2  +  w2  +  ?e/a  +  v'*  +  wz, 
dfldQ,'  =  dudvdwdu'dv'dw', 
pq  cos  </>  =  uu'  +  vv'  +  ww'. 

Therefore,  since,  with  these  substitutions,  the  expression  does 
not  contain  <f>,  it  gives  the  chance  that  two  molecules,  both 
within  the  R  sphere,  shall  respectively  have  velocities  u...u+dut 
&c.,  in  the  form 

+  ^-  X''2  (UU   -f  VV'  +  WW')  \  , 

6  J 
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that  is,  assuming  X'2,  &c.,  to  be  small, 

A*e-**du...dw', 
in  which 

4A 

$  =  u2  +  v2  +  w2  +  u'2  +  v'2  +  w'2  -  -s-  X'2  (MM'  +  tw'  +  W). 

o 

Evidently,  if  we  take  F'  and  Z'  into  account,  we  shall  have  to 
substitute  X'2  +  Y'2  +  Z'2  for  X'2  in  the  last  expression. 

We  see  then  that  in  the  case  now  treated  the  coefficient  of 
correlation  b  has  the  form 


or,  since  Z~2+  F2"  +  Z*  = 

X,  F,  Z  denoting  the  velocities  of  the  centre  of  inertia  of  the  n 
molecules, 

b  =  -  -  (X'2  +  F2  +  Z'^KX2  +T2  +  Z~2)   ......  (31). 

82.  Next  let  us  employ  the  converse  method,  namely,  the 
coefficients  b  being  supposed  given,  to  find  the  ratio  f'2:f2. 
And  first  let  us  consider  a  finite  number  of  molecules,  and  all 
the  b  coefficients  equal,  and  very  small,  so  that  62,  b3,  &c  ,  may 
be  neglected. 

We  have  generally,  referring  to  the  determinant  of  Art.  63, 
f  2  =  n-,  ,  and  therefore 


-  fca  ,  3       * 
d' 


3f  2  being  the  variation  of  f  2  due  to  the  introduction  of  the  b 
coefficients.  Now  it  is  proved  in  the  Appendix  (n)  that  in  this 
case  9D  =  0.  Therefore 
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Also     A,-T+7?,        /,/,,        /,/„... 


/l/3>  /2/3>  1    +/32>   ". 

and  (see  Art.  83)  9  (1  +/>»)  =  -  6/lf     3(1  +/,»)  =  _  6/2>  &c., 


If  we  were  now  to  make  every  /=  1,  we  find  that  every 
constituent  in  Dn  as  it  originally  stood  is  in  the  varied  form  of 

Ai  multiplied  by  .     Hence 


But  on  the  same  assumption  that  every  f=  1, 


_ 

2nh' 

And  2  =  a 


n 


whence 


ri  being  large,  or 


as  we  obtained  in  a  different  way  in  Art.  81. 

83.  Next,  let  us  assume  6,  the  coefficient  of  correlation 
between  the  velocities  of  two  molecules,  to  be  a  function  of  the 
distance  r  between  the  molecules  to  which  it  relates.  And  let 
it  be  required  to  find,  on  this  hypothesis,  the  ratio 


For  the  system  of  finite  spheres  we  shall  have 

6—2 
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where  bpq  is  the  value  of  b  for  the  two  molecules  whose  velocities 
are  ap,  aq,  &c.,  whence,  substituting  for  an  as  in  Art.  63, 

Q  =  (1  +/<»  -  b^fj  ttl2  +  (1  +/22  -  bmf,)  «22  +  &C. 


C.     .  .  .(32). 

Then  + 


or 


every  6  being  small,  so  that  its  square  and  higher  powers  may 
be  neglected,  and,  for  the  same  reason,  the  variation  of  D  may 
be  neglected,  as  will  be  easily  found.     Appendix  (n). 
Now  in  this  case  [see  Appendix  (o)] 


in  which  p  and  q  respectively  take  all  values  from  1  to  n,  there 
being  n  molecules  in  the  system.  Also  if  we  assume  the 
molecules  to  be  scattered  through  space  without  regard  to  the 
position  of  p,  with  density  p,  ^bpqfpfq  may  be  represented  by 
the  integral 


and  if  bpq  is  very  small  except  for  very  small  values  of  rpq  (see 
Appendix  (p)),  we  find 


&^r   ......  (33), 

••  c 

,  if  B  =  f 


and  therefore 


because  2/p8  =  Du  ;  and  if  6  be  so  chosen  that  ^b  =  —  K, 

¥*  =  *?    ........................  (34). 

84.     Again,  it  can  be  shown  that  under  the  same  conditions 

(35), 


and  this  is  the  required  ratio. 
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For  since  f  =  ax  g,+  ou,  ^.+  &c., 

d%          d  f,  d  /2 

^   =   ^l    J~    ^>+   «2  :rx6- 

cfo          dzzf         dz  zf 


if  *-**- 

1          <2f  ft  ft 

whence  an  =  ^-  .  -^  -  ^  ^-  -  02  ^-  -  &c., 

t/7l       ttZ'  17  7j,  un 

and  Q  =  a2  +  a2-l-  ...  +  a2 


Now  the  determinant  of  this  expression  differs  from  that  of 
Art.  63  only  in  the  form  of  the  coefficients.     Therefore 

f*fy  -  — ii  -  JL  2#2 

7          J         "^    «Tk  1     T\    ~™ "    iTfc  T       ^*vr     • 

air/ 


Again,  when   the   b   coefficients  are  introduced,  we  have, 
substituting  for  CL,I  in 


and  by  the  same  process  as  before 


" 
It  follows  from  the  above  that  any  form  of  6  which  makes 

=  -K,  or  £  I    4>7rr*pbdr  =  —  /c, 

^0 
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and  satisfies  the  equation  (33),  gives  a  solution  of 


d-z  dz    ' 

and  the  solution  is  independent  of  the  form  of  the  function/, 
so  long  as  that  function  satisfies  the  conditions  of  Art.  61. 

The  value  of  26  is  thus  determinate.     But  the  particular 

form  of  b  which  will  make  J  /    4nrr2pbdr  =  —  K  is  not  determi- 

J  c 

nate.     For  instance,  let 

/*-3c* 
~~       ' 


That  is  a  solution  of  (33)  and  (35)  for  all  values  of  //,  greater 
than  3. 

A      •     i  ,        ,         /*  —  3  c3     , 

Again,  let        b  =  —  c        —  when  r  <  a, 

p  -  3  c3  a* 

o=  --  -  when  r  >  a. 

fju     a?  r* 

That  also  is  a  solution  for  all  values  of  fj,  greater  than  3. 
Here  a  may  be  the  radius  of  a  sphere  which  on  average  contains 
one  molecule,  as  in  Art.  61. 

85.     The  law  of  distribution  of  velocities  being  €~hQ,  where 

Q  =  C^tti2  -f  &ia«iff2  +  a.2«22  +  &C., 

let  a,  /3,  7  be  the  velocities  of  a  molecule  of  mass  m,  a,  ft  ,  y' 
those  of  any  other  molecule.  Then,  b  being  small,  it  is  easy 
to  prove  that 


that  is,  ^m2  (OLOL  +  ftp'  +  yy)  is  equal  to  the  Virial  (taken  as 
positive)  of  the  intermolecular  forces,  represented  in  case  of 
elastic  spheres  by  collisions.  In  Art.  58,  dealing  with  finite 
intermolecular  forces,  we  advanced  a  step  further,  proving  name- 
ly that  the  value  of  \m  (cu7  +  ft  '{3'  +  77')  for  any  given  r  is  equal 
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to  a  particular  part  of  the  Virial  specially  related  to  that  value 
of  r.  If  we  are  to  make  a  similar  step  in  advance  in  the  case  of 
elastic  spheres,  it  must  be  by  dividing  the  Virial,  if  that  be 
possible,  into  distinct  portions,  each  specially  related  to  a 
particular  r.  Instead  of  the  Virial  of  the  definite  force  which 
ra'  exerts  on  m  at  distance  r,  let  us  calculate  the  chance  that 
ra',  so  situated,  shall  be  the  next  sphere  to  collide  with  m, 
multiplied  by  the  mean  Virial  of  that  collision  if  it  occurs. 
That  gives  us  a  part  of  the  whole  Virial  which  stands  in 
special  relation  to  r. 

Consider  a  molecule  m.  The  number  of  collisions  between 
m  and  m  with  relative  velocity  V  in  time  dt  is  proportional  to 
Trc-pVdt.  Suppose  the  whole  system  to  have  attributed  to  it  a 
velocity  equal  and  opposite  to  that  of  m  for  the  time  being,  so 
as  to  reduce  m  to  rest.  Then  the  sphere  m'  which  collides  with 
m  with  relative  velocity  V  must  have  described  a  certain  free 
path,  which  must  have  commenced  at  some  distance  r  ...  r  +  dr 
from  m. 

The  chance  that  it  shall  have  commenced  at  distance 
r...r  +  dr  is  of  the  form  q€~vdr,  where  q  is  independent 
of  r,  but  a  function  of  V.  But  m'  starting  from  collision  at 
distance  r  from  m,  is  as  likely  to  be  moving  in  any  direction  as 
in  any  other;  we  must  therefore  divide  e~^  by  47rr2.  Again, 

2      Vc 

the  Virial  of  the  collision  if  it  takes  place  is  =  m  —  -  by  Art.  28. 

o       4 

f* 
Let  us  then  introducing  the  factor  —  ,  assume  for  given  V 

7TC2FC  2          VC 

- 


and  therefore,  given  V, 

I    4>7rr*pbdr  =  -  7rC*p  ^7—  I     — 
2          mV* 


> 
and  taking  means  for  all  values  of  V 
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To  express  b  as  a  function  of  r,  we  must  perform  the 
integration  according  to  V,  and  not  according  to  r.  That  is 

1,2  _9TC*_ 

23  47rr4 
in  which  q  is  a  function  of  V,  and  C  is  the  usual  constant. 

That  appears  to  satisfy  all  necessary  conditions  for  b.  But 
it  is  not  a  complete  solution  of  the  problem,  because  the  intro- 
duction of  the  factor  — -  is  not  proved  to  be  necessary,  except 

to  satisfy  (33).  It  is  possible  that,  as  in  Art.  57,  the  corre- 
lation, that  is,  6,  may  depend  on  the  angles  which  the  velocities 
of  two  molecules  make  with  r,  the  line  joining  their  positions, 
as  well  as  on  r. 


CHAPTER  VII. 

FINITE   INTERMOLECULAR   FORCES.      BINARY   ENCOUNTERS. 

86.  IT  is  worth  while  to  apply  the  method  of  Chapter  v.  to 
molecules  between  which  finite  forces  act.     Nothing  is  assumed 
concerning  the  force  acting  between  two  molecules,  except  (1) 
that  it  acts  in  the  line  joining  their  centres,  (2)  is  a  continuous 
function  of  the  distance  between  their  centres,  (3)   becomes 
repulsive  and  infinite  as  that  distance  is  indefinitely  diminished, 
(4)  is  evanescent  at  all  distances  greater  than  a  certain  line 
c  which  is  very  small  compared  with  the  dimensions  of  the 
space  in  which  the  system  of  molecules  is  moving.    A  sphere  of 
radius  c  described  about  a  molecule  shall  be  called  its  sphere  of 
action. 

It  seems  necessary  to  assume  (3)  that  for  sufficiently  small 
distances  the  force  becomes  repulsive  and  infinite,  if  we  are  to 
hold  that  there  is  a  very  small  distance  within  which  two 
molecules  cannot  approach  each  other. 

87.  With  this  definition  of  a  molecule,  let  us  assume  that 
the  number  of  molecules  scattered  through  a  given  space  is  so 
small,  or  the  aggregate  volume  of  their  spheres  of  action  bears 
so  small  a  ratio  to  the  space,  that  no  molecule  is  ever  within 
the  sphere  of  action  of  more  than  one  other  molecule  at  the 
same  instant.     In  other  words,  the  encounters  are  binary.     Or 
that  the  number  of  such  complex  encounters  which  take  place 
in   unit   of  volume   and  time  bears  so  small  a  ratio  to  the 
number  of  binary  encounters,  that  we  may  without  appreciable 
error  assume  all  the  encounters  to  be  binary. 
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88.  Suppose  first  a  molecule  P  fixed  in  space.  Let  PX 
be  any  line  through  P. 

Let  Q  be  a  point  just  outside  of  the  sphere  of  action  of  P 
and  distant  I  from  PX.  And  let  another  particle  at  Q  of  mass 
m  move  with  velocity  q  parallel  to  XP.  An  encounter  takes 


place,  and  since  the  encounters  are  assumed  to  be  binary, 
m  describes  a  symmetrical  curved  path  QRQ',  R  being  the 
apse.  Let  the  angle  XPR  =  ^.  The  whole  effect  of  the 
encounter,  so  far  as  regards  the  change  of  direction  of  the 
motion  of  m,  is  the  same  as  if  m  received  an  impulse  2mq  cos  -^ 
in  direction  PR.  Then  also  P  experiences  an  impulse  2mqcosifr 
in  direction  RP. 

Let  the  number  per  unit  of  volume  of  molecules  which  at 
or  near  Q  have  velocity  q  ...  q  +dq  in  direction  XP  be/(</)  dq. 
Then  the  sum  of  the  impulses  which  P  experiences  in  unit  of 
time  from  all  the  molecules  which  so  move  is  2mf(q)  q*dq  cos  i/r, 
and  the  mean  force  per  unit  of  time  on  P  due  to  all  molecules 
moving  from  Q  in  direction  X P  is 


2m  lf(q)  qz  cos  ^r  dq, 


o/r  being  a  function  of  q,  and  the  integration  including  all  values 
of  q.     If  I  also  varies  ty  is  a  function  of  q  and  I. 

Now  suppose  a  circle  described  about  X  with  radius  I  in  a 
plane  at  right  angles  to  XP.  Let  m,  moving  parallel  to  XP, 
be  before  encounter  anywhere  on  that  circle,  and  let  the  angle 
made  by  the  plane  QXP  with  the  plane  of  the  paper  be  e.  R 
will  evidently  lie  on  a  circle  parallel  to  the  circle  described 
about  X. 
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If  PZ  be  drawn,  at  right  angles  to  PX  in  the  plane  of  the 
paper,  the  mean  force  on  P  resolved  in  direction  PZ  due  to  the 
molecules  moving  from  Q  in  direction  parallel  to  XP  is 

2m  cos  e  I  f(q}  (f  cos  -vfr  sin 


I  /(<?)  <f 


And  if  now  I  also  vary,  the  mean  force  on  P  in  direction  PZ 
due  to  all  molecules  moving  in  direction  XP,  whatever  the 
value  of  I,  is 

2m  I  cos  ede  I   Zirldl  \  f(q)  q2  cos  T|T  sin  ^frdq, 

.'  o  J  Q  J 

i/r  being  now  a  function  of  I  and  q.     It  is  of  course  zero. 

89.  If  P,  instead  of  being  at  rest,  be  moving  with  velocity 
u  in  direction  PX,  we  must  write  u  +  q  for  q  in  the  above 
expressions.  The  force  on  P  for  any  given  value  of  e  will 
then  be  in  direction  RP,  and  equal  to 


2m     ^wldl     dqf(q)  (q  +  w)2  cos  ^. 
And  the  mean  force  on  P  in  direction  PZ  will  be 

2m  I  cos  ede  \    Inldl  lf(q)  (q  +  u)2  cos  \/r  sin  ^rdq. 

J  o  Jo  J 

Now  let  the  plane  of  the  paper  be  the  plane  of  XZ,  XP  the 
direction  of  x.     Let  the  function  f  have  at  P  the  value  f0,  and 

at  any  point  whose  z  ordinate  is  dz  the  value  f0  -f-  -^  9#.     Then 

Ct<2 

the  velocity  of  P  relative  to  Q  will  be  q  +  u  +  ~  £  cos  e. 


And  the  mean  force  on  P  in  direction  PJ£  is  now 


cos  ir  sin 


4t7rm  I   cos  ede  I  Idl  \  dqf(q)  (q  +  u  -\-  I  ~-  cos  e  J 

That  is,  if  u,  v,  w  be  the  component  velocities  of  P, 
-y.  =  —  4-Trm  I  cos  ede  I  Idl  I  dqf(q)  I  q  +  u  +  I  -p  cos  e  J  cosi/rsin>|r 

=  —  4-Trm  I   Wi  I  dqf(q)  (q  +  u)  l-^ 


cos  \r  sn 
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And  u  being  negative  let  us  write  —  u  for  u.     Then 

u  -=-  =  +  4?rm  I   Idl  I  dqf(q)  (qu  +  it2)  £  |5  cos  i|r  sin  -»|r 
etc  J  o      J  ctz 

=  +  4-Trm  I  £2e^  /  dqf(q)  uz  cos  i|r  sin  ^  ~  , 


~ 


on  average  of  all  values  of  ^  ;  or  writing  ^ry  for 


where  «;  =  4-7T  I  I2dl  I  dqf(q) 


cos      sn 


By  symmetry  w-j-  has  the  same  value  mutatis  mutandis, 
at 

and  therefore  -y  uw  - 


and 


90.  This  agrees  in  form  with  the  result  obtained  for 
elastic  spheres  in  Art.  76.  But  K,  instead  of  being  a  known 
constant,  now  depends  on  the  unknown  law  of  force  and  is  a 
function  of  (f.  Further,  in  the  expression  for  K,  COST/T  is 

necessarily  positive,  because  -fy  lies  between  zero  and  —  ,  and 
also  cos  T|T  sin  -fy  is  necessarily  positive,  and  K  is  necessarily 
positive.  Now  if  ty  >  -r  ,  cos  ty  sin  -fy  diminishes  as  -\Jr  increases ; 

d^l/*  IT 

that  is  as  q  increases,  because  -f-  is  positive  ;  if  i/r  <  —  ,  cos  -\Jr 

dq  4 

sin  T|T  increases  as  q  increases. 

Also  the  condition  that  -\/r  shall  be  less  than  —  is  that  I  be 

4 

less  than  a  certain  magnitude,  which  we  may  call  10,  and  the 
greater  q  is  the  smaller  must  10  be.  It  follows  that  as  q2,  or 
the  mean  kinetic  energy,  increases,  it  becomes  less  probable 

that  >/r  shall  be  less  than  —  ,  and  therefore  more  probable  that 
cos  -\Jr  sin  i|r  shall  diminish  as  q  increases. 
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Therefore  if  c£  be  great  enough,  cos  i|r  sin  ^  taken  on  average 
of  all  values  of  I  and  q,  will  ultimately  diminish  as  q2  increases. 
And  therefore  K  must  for  sufficiently  high  values  of  (f,  that  is 
of  the  mean  kinetic  energy,  or  temperature,  if  the  two  are 
identical,  either  diminish  as  the  temperature  increases,  or  in- 
crease less  rapidly  than  the  temperature. 

dM 
91.     We   have  now  to  consider  the   other   term   in 

at 

namely, 


Let  us  again  assume  the  encounters  to  be  binary.  As  in 
Art.  68  describe  a  sphere  of  radius  a  about  any  point  0,  and  let 
PP'  be  a  double  ordinate  parallel  to  z,  P  and  P'  being  on  the 
surface  of  the  sphere. 


Then 


Suppose  two  molecules  to  encounter  one  another,  the  point 
of  contact  of  their  spheres  of  action  being  at  the  centre  of  the 
sphere.  Let  VX)  Vy,  Vz  be  the  components  of  their  relative 
velocity  before  encounter.  We  may  suppose  their  common  centre 
of  inertia  at  rest,  since  its  motion  will  not  affect  f  P  on  average. 

Let  us  consider  the  change  of  -^  due  to  this  encounter 
only,  that  is 

4?ra3 


First  to  find  -^ ,  let  /  be  referred  to  P  as  centre,  and  let 

/0  denote  its  value  at  the  centre  of  inertia  of  the  two  spheres  of 
action.  Let  /  be  the  value  of  f  at  the  centre  of  either  sphere 
when  the  encounter  begins.  After  the  encounter  has  ceased  / 
will  have  become/',  corresponding  to  the  final  relative  position 
of  the  two  spheres  of  action.  Let  V.x,  Vz  have  become  Vx',  Vz. 
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If  d%P  be  the  change  in  fp  due  to  the  encounter,  3fP  is, 
itting  the  factor  =-*  which  is  o 
tional  to,  and  has  the  same  sign  as 


omitting  the  factor  =-*  which  is  on  average  constant,  propor- 


that is,  if  we  write  VX  =  \V,  Vxf  =  X'F,  &c.,  the  same  sign  as 

x/'-V- 

Now  let  0  be  the  angle  between  the  relative  velocity,  whose 
direction  cosines  are  X,  //.,  v,  and  the  line  of  centres  at  the 
commencement  of  encounter. 

Then  2i|r  —  6  is  the  corresponding  angle  at  the  end  of  the 
encounter,  ty  being  the  angle  so  denoted  in  Art.  88. 

Q 
Then  we  have,  ^  being  the  radius  of  the  sphere  of  action, 


/=/0  +  -  (-  v  cos  6  +  sin  0  cos  0)      , 

/'=/0  +  I  (-  v  cos  2«f  -  6  +  Vl-^sin  2f  -  <9  cos  </>)  -f  ; 

A  (Hz 

also  as  shown  in  Art.  72, 

X'  =  —  \  cos  2i|r  --  sin  2i/r  cos  0, 

vl  —  v2 

and  therefore,  omitting  cos  </>  and  writing  \TT  for  cos2  <^>  as  the 
result  of  integrating  according  to  </>, 


\'f  =  X'/0  +    Xz^  .  (cos  2-»/r  cos  2i|r-0  -  \  sin  2-f  sin 


•  cos 


—  cos  6  2  cos2  -\Jr  sin2  ^ 

+  sin  0  2  sin  i|r  cos  i|r  (2  cos2  ^  —  1) 

+  sin  0  sin  ^  cos  i|r  (2  cos2  -v^  —  1)}. 
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Also         \f=     \f0 


+  \     (-  v  cos  0  +  Vr^  sin  6  cos  <£) 

C  /i^/* 


on  average.     Also  for  the  pair  of  spheres  (X'  —  X)/0  =  0,  and 

r%  cos  e  sin  eae  (x'y  -  x/) 

Jo 

c  df  f  i* 
=  \v-^-\     2  cos20  sin  6  d0  (6  cos4  -^  -  6  cos2  i/r  +  2) 

2  Gk£  J  o 


which  can  be  easily  shown  to  be  positive,  6  being  less  than  i/r, 

and  >/r  less  than  ATT.     Therefore  37  -r-  (  ~  +  ^  )  is  positive. 

1  rf^  \c^^      dx) 

92.  If  instead  of  being  repulsive  as  hitherto  assumed,  the 
force  between  the  molecules  P  and  m  be  attractive,  it  is  possible 
that  the  approach  of  the  two  may  result  in  a  permanent  union, 
each  molecule  revolving  in  a  closed  curve  round  the  common 
centre  of  inertia.  Leaving  this  case  out  of  consideration,  the 
form  of  the  path  described  by  m  with  reference  to  P  regarded 
as  fixed  will  for  attractive  force  be  as  in  this  figure. 
Q  _ 


The  angle  ^  is  now  between  ^  and  TT  instead  of  between 

zero  and  ^  as  in  the  case  of  the  repulsive  force.     The  impulse 

on  m  is  represented  as  before  by  2mq  cos  i/r,  cos  t/r  being  now 
negative.  The  mean  force  on  P  due  to  a  succession  of  en- 
counters is  in  direction  PR  instead  of  RP.  The  reasoning  of 
Art.  89  will  apply  equally  well  to  this  case,  except  that  in  the 
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result  u-jr  will  be  of  the  opposite  sign  to  -^-  instead  of  the 

, (d£     d£\  d 
same  sign,  and  I  -^-  +  -r-  )  -r-uw  will  be  negative. 

\d/z     cLx  /  (tt 

That  leads  in  stationary  motion  to  a  positive  sign  for  the 
coefficients  6.  It  is  conceivable  that  the  force  between  two 
molecules,  being  a  function  of  the  distance  r  between  them, 
may  change  sign  at  some  value  of  r.  If  however  there  be  a 
limit  of  distance  within  which  the  centres  of  inertia  of  two 
molecules  cannot  approach  each  other,  the  force  must  ulti- 
mately, as  r  is  indefinitely  diminished,  become  repulsive.  And 
as  the  density  increases  the  near  approaches  becoming  more 
frequent,  the  repulsive  forces  must  ultimately  predominate. 
The  effect  of  the  attractive  forces,  if  such  exist,  will  be,  not  to 
make  the  6's  positive,  but  to  make  them  have  less  negative 
value  than  they  otherwise  would  have. 


93.  So  far  we  have  assumed  the  encounters  to  be  binary, 
so  that  for  each  encounter  either  molecule  describes  a  complete 
symmetrical  orbit  with  reference  to  the  other,  complete  that  is 
in  the  sense  that  after  the  two  have  separated  to  the  very  small 
distance  c  from  each  other,  the  motion  of  either  is  sensibly 
rectilinear  in  the  asymptote.  If  we  remove  the  restriction  to 
binary  encounters,  we  may  reason  as  follows. 

Assume  -^  to  be  positive  at  the  origin  P.   Then  a  molecule 

at  P  moving  with  u  negative  has  greater  velocity  relative  to 
molecules  on  the  side  of  positive  z  than  to  those  on  the  side  of 
negative  z.  On  the  average  therefore  the  potential  of  the 
mutual  action  of  P  and  the  other  molecules  will  be  greater 
with  z  positive  than  with  z  negative.  Let  p  be  that  potential. 

Then  jg  is  positive. 

Therefore   if  u,  v,  w   be  the   component  velocities   of  P, 

dw         dp  .  ..  f          dw  . 

-=-  =  — fr  is  on  average  negative.  And  therefore  u  - -  is  on 
dt 
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average  positive,  or  has  the  same  sign  with  J5  .     Also  its  sign 

cLz 

70 

has  no  relation  to  -j-  . 
dx 

Similarly  w  -j-  is  on  average  of  the  same  sign  as  -=^  ,  and 
independent  of  the  sign  of  -~  .     And  therefore 

d  —       d       d\      dw        du 


is  necessarily  positive. 

And  therefore  in  a  system  of  molecules  of  the  kind  now 

considered  —  ^—  is  positive,  and  the  motion  not  stationary  if 

f2,  i?2,  f2  of  Art.  61  have  the  same  mean  values  as  they  would 
have  in  the  medium  of  material  points.  To  effect  stationary 
motion,  £2,  rf,  f2  must  be  increased  as  we  found  to  be  the  case 
with  elastic  spheres. 


CHAPTER  VIII. 

GENERAL   THEORY    OF   THE   STATIONARY   MOTION. 

94.  HAVING  established  the  results  contained  in  Chapters 
IV.  and  v.,  I  shall  now  assume  for  the  law  of  distribution  of  the 
coordinates  and  velocities  among  an  infinite  number  of  mole- 
cules in  stationary  motion  the  following,  namely :  The  chance 
that  at  any  instant  the  coordinates  of  the  molecules,  n  in 
number,  shall  lie  respectively  between  the  limits 

a?    and 


z\     „     z 
xn  and  xv 


and  their  component  velocities  between  the  limits 
«!  and  «!  4-  da^ 

A    „    &  +  <%, 


and 


s 


(-}  *JDA  e-*«  dx,  .  .  .  dznda,  ...djn. 
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Here  Q  =  a,  (af  +  &2  +  7i2)  +  a*  (***  +  A2  +  722)  +  &c. 

+  &c. 


The  '  a  '  coefficients  are  independent  of  the  positions  of  the 
molecules.  Every  6,  as  bpg,  is  a  function  of  the  distance,  r^,  at 
the  instant  in  question  between  the  molecules  to  whose 
velocities  the  suffixes  pq  relate,  such  functions  satisfying  the 
conditions  of  Chap.  v.  D  is  the  determinant  of  the  coefficients 
inQ. 

Again,  A  may  be  explicitly  a  function  of  the  coordinates 
x,  y,  z,  &c.,  or  may  be  a  constant,  >J~D  is  a  function  of  the 
coordinates,  only  as  contained  in  the  coefficients  b.  Then 
evidently 


Instead  of  using  the  word  chance,  we  may  say  that  the 
above  expression  represents  the  time  during  which  on  the 
average  of  any  very  long  time  the  coordinates  and  component 
velocities  will  be  found  respectively  within  the  limits  aforesaid. 

95.  The  function  Q  may  conveniently  be  put  in  the  form 

Q  =  a^2  +  ^uftiOa  +  i&isaiOj  +  &C. 

+  a^2  +  i^aaCCx  +  JfraC^as  +  &c. 

+  &C., 

in  which  each  line  is  appropriated  to  a  particular  molecule. 

96.  According   to   the   results   obtained  in   Chapters   IV. 
—  VI.  it  is  necessary  for  stationary  motion  that   f  /2  +  ?/2  +  f  /2 
(Art.  79)  shall  have  a  certain  value.     That  necessary  value  it 
will  have  if  suitable  values  be  given  to  the  coefficients  b  as 
functions  of  r.     The  coefficients  b  being  so  determined,  the 
motion   is   stationary,  so   far  as  not  to  be  disturbed  by  the 
collisions  or  mutual  forces  of  the  molecules.     I  assume  now 
that  in  Q  the  b  coefficients  have  the  values  so  determined. 

7—2 
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97.  We  can  now  deduce  the  fundamental  property  of  this 
distribution,  namely,  that  for  any  given  set  of  positions  of  the 
molecules,  or  as  we  shall  call  it,  for  any  given  configuration, 

dQ         dQ 


For  --^  =  20^!  +  612a2  +  613a3  +  &c. 


If,  with  «!  constant,  «2,  «3  &c.,  assume  all  possible  values 
consistent  with  the  conservation  of  energy,  we  find  (see 
Appendix  (/&))  that  the  mean  value  of 

-^  ,  or  2o1«1  +  612«2  +  &c., 

is  -=-  «!,  D  being  the  determinant  of  the  coefficients  a,  b,  &c. 
*+a. 

in  Q,  and  Dn  its  coaxial  minor  obtained  by  striking  out  the  row 
and  column  containing  alf     Therefore  also 

dQ       D 

«!  -j—  =  -=r-  otj    on  average. 
dot}     JJn 

But  it  is  also  shown  (Appendix  (g))  that 


whence  on  average 

~dQ      D 


Similarly          a2       =  -  &c. 


j  dQ          dQ 

and  «^  =  « 


98.  Inasjmj£b_a^the_system  is  infinite,  and  there  is  no  direct 
action  between  any  molecule  and  aloyr"oTirleT~nrDiecule  except  at 
distances  very  small  compared  with  the~~drHiensions  of  the 
system,  not  only  is  the  sum  of  the  potential  and  kinetic 
energies  constant,  but  the  kinetic  energy  of  the  whole  system 
is  itself  separately  constant.  Let  the  kinetic  energy  be 
T=  m  a2  + 
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Here  if  ^  be  velocity  of  translation  m  is  mass,  if  a  be 
velocity/of  rotation  m  denotes  moment  of  inertia,  and  so  on. 
Then/we  have  by  the  constancy  of  T, 

TTljOfj  ~~  +  m^Lz  -T-  +  &C.  =0    ...............  II. 

For  the  same  reason  Q  is  separately  constant,  or 


+  &c=0 
da.1  dt      dc^  dt 

Now  II.  and  III.  suggest  the  inference 
dQ 


An  obvious  solution  of  this  is  Q  =  £2ra  (a2  +  f&  +  T2).  But, 
as  proved  in  Chapters  iv.  and  v.7this  does  not  give  stationary 
motion  when  the  molecules  have  finite  dimensions,  except  in 
the  limiting  case  of  infinite  rarity. 

99.     Concerning  the  Maxwell-  Boltzmann  Law  ^/ 
ra^!2  =  tn&?  =  &c. 

This  follows  at  once  from  Art-  97,  and  IV.f  if  IV.  be  accepted. 
It  must  however  be  noted  that  IV.,  although  a  consistent,  is  not 
a  necessary,  consequence  of  II.  and  III.  For  while  ttj,  Oj,  &c.  are 
unchanged,  let  all  the  molecules  undergo  small  displacements. 

Then  -y-1  ,  -^  &c.,  depending  as  they  do  on  the  intermolecular 
dt      dt 

forces,  will  receive  variations,  which  may  be  arbitrary,  as  the 
displacements    are   arbitrary.     If    under   these   circumstances 

-j-  ,  &c.  remained  unchanged,  IV.  would  be  a  necessary  conse- 

quence of  II.  and  III.     But  -=-  &c.  do  not  remain  unchanged. 

So  IV.  is  not  proved  by  II.  and  III. 

If  we  assume  condition  A,  we  assume  in  effect  that  the 
b  coefficients  are  all  zero,  and  Q  =  Swia2.  It  follows  then  at 
once  from  the  law  of  distribution  e~AQ  that  m^a*  =  w2a22  = 
The  law  is  therefore  true  whenever  condition  A  can  legitimately 
be  assumed.  It  is  true,  that  is,  in  the  limiting  case  of  infinite 
rarity. 
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/    I  am  not  aware  that  it  has  ever  been  proved  in  any  other 
^  case,  or  by  any  assumption  not  equivalent  to  condition  A. 

If  it  were  true  in  the  general  case,  we  should  have,  comparing 


II.  with  III, 


and  therefore  by  Art.  97, 

D  D 


or  m^  Dn  =  mJ)Z2  =  &c. 

But  Du  Z)22,  &c.  are  functions,  not  only  of  the  masses  of  the 
molecules  ml  ra2,  &c.,  but  (if  they  be  elastic  spheres)  also  of 
their  diameters,  (if  they  be  centres  of  force)  of  their  effective 
volumes.  It  seems  therefore  to  follow  that  the  law 

cannot   hold   universally.      It    can   be   accepted   only   on    the 
authority  of  the  great  physicists  by  whose  name  it  is  known. 

100.     If  therefore  we  write 

dQ  / 


=yLt202,  &C., 


fa,  fa,  &c.  will  in  general  have  different  values  for  different 
kinds  of  molecules  in  our  system.     Only  in  the  case  where  the 

molecules  are  all  of  the  same  kind,  may  we  write  ~=  \  m  =  LLO., 

da  I 

where  //-  is  constant. 

In  all  cases  however  in  which  the  b  coefficients  are  very 
small  fa,  fa,  &c.  will  be  very  approximately  equal  to  each  other. 
For  in  the  limiting  case  when  the  6's  are  all  zero 


and  we  know  that  ax  =  ra^  a2  =  m2,  &c.  When  the  b's  are  not 
zero,  every  coaxial  minor,  as  Dn,  consists  of  the  product  of  all 
the  axial  constituents  m2m3 . . . ,  and  of  other  terms  each  of  which 
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contains  the  product  of  not  less  than  two  6's.  The  quantities 
??i1a12,  m&f,  &c.  will  therefore,  as  long  as  the  6's  are  small,  differ 
from  each  other  only  by  small  quantities  of  the  second  order. 

Field  of  external  force. 

101.     We  may  now  equate  to  zero  the  time  differential  of 
,  and  so  obtain,  applying  Boltzmann's  general  method, 


But 


dx  dy      0      dz 

8        - 


Again  let  the  forces  acting  on  the  system  have  a  potential 
Then 


da?         dy    0 
m2-^  =  -^,&c. 
dt         dxz 

And  the  equation  becomes,  arranging  the  terms, 


m  dx  daj  \     dy       m 


dz       m  dz  dy 

1      ,v    dDdb 
+  A  Z-TK  S2<a  -jr  TT^  +  &c- 


Let  us  first  consider  the  terms  which  involve  differentiation 
of  the  6's.  Since  every  6  is  a  function  of  the  corresponding 
r,  and  of  that  alone, 

db  _  d!6  dr  _  dft  x-x'  . 
dx     dr  dx     dr      r 

where  r2  =  (a?  -  #')2  +  (y  -  2/')2  +  («  -  ^)2- 
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dDdb  dDdbdr 


„„ 

Therefore  -35-  j~  ^r-  —  -7- 

db  ax  do  dr  dx 

dD  dblz  Xi—x<i         dD  dbls  x^  — 

=  Ofi  -77  --  1  -      --  1"  #1  ~ri  --  1  - 


~ri  --  1  - 

dblsdr13     rl3 
dD  dbl2  x1—x2     Q 


,  12  12  0 

because  -~a=—  -j-=  &c. 

dx1         cLx?, 

The  above  expression  can  by  taking  the  terms  in  pairs  be 
resolved  into  a  series  of  terms  of  the  form 


db  dr  r 

and  is  therefore  zero  on  average. 
Similarly 


can  be  resolved  into  a  series  of  terms  of  the  form 


and  is  therefore  also  zero  on  average. 
102.     We  have  then 


dx       m 

+  &c.  =  0, 
or  dividing  by  A, 


And  taking  mean  values  of  -~  from  Art.  97, 

da. 

fdlogA         dx     D    \          fdlogA      .  dy     D    \ 
ttl    ~^f"     ~  /l  T^  —  rr    +  Oa    "7-  -----  A-  -^  —  ^r-    +  &c.  =  0 

V    dxl  dxi  mj)^)         \    dxz  dx2  mj)*) 
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Andsinoe 


dx 


_  _ 


This  is  satisfied  by  assuming,  if  the  molecules  be  all  of 
the  same  kind, 

A  =€-*»*  ...........................  (37). 

We  have  then  for  the  law  of  distribution  of  coordinates  and 
velocities  in  a  field  of  external  force  whose  potential  is  % 


(-} 


-^  dxl  .  .  .  dznd^.^dyn  ......  (38). 

XT  — 

Now 


hD 

and  nT=2m 


the  system  consisting  of  n  molecules  of  the  same  kind. 
But  Q  =  a^2  +  612a^a2  +  &c. 

_  C^Al  +  ^18^12  +  ^2^2 

hD 

_3n 

Therefore  Q  =  pnT  ........................  (39). 


103.  With  regard  to  the  inclusion  of  the  intermolecular 
forces  under  the  potential  %  the  same  considerations  apply  as 
in  Chap.  ill.  I  shall  assume  that  we  may  use  the  form 


As  an  example  let  us  now  again  consider  the  system  of 
equal  elastic  spheres  in  vertical  column,  Arts.  29  and  45.  In 
that  case  ^  is  the  potential  of  the  constant  vertical  force.  Let 
us  suppose  K,  and  therefore  all  the  b  coefficients,  very  small.  If 
the  x  velocity  of  any  molecule,  as  ul}  be  given,  the  mean  value 


OF  THB 
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of  u2  is  —  •=—  ult  where  D  is  the  determinant 


D 


2          bl2        bls... 
b12        2          633... 


for  in  this  case  the  a  coefficients  may  each  be  replaced  by  unity. 

Now  Z>!2  contains  the  term  —  612  multiplied  by  axial  coefficients 
only,  and  other  terms  each  of  which  contains  the  product  of  at 
least  two  6's.  Therefore  if  each  b  be  small  enough,  the  mean 
value  of  u2  is  ultimately  —  ^612t^. 

Similarly  us  =  —  J613^!,  &c., 

and  therefore  -~-  —  ^m^u^  =  (2  —  J  26-)  ult 


or 

It  thus  appears  that  €~h™  differs  from  e~2/lx,  when  K  is  small, 
only  by  small  quantities  of  the  second  order.  Now  the  factor 
€-(K-KQ)f  or  €-*^  m  the  expression  for  the  density  in  Art.  29 
corresponds,  not  to  e"^*,  but  to  e"27^.  That  is,  if  -\/r  be  the 
potential  of  the  intermolecular  forces,  e~K  is  what  e~2/^  or  e"^^ 
becomes  in  the  limit  when  the  forces  become  infinite  for  the 
distance  r  =  c. 

And  K  being  small  the  solution  agrees  with  that  of  Art.  29. 
It  does  not  however  agree  with  that  of  Art.  29  when  K",  &c. 
cannot  be  neglected.  Which  result  then  are  we  to  accept 
in  this  case  ?  Evidently  the  result  e~^x,  because  Art.  29  is 
based  on  the  tacit  assumption  of  condition  A,  and  therefore 
holds  only  for  small  values  of  K.  To  sum  up  the  discussions  of 
this  problem. 

o 

If  K  be  neglected,  the  law  of  density  is  e~2hx,  with  T  =  -^  . 
If  K3,  &c.,  but  not  K,  be  neglected,  the  law  is  e~2ftxe~'c,  with 


If  /c2,  &c.,  are  not  to  be  neglected,  the  law  is  €-h^6~K,  with 

T—    ^ 

~ 
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Boltzmanris  Minimum  function. 

104.     The  following  theorem   is   not   to   be   compared   in 
egance  with  the  theorem  founded  on  condition  A,  Art.  38, 
but  has  the  advantage  that  condition  A  is  not  assumed. 

Let  f(ul  .  .  .  unxl  .  .  .  zn),  or  /,  be  any  positive  function  of  the 
coordinates  x^...zny  and  velocities  u^...un,  of  a  system  of 
molecules,  and  of  certain  parameters  a,  6,  &c.  Let  df  be  a 
small  variation  of/,  a,  b,  &c.,  being  constant.  Let  us  write 


Let  the  variation  df,  or  qft  be  made  subject  to  the  following 
conditions,  viz. 


II.. .f  1  +  qdxl . . .  dzndui . . .  dun 

=  11..  .fdxl . . .  dzndui . . .  dun, 
or  writing  dsdcr  for  the  continued  product  of  differentials, 

...fqdsd(7  =  0  (1). 

Also  ff. .  ./I  +  q  Cdsda-  =  ff . .  .fCdsdo; 

or  [j...fqCdsd<r  =  Q  (2), 

where  C  is  any  function  the  constancy  of  which  is  prescribed. 
Again,  let 
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We  have  then 

H  +  dH  =  I  j.  .  .  f  (log/-  1  )  dsdo- 

(/.../log  1  +  qdsda- 
jj...fqlogfdsda- 

1  1  •  •  •/<?  l°g  !  +  9  dsda- 
—  \\...fqdsd<T 


+  jj...fq\ogfdsd(T, 

because  (f  is  to  be  neglected. 

In  order  that  dH  may  be  necessarily  positive,  or  H  may  be 
least  possible,  subject  to  (1)  and  (2),  it  is  necessary,  and  it 

is  sufficient  that  1  1  fq  logfdsdv  =  0,  that  is  by  (2)  that  log  / 
shall  be  a  function  of  C,  that  is  that/=  Ae~hc,  where 

(L.fCdsdo- 


is  constant. 


We  may  write  f— Y  A*JD  for  A  and  Q  for  C.     The  actual 

\7T/ 


value  of  H  when  minimum  is 

n ,      h 
2l0^' 
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Physical  effect  of  the  b  coefficients. 

105.  From  the  law  of  distribution  at  which  we  have  now 
arrived  we  can  find  the  mean  value  of  any  of  the  products  of 
the  form  act',  &c.  For  instance,  a^  =  D^/hD,  D  being  the 
determinant  of  the  coefficients.  And  all  the  6's  being  negative, 
D12  is  necessarily  positive.  If  we  have  only  two  variables  x 
and  y,  and  if 

Q  =  aa?  +  bxy  +  cy*, 
b  being  negative,  we  find  D12  =  —  bt  which  is  positive. 

So  with  three  variables,  x,  y,  z,  and 

Q  =  arf?  +  a$f  +  a^z* 

+  \&y  +  bvapcz  +  bzgjz, 
we  find  -D12  =  -  2a3612  +  bj)^, 

which,  every  6  being  negative,  is  positive.    (See  Appendix  (6) 
as  to  the  sign  of  these  minors.) 

The  conclusion  to  be  drawn  from  this  is  that  if  two  mole- 
cules are  so  near  each  other  that  6  is  not  negligible,  they  are 
on  average  moving  in  the  same  direction.  There  is  at  every 
point  on  average  a  determinate  energy  of  stream  motion.  It 
follows  further  from  this  uniformity  of  direction  that  if  a#,  OLJ, 
/3y,  be  the  mean  values  of  those  functions  for  all  molecules 
within  any  small  finite  space,  (a/?)2,  &c.  have  greater  value 
than  they  would  have  if  the  molecules  were  material  points. 

Another  consequence  of  the  b  coefficients. 

106.  But  we  might  also,  treating  the  velocities  as  con- 
stants, proceed  as  follows. 

Let  for  a  moment  the  chance  that  two  molecules  distant  r 
from  each  other  shall  have  velocities  in  x,  u  ...n  +  du  and 
u'  .  .  .  u'  +  du  be  proportional  to 


cn 
in  which  b  may  be  q  —  ,  and  n  >  3. 
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It  is  convenient  here  to  treat  b  as  positive,  and  prefix  to  it 
the  minus  sign  in  the  index.  That  being  the  case,  if  b  be  not 
negligible,  that  is  if  r  be  small,  uu  is  more  likely  to  be  positive 
than  negative.  If  r  be  great,  and  therefore  b  negligible,  it  is 
indifferent.  It  follows  by  the  theory  of  inverse  probabilities 
that,  it  being  given  that  uu'  is  positive,  r  is  more  likely  to  be 
small  than  if  nothing  were  known  of  the  sign  of  uu'.  If  on  the 
other  hand  it  be  given  that  uu'  is  negative,  r  is  less  likely  to  be 
small. 

Suppose  a  molecule  at  P  with  velocity  u,  and  an  element  of 
volume  da)  at  P',  distant  r  from  P.  Let  the  chance  that  there 
shall  be  within  dw  a  molecule  moving  with  velocity  uf  be  ^dco, 
when  PP,  or  r,  is  great,  and  6  =  0.  Then  when  b  =(=  0,  the 
chance  becomes 


which  is  greater  or  less  than  -^dco,  according  as  uu'  is  positive 
or  negative. 

Further,  the  chance  that  there  shall  be  a  molecule  within 
da)  with  velocity  either  u  or  -  u  is,  since  a  priori  4-  u'  and  —  u 
are  equally  probable,  and  the  system  is  supposed  infinite, 
%^da>  (ehuu'b  +  e-huu'b), 

uu'  being  now  treated  as  essentially  positive.  That  expression 
is  equal  to  -v/rcta  if  62  and  higher  powers  of  b  are  neglected.  But 
if  62,  &c.  be  taken  into  account,  it  is  greater  than  ^dco. 

107.  It  follows  from  this  result  that  under  the  law  of 
distribution  of  coordinates  and  velocities  at  which  we  have 
arrived  in  this  chapter,  if  at  any  instant  we  calculated  for  each 
of  the  molecules  the  value  of  26,  or  the  sum  of  the  b  coefficients 
connecting  the  u  velocities  of  that  molecule  with  the  velocities 
of  each  of  the  others,  the  average  value  of  26  would  be 
greater  than  it  would  be  were  all  the  molecules  scattered 
through  space  at  random.  That  signifies  that  the  molecules  of 
the  entire  system  would  be  collected  into  denser  and  rarer 
masses  —  the  effect  of  which  evidently  is  to  increase  the  mean 
value  of  26  for  all  molecules.  Further  the  ratio  f  '2  :  £2  would 
have  a  greater  value  for  the  denser  masses  than  for  the  rarer 
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ones,  or  the  stream  motion  would  be  greater  for  the  denser 
masses  than  for  the  rarer  ones.  So  that  the  system  would  tend 
mOre  and  more,  with  increasing  number  of  molecules  in  a  given 
i  /space,  to  assume  the  form  of  a  number  of  denser  aggregates,  say 
clouds,  moving  through  a  comparatively  rare  medium. 

108.  It  is  not  essential  to  my  argument,  but  I  will  here 
give  what  appears  to  me  the  best  possible  definition  of  density 
at  a  given  point  at  a  given  instant,  for  molecules  of  finite 
dimensions. 

I  would  define  namely  as  follows.  The  density  at  any 
point  P  at  any  instant  is  —  ^b,  where  2$  is  the  sum  of  the 
b  coefficients,  all  being  negative,  of  all  the  molecules  of  the 
system  referred  to  the  point  P — that  is,  b  for  any  molecule  is 
a  function  of  the  distance  of  that  molecule  from  P. 

109.  It  is  assumed  throughout  this  investigation  that  the 
quadratic  function  Q  is  positive,  that  is,  that  the  coefficients 
a,  b,  are  such  that  Q  cannot  be  made  negative,  whatever  values 
we  assign  to  the  velocities.     It  is  assumed  that  the  distribution 
of  the  molecules  in  space,  and  therefore  the  coefficients  b,  are 
on  the  average  such  as  would  be  calculated  by  the  method  above 
indicated,  the  whole  number  of  molecules  of  the  system  in  a 
given  finite  space,  i.e.  p,  being  known.     Should  that  condition 
ever  fail,  the  motion  of  the  system  can  no  longer  be  represented 
by  the  law  considered   in  this   chapter.     It  is  reasonable  to 
expect   that   the   change   of  the   mathematical   formulae   will 
coincide  with  a  change  in  the  physical  state  of  the  system. 

110.  The  condition  that  Q  shall  be  necessarily  positive, 
whatever  values  be  assigned  to  the  velocities,  is  that  the  deter- 
minant D,  of  the  coefficients  in  Q  and  all  its  coaxial  minors  must 
be  positive.    See  Appendix  (p).    As  the  b  coefficients  increase  in 
absolute  magnitude,  D,  being  positive,  diminishes.     But  as  the 
density  p  increases  cceteris  paribus,  the  b  coefficients  do  generally 
increase  in  absolute  magnitude,  and  D,  with  all  its  coaxial 
minors,  generally  diminishes.     The  b  coefficients  are  generally 
also  functions   of  the  mean   kinetic  energy,  or   temperature, 
as  well  as  of  the   density,  and  if  the   temperature   be   high 
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enough,  probably  diminish  as  that  increases  (Art.  90).  As  the 
density  increases,  or  temperature  diminishes,  or  both,  the 
system  may  ultimately  arrive  at  a  state,  in  which  D  or  some  of 
its  coaxial  minors  becomes  zero  or  negative,  and  then  the  change 
of  physical  state  ensues,  whatever  its  true  description  be.  Now 
to  any  group  of  molecules  we  choose  to  select,  and  therefore  to 
all  the  molecules  forming  one  of  the  denser  masses  or  clouds 
above  mentioned,  belongs  a  particular  coaxial  minor,  that 
namely  obtained  by  striking  out  all  the  axial  constituents  of  D 
except  those  which  belong  to  molecules  of  the  group.  Gener- 
ally therefore  the  condition  may  be  expected  to  fail  for  some  of 
these  coaxial  minors  before  it  fails  for  the  complete  determinant 
D.  So  that  as  density  and  temperature  vary  continuously,  the 
change  of  state,  or  liquefaction  if  such  it  be,  will  take  place 
partially  and  gradually,  some  portions  of  the  system  having 
passed  into  the  changed  state,  while  other  portions  remain  in 
the  original  state.  Further,  the  higher  the  temperature,  at  all 
events  after  a  certain  point  is  reached,  the  greater  the  conden- 
sation necessary  to  make  D,  or  any  coaxial  minor  vanish,  and 
therefore  to  produce  the  physical  change  in  question. 

Such  is  the  process  which  our  analysis  leads  us  to  expect. 
Physicists  may  consider  how  far  it  corresponds  with  what  is 
known  to  take  place  in  gases  under  condensation,  or  on  what 
(if  any)  farther  hypothesis  it  may  be  made  to  correspond 
with  it. 


CHAPTER  IX. 


OF  MOLECULES  AS   CARRIERS. 

\/lll.  THE  molecules  being  elastic  spheres,  and  condition  A 
being  assumed,  to  find  the  mean  free  path  for  a  sphere  issuing 
from  collision  with  velocity  &>.  As  the  mean  free  path  has 
received  more  than  one  definition,  I  define  it  as  follows.  Sup- 
pose a  sphere  m  to  start  from  a  collision  with  velocity  G>,  and  to 
undergo  n  successive  collisions.  After  each  collision  let  each  of 
the  colliding  molecules  have  restored  to  it  the  velocity  which  it 
had  before  collision.  Then  our  sphere  will  describe  n  paths 
Z,  ,  £2  .  .  .  ln  between  its  collisions,  each  with  the  same  velocity  o>. 
The  mean  free  path  for  velocity  co  is  then 


when  n  becomes  indefinitely  great. 

If  i|r  be  the  velocity  before  collision  of  a  sphere  m'  which 
collides  with  m,  E  the  angle  between  the  directions  of  motion 


of  the  two  spheres  before  collision  as  in  the  figure,  the  relative 
velocity  of  approach  is 

R=  Vo>2  +  Tjr2  -  2onjr  COS  E. 

Let  c  be  the  sum  of  the  radii  of  m  and  m'.     Let/(i/r)  d^r  be 
the  number  of  m  spheres  in  unit  of  volume.    Then  the  number 
B.  8 
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of  collisions  which  m  undergoes  with  the  m'  spheres  per  unit  of 
time  is 


i  j    ( 


sin  EdE=>rrciNl  suppose, 

and  the  mean  free  path  for  an  m  moving  with  velocity  CD,  if  no 
collisions  occur  except  with  spheres  of  the  m'  class,  is  Xw  =  —  ^-  . 

If  there  be  two  classes  raj  and  ra.2,  and  m  can  collide  with  another 
sphere  of  either  class,  we  shall  have  a  corresponding  number  of 
collisions  per  unit  of  time  between  m  and  mz,  which  we  will 
denote  by  7rc2N2.  Then  tlje  mean  free  path  for  velocity  &>  is 

Xw  =  —  —  -^  -  vr-  and  so  on.    To  calculate  the  mean  free  path 

7TC2  (JV,  +  Nz) 

is  to  calculate  ^  -f  N2,  &c. 

112.     Another  method  of  obtaining  this  result  is  as  follows. 

Let  S  be  the  whole  space  in  which  our  spheres  are  moving. 
Consider  a  sphere  m  moving  with  velocity  o>,  and  another  sphere 
raj  whose  velocity  and  direction  of  motion  are  comprised  re- 
spectively within  the  limits  ty  ...  ty  +  dty  and  E  .  .  .  E  +  dE,  as 
in  the  figure. 

Now  ml  may  be  anywhere  within  8.  If  its  centre  be  at  this 
instant  within  the  cylinder  whose  base  is  a  circular  area  of  centre 
m  and  radius  c,  the  sum  of  the  radii  of  m  and  ml)  at  right 
angles  to  Rt  and  whose  height  is  Rdty  then,  but  not  otherwise, 
a  collision  will  occur  between  m  and  m^  within  the  time  dt  after 
this  instant. 

The  chance  that  such  collision  shall  not  occur  is  then 

7TC2Rdt 

S 

Now  the  number  per  unit  of  volume  of  spheres  ml  whose 
velocity  and  direction  of  motion  are  comprised  within  the  limits 
aforesaid  is  /(^)  d-fy  |  sin  EdE.  And  therefore  the  number 
within  S  is  Sf(ty)  dty  J  sin  EdE.  Let  us  now  assume  that  the 
aggregate  volume  of  all  the  ml  spheres  within  S  is  negligible 
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compared  with  S.     Then  the  chance  that  during  dt  m  shall  not 
collide  with  any  one  of  the  class  ml  is 

/-          7TC*Rdt\ 

(l     ~s~) 

raised  to  the  power  of 


that  is,  in  the  limit  when  dt  becomes  infinitely  small, 
1  -  7rc2Rdtf(^)  d^jr  %  sin  EdE, 

and  the  chance  that  m  shall  not  during  dt  collide  with  any  mly 
whatever  be  the  values  of  ty  or  E,  is  by  the  same  reasoning 

I    {'/ft)  d^^R  sin  EdE  =  1  -  ir(fdtN1  suppose. 


1  - 

Similarly  if  there  be  other  classes  of  spheres,  m?,  m^  &c., 
with  which  m  may  collide,  the  chance  that  it  shall  undergo  no 
collision  in  the  time  dt  is 

N2  4-  &c.)  =  1  - 


By  the  same  reasoning  the  chance  that  m  shall  undergo  no 
collision  in  n  successive  intervals  of  time  dt  is 

1  -  7TC*Nndt  =  1  -  w&Nt, 
if  t  =  ndt. 

Also  if  X  be  the  space  described  in  time  t 

\  =  a>t,   or   t  =  -  . 
o> 

Therefore  the  chance  that  m  shall  survive  for  a  distance  X 
without  undergoing  any  collision  is  1  —  7rczN  -  . 

Let  l- 

Then  <f>  is  the  proportion  of  the  whole  number  of  spheres  m 
which,  starting  with  velocity  o>,  traverse  the  space  X  without 
collision.  Also  <f>  (X  +  dX)  =  </>  (X)  <f>  (d\).  It  follows  that 


7    ,       -,    7 

=-&</>,  if  k= 

ay 

8—2 
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and  therefore  0  =  Ce~k\     Also  </>  =  1  if  X  =  0.     Therefore  0=1 
and  </>  =  e~k\     Whence  if  X  be  the  mean  value  of  X  given  o> 


-00 

=*/ 

Jo 


ft)  CO 

^N^ 


0 

which  is  the  same  result  as  we  obtained  above  (111). 
The  expression 

7TC2  (  d^lr  ("sin  EdEf(^r)  d^r  Vo>2  +  ^  -  2co^  cos  E 

Jo        Jo 

can  be  calculated  numerically  in  terms  of  ?rc2  and  h  if 


The  fundamental  equation. 

113.  Let   G  be  any  quality,  e.g.   mass,  colour,  charge  of 
electricity,  momentum  in  given  direction,  &c.  which  a  molecule 
may  possess  in  greater  or  less  degree,  and  which  it  can  carry 
with   it   unchanged   by   its   own   free   motion,   and  which,  as 
regards   the   aggregate   possessed  by  two   colliding   molecules 
is  not  changed  by  collision.     Suppose  a  quantity  of  gas  con- 
tained in  a  vertical  cylinder  of  which  the  height  is  very  small 
compared  with  the  diameter  of  the  base.     Every  molecule  that 
enters  the  cylinder  through  its  upper  flat  surface  shall  enter 
charged  with  the  quantity  GI  on  average  of  our  supposed  quality 
G.     Every  molecule  that  enters  the  cylinder  through  the  base 
shall  enter  charged  on  average  with  the  quantity  G0.    When 
the  motion  has  become  stationary,  every  molecule  within  the 
cylinder  at  a  height  z  above  the  base  will  be  on  average  charged 
with  a  quantity  of  the  thing  in  question  between  G0  and  GI. 
Let  it  be  called  G(z).     G(z)  in  stationary  motion   will  be  a 
function  of  z. 

114.  Now  consider  a  layer  Pl  of  molecules  in  the  cylinder 
between  the  planes  z1  and  zl  +  dz1.     If  the  molecules  within 
that  layer  were  enclosed  between  two  elastic  planes   z1   and 
zl  +  dzl  for  a  finite  time,  the  distribution  of  velocities  among 
the  molecules  within  it  would  become,  under  the  circumstances 
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supposed,  completely  independent  of  the  quantity  of  G  which 
they  happened  to  carry.  It  would  become  in  the  strictest  sense 
of  the  term  "  molecular  ungeordnet."  In  other  words  condition 
A  would  prevail. 

Collisions  would  occur  within  that  layer,  and  any  molecule 
emerging  from  such  collision  would  be  as  likely  to  be  moving  in 
any  one  direction  as  in  any  other.  The  chance  of  its  having 
any  particular  direction  would  be  independent  of  its  charge  G. 

Consequently  if  any  such  molecule,  having  undergone  col- 
lision within  the  layer  PJ  arrives  without  further  collision  at 
the  layer  P,  the  mean  quantity  of  G  which  it  brings  with  it  is 
the  mean  quantity  of  G  per  molecule  for  all  the  molecules 
within  the  layer  at  Px.  This  is  the  fundamental  assumption. 

115.  Now  let  Gz  be  the  mean  value  of  that  quantity  for 
all  the  molecules  within  the  layer  z  .  .  .  z  +  dz  at  P.  Suppose  a 
molecule  to  arrive  at  that  layer  from  above  with  velocity 
between  v  and  v  +  dv,  and  in  direction  making  angle  between 
0  and  6  +  d6  with  the  vertical.  Let  X  be  the  distance  which 
on  its  arrival  at  P  it  has  traversed  since  its  last  collision.  Then 
such  last  collision  must  have  occurred  at  a  height  X  cos  6  above 
P.  Therefore  the  quantity  of  G  with  which  the  molecule  is 
charged  is  the  mean  quantity  of  G  for  the  height  z  +  X  cos  6, 
that  is 


because  X  cos  0  is  small. 

The  number  of  molecules  which,  having  velocity  and  direction 
between  those  limits,  pass  through  the  horizontal  plane  P  per 
unit  of  area  and  time  is  vf(v)  £  sin  9  cos  Odddv. 

Hence  the  quantity  of  G  carried  through  the  layer  P  per 
unit  of  area  and  time  by  molecules  coming  from  above  is 

7T  ff 

f2  rz  fin. 

vf(v)  dv%     sin  6  cos  0Q  (z)  dO  +  ±1   X  cos2  6  sin  OdO  ^  f(v)  vdv 
J  o  Jo  &2 
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and  taking  the  mean  of  all  values  of  X  and  v  it  is 


Here 


,oo 

=  I   f(v)v\dv. 
Jo 


By  the  same  reasoning  the  quantity  of  G  carried  through 
the  layer  P  per  unit  of  area  and  time  by  molecules  coming 

from  below  is  \G  (z)  —  J  X  -7-  . 

And  the  quantity  passing  through  P  downwards  exceeds 
the  quantity  passing  through  P  upwards  per  unit  of  area  and 

time  by  JX. 


116.  This  is  the  fundamental  equation  from  which  Boltz- 
mann  calculates  the  rate  of  diffusion,  viscosity  &c.  in  gases  as- 
suming condition  A  to  prevail.   It  would  represent  with  complete 
accuracy  the  initial  motion  of  the  system,  if  the  quality  G  were 
suddenly  attributed  to  the  molecules  in  the  degree  proper  to  z, 
but  without  regard  to  their  velocities. 

It  takes  no  account,  as  he  points  out,  of  any  variation  in  the 
distribution  of  velocities  which  may  be  caused  by  the  quality 
G  itself,  for  if  for  instance  G  be  momentum  in  given  direction 
at  right  angles  to  the  axis  of  our  cylinder,  such  a  change  in 
distribution  must  exist.  It  is  assumed  that  the  effect  of  this 
change  in  distribution  on  the  mean  free  path  X  is  so  small  as  to 
be  negligible,  a  condition  which  must  be  secured  if  the  given 
terminal  conditions,  e.g.  GQ  and  6rl5  do  not  differ  very  widely 
per  unit  of  distance. 

117.  The  method  also  takes  no  account  of  the  stream.   The 
effect  is  to  set  up  a  stream,  so  to  speak,  of  G  from  the  upper  to 
the  lower  parts  of  the  cylinder.     And  that,  as  will  be  found, 
causes  a  small  error  in  the  result.     The  error  may  indeed  be 
negligible  so  long  as  the  stream  velocity  is  very  small  compared 
with  molecular  velocities,  which  will  be  the  case  if  the  terminal 
conditions  do  not  differ  widely  per  unit  of  distance. 
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It  therefore  does  not  sensibly  affect  the  results  calculated 
from  Boltzmann's  formula  in  any  case  to  which  as  I  understand 
he  proposes  to  apply  that  formula.  It  is  worth  while  however 
to  consider  the  sign  and  order  of  magnitude  of  this  deviation 
from  complete  accuracy. 

We  are  asked  to  suppose  that  in  stationary  motion  molecules 
carrying  positive  G,  have  a  mean  velocity  u  downwards ;  mole- 
cules carrying  negative  G  have  a  mean  velocity  u',  u  and  u 
being  functions  of  z.  Hydrodynamical  considerations  make  Gu 
and  Gu'  sensibly  constant  throughout  the  cylinder  in  stationary 
motion. 

But  when  collisions  take  place  in  the  layer  P1}  condition  A 
is  supposed  to  prevail  so  that  the  molecules  issuing  from 
collision  are  as  likely  to  be  moving  in  any  one  direction  as  in 
any  other,  and  to  carry  with  them,  in  whatever  direction  they 
do  move,  the  quantity  of  G  due  to  the  height  z. 

In  fact  molecules  with  positive  G,  having  before  collision 
the  mean  velocity  u  downwards,  will  issue  from  collision  with 
some  mean  velocity  downwards,  less  indeed  than  u  but  com- 
parable with  it,  and  therefore  will  not  be  moving  in  all  directions 
indifferently,  but  on  the  average  downwards.  And  molecules 
moving  downwards  will  possess  on  average  a  greater  quantity 
of  G  than  that  due  to  the  layer  in  which  their  last  collision 
took  place.  For  let  co  be  the  absolute  velocity  of  a  molecule 
ra  before  collision.  Let  ty  be  the  velocity  before  collision 
of  the  other  molecule  m  concerned,  E  the  angle  between 
their  directions  as  in  the  figure  (Art.  111).  Then  the  velocity  of 
m  after  collision  is  the  resultant  of  (1)  the  common  velocity  of 
m  and  m,  (2)  their  relative  velocity, in  whatever  direction  it  be, 
after  collision.  But  for  the  relative  velocity  after  collision  all 
directions  are  equally  probable.  The  relative  velocity  therefore 
contributes  on  average  nothing  to  the  resultant. 

The  resultant  velocity  after  collision  is  then  on  average 
the  common  velocity.  And  the  velocity  after  collision  re- 
solved in  its  direction  before  collision  is  the  common  velo- 
city resolved  in  direction  co,  that  is  if  m  =  m'  co  +  i/r  cos  E. 
We  have  then  for  the  mean  value  of  the  velocity  after 
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collision   resolved    in   the   direction   before   collision   the   ex- 
pression 


I    d^e-h**  W  I  \/a)z  +  -»/r2  -  2ft>^  cos  E  (&>  4-  ^  cos  E)  sin  EdE. 

Jo  Jo 


A 


That  can  be  calculated  as  a  function  of  &>.  As  however  we 
intend  to  treat  it  as  negligible,  it  is  not  worth  while  to  calculate 
it  here.  But  we  can  easily  prove  that  it  is  necessarily  positive. 
We  may  therefore  conclude  that  the  mean  velocity  after  col- 
lision of  a  molecule  carrying  positive  G  resolved  in  its  direction 
before  collision  is/to  on  average,  where/"  is  some  positive  quantity 
less  than  unity. 

It  follows  that  the  mean  velocity  downwards  of  the  mole- 
cules carrying  positive  G,  which  before  collision  was  u,  is  after 
collision  fu. 

We  have  thus  proved  that  the  molecules  carrying  positive 
G  issue  from  collision  in  the  layer  Pj  with  mean  velocity  fu 
down  the  cylinder.  For  the  same  reason  the  molecules  carrying 
negative  G  issue  from  collision  with  mean  velocity  fu0  upwards. 
Therefore  of  the  molecules  which  reach  the  layer  P,  having  had 
their  last  collision  in  the  layer  Plf  the  proportion  which  carry 
positive  G  is  not  the  same  as  the  proportion  of  all  the  molecules 
within  the  layer  Pl  which  carry  positive  G,  as  it  was  assumed 
to  be  in  forming  the  fundamental  equation,  but  is  greater. 

118.  The  effect  of  the  introduction  of  the  b  coefficients  in 
this  case. 

It  has  been  shown  in  Chapters  IV.  —  VI.  that  molecules  near  to 
one  another  are  on  average  moving  in  the  same  direction.  The 
energy  of  their  relative  motion  is  therefore  less  than  it  would 
be  were  condition  A  assumed,  and  therefore  their  encounters 
pro  tanto  less  frequent.  That  tends  to  increase  the  mean  free 
path. 

On  the  other  hand  under  assumption  B  the  molecules  tend 
to  become  aggregated  into  denser  and  rarer  masses.  That 
tends  to  increase  the  frequency  of  their  encounters  and  there- 
fore to  diminish  the  mean  free  path. 
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It  is  probable  that  for  ordinary  gases  under  ordinary  con- 
ditions the  b  coefficients  will  be  very  small.  And  if  so  the 
effect  of  the  b  coefficients  on  the  phenomena  of  diffusion  &c. 
will  be  of  the  same  order  of  magnitude  as  those  sources  of  error 
mentioned  in  Arts.  116,  117,  which  as  we  have  said  Boltzmann 
and  others  have  agreed  to  neglect  in  calculations  concerning 
these  phenomena.  I  shall  not  therefore  here  further  consider 
the  effect  of  the  use  of  assumption  B  instead  of  A  on  such 
calculations. 


CHAPTER  X. 

ON   THEJYPOTHESIS  THAT   HEAT   OR   TEMPERATURE   IS   REPRE- 
SENTED  BY  THE   KINETIC   ENERGY   OF   MOLECULAR   MOTION. 

Natanson's  Theorem. 

119.  M.  Ladislas  Natanson  (Interpretation  cinetique  de  la 
fonction  de  dissipation — Bulletin  de  I'Academie  des  Sciences  de 
Cracovie,   Decembre   1893)  considers  a  medium  composed    of 
molecules  in  motion.    The  components  of  the  molecular  velocity 
of  a  molecule  are  u,  v,  w,  those  of  the  apparent  or  stream 
velocity  of  an  element  of  volume  are  f ,  77,  f.     I  have  here  inter- 
changed Natanson's  symbols  to  make   them  agree  as   far   as 
possible  with  my  own,  and  so  avoid  confusion.  Evidently,  taking 
mean  values, 

u  =  v  =  w  =  0. 

Let  K  =  i  [([ p  (f 2  +  772  +  ? 2)  dxdydz, 

p  denoting  density,  and 

E  =  J  1 1 1  p  (u2  +  ?j2  +  w2)  dxdydz. 

The  integrations  are  throughout  all  space  occupied  by  the 
system. 

120.  M.  Natanson  gives  no  definition  of  f,  77,  f,  but  is,  as  I 
understand,  willing  to  accept  Boltzmann's  definition — that  f  at 
P  is  the  momentum  in  x  of  the  centre  of  inertia  of  all  the 
molecules  contained  in  a  sphere  of  radius  r  described  about  P 
as  centre  divided  by  the  volume  of  that  sphere,  in  the  limit 
when  r  becomes  infinitely  small. 

I  should  propose  myself  to  define  it  as  in  Art.  61  by  the 
function  /,  so  that  f ,  77,   f,  Natanson's  components  of  stream 
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velocity,  agree  nearly  with  the  f,  77,  f,  of  Chapter  v.  But 
M.  Natanson  does  not  on  his  own  assumption  require  any 
further  definition. 

The  normal  pressures  per  unit  of  surface  are 

pxx  =  pu*,    pyy  =  P»*>    Pzz 
The  tangential  pressures  are 


121.  He  then  gives  the  following  fundamental  equation, 
viz.  q  being  any  function  of  u  4-  f  ,  v  4-  rj,  w  4-  f,  the  time  varia- 
tion of  q  is  found  from 

da      d  ,   —  x      d  ,  —  v      d 
- 


r  dt 

Here  X,  Y,  Z  are  the  components  of  external  force,  and  -^ 

da  dt 

the  change  of  q  due  to  encounters,  -^  the  total  change  of  q 

with  the  time. 

Writing  u  +  f  for  qy  he  finds 


And  again,  putting  q  =  (u  +  f  )2  -4-  (v  +  rjf  +  (w  +  f  )2,  he  obtains, 
neglecting  terms  of  the  third  order, 


T-  (2f/3^  +  Zrjpwv  4-  2fp  ww 


-j- 
ay 

4-  T-  (Zgpuw  +  fypvw  4-  2£pw-) 


(5). 
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Combining  (4)  and  (5)  we  get 


2  +     +     )    ........................  (6). 

Let  f  =  a,       f  =  b,       f=c, 

dx  dy  dz 

dv  +  ^-A     *  +  *-/?      d*  +  dr>    r 

dz^  dy~  dx+~dz-">       dy  +  dx~ 

Add  to  the  first  member  of  (6) 


which  is  zero.     Also  write 


And  then  in  the  equation  so  obtained  integrate  for  vc,  y,  z, 
and  make  1%  +  mrj  +  ?if=0  at  the  bounding  surface,  Z,  m,  n 
being  direction  cosines  of  the  normal  to  that  surface.  Assume 
further  that  X  =  Y  —  Z  =  0,  or  no  external  forces  act. 

The  result  is 
TT  +  1  1  \(pu?a  +  pv*b  +  pw*c  +  pvwA  4-  puwB  +  puvC)  dxdydz  =  0. 

Treating  equation  (5)  in  the  same  way,  we  obtain 
-j-  —  III  (pu2a  +  pv2b  +  pw*c  +  pvwA  +  puwB  +  puvC)  dxdydz  =  0. 

And  therefore,  there   being   no    external   forces,  —  r-  =  —  -r-  , 

or  the  kinetic  energy  of  stream  motion  can  increase  or  diminish 
only  at  the  expense  of  molecular  energy. 
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122.     Let  us  now  write 


=  pu* 


F  =  (p-  pu?)  a  +  (p  -  pv*)  b  +  (p  - 

—  /jvw^L  —  puwB  —  pttvC. 
By  a    further  use  of  the   fundamental   equation,  putting 


q  =  (u  +  f)2,  and  then  q  =  v  +  ijx  w  +  £,  we  find 

JO  

4p3£ 

+  2  (/ow;)2  + 
And  finally 


where  ,_a  +  6  +  fl_       +       +      . 

123.  It  appears  then  that  ^.  the  molecular  or  heat  energy, 
gains  or  loses  at  the  expense  of  the  energy  of  stream  motion, 
according  to  the  law  of  encounters  "  loi  des  chocs  moleculaires." 
If,  that  is,  as  the  result  of  encounters  the  tangential  pressures 
are  increasing,  E  diminishes,  and  K  increases.  Now  this  is 
precisely  what  happens  according  to  Chap.  v.  ante,  when, 
for  instance,  instead  of  being  material  points,  the  molecules 
become  elastic  spheres  of  finite  diameter.  As  the  result  of 
collisions  in  that  case  f2  in  my  notation  becomes  f2  +  f  /2,  &c. 
And  as  we  have  seen,  Art.  105,  («7)2,(a/3)2,  and  (@y)2  increase, 
but  only  up  to  a  certain  point,  at  which  namely  the  ratio  f  /2/f* 
has  the  proper  value  for  stationary  motion.  The  relation 
between  the  energy  of  visible  motion  K  and  that  of  molecular 
motion  E  is  in  Natanson's  theorem  the  same  as  that  between 
K  and  heat  in  Thermodynamics. 

i/ 
Messrs  Bryan  and  Boltzmanris  method. 

124.    In  order  to  confirm  the  hypothesis  that  the  temperature 
of  a  gas  or  other  substance  is  proportional  to  the  mean  kinetic 
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energy  of  its  molecules,  Messrs  Bryan  and  Boltzmann  (Vienna 
$itzungsberichte  math,  natumu.  Classe,  Band  III.,  Abtheil  II. 
/Dec.  1894,  also  referred  to  by  Boltzmann,  Gastheorie,  p.  136,  et 
seq.)  treat  the  following  case.  X  and  Y  are  two  infinite  parallel 
planes.  In  the  space  Si  to  the  left  of  X  is  a  gas  A.  In  the 
space  $2  to  the  right  of  Y  another  gas  B.  In  the  space  S 
between  X  and  Y  there  acts  on  molecules  of  gas  A  a  force 
directed  from  Y  whose  potential  becomes  infinite  on  Y,  and 
there  acts  on  molecules  of  gas  B  a  force  directed  from  X  whose 
potential  becomes  infinite  on  X.  These  conditions  insure  that 
no  molecule  of  gas  A  will  be  found  in  $2,  and  no  molecule  of  gas 
B  will  be  found  in  $lt  But  in  the  space  S  there  will  be  a 
mixture  of  A  and  B,  and  encounters  will  take  place  between 
their  respective  molecules. 

Messrs  Bryan  and  Boltzmann  use  generalised  coordinates, 
but  it  will  be  sufficient  in  this  notice  to  deal  only  with  velocities 
of  translation,  as  though  both  A  and  B  molecules  were  elastic 
spheres. 

It  is  assumed  that  the  mean  kinetic  energy  of  A  molecules 
at  any  point  within  S  is  the  same  as  their  mean  kinetic  energy 
within  8it  and  that  the  mean  kinetic  energy  of  B  molecules 
within  8  is  the  same  as  their  mean  kinetic  energy  within  $2. 

It  is  assumed  further  that  as  the  result  of  encounters  taking 
place  within  S,  between  A  and  B,  the  molecules  of  A  acquire,  if 
they  have  it  not  already,  the  same  mean  kinetic  energy  as  those 
of  B.  That  result  seems  to  me  to  be  questionable  or  to  be  only 

VV,  . 

'Vapproximately  true,  if  the  molecules  have  finite  diameters  or 
spheres  of  action.  But  I  state  this  with  diffidence  when  such 
authorities  are  against  me. 

Making  however  these  two  assumptions,  the  conclusion 
follows  that  the  mean  kinetic  energy  for  the  A  molecules  in  S1 
must  become  equal  to  that  of  the  B  molecules  in  $2.  And 
therefore  that  if  the  mean  kinetic  energy  be  originally  different 
for  A  and  for  B,  a  finite  quantity  of  energy  will  pass  from  Si  to 
$2 ,  or  vice  versa. 

The  condition  that  no  energy  shall  pass  is  that  the  mean 
kinetic  energy  shall  be  the  same  for  A  as  for  B.  But  we  may 
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regard  the  two  gases  as  two  bodies  in  contact,  which  condition 
is  more  nearly  realised  as  S  is  diminished  indefinitely.  The 
known  experimental  condition  that  no  energy  shall  pass  between 
fc/o  bodies  in  contact  is  that  both  shall  be  at  the  same  tempera- 

i  Xure.     Hence  in  the  case  of  the  two  gases  at  all  events  kinetic 

"  energy  plays  the  same  part  as  temperature. 

Professor  J.  J.  Thomsons  method. 

125.  The  same  hypothesis  is  tested  in  a  different  way  by 
Professor  J.  J.  Thomson  (Application  of  Dynamics  to  Physics  and  Sf 
Ghemistry,  p.  91).  He  puts  the  following  case.  Suppose  a  tube  in 
which  is  moveable  a  piston  of  very  great  mass  M.  On  one  side 
of  it  is  a  set  A  of  material  particles  each  of  mass  m^  in  motion. 
On  the  other  side  a  set  B  of  material  particles  each  of  mass  m2  in 
motion.  Considering  the  two  systems  of  particles  as  two  sub- 
stances, and  the  piston  as  a  conductor  of  heat,  the  condition 
that  no  heat  shall  pass  through  the  piston  from  one  substance 
to  the  other  is  that  the  two  substances  shall  be  at  the  same 
temperature.  If  temperature  be  represented  by  kinetic  energy 
of  translation,  then  the  condition  that  neither  set  of  particles 
shall  on  average  either  lose  or  gain  energy  by  collisions  with 
the  piston  ought  to  be  that  the  mean  kinetic  energy  of  transla- 
tion is  the  same  for  the  A  as  for  the  B  set. 

Let  U  be  the  velocity  of  the  piston,  v^  the  velocity 
normal  to  the  surface  of  the  piston  of  a  particle  m1  before 
impact.  It  is  assumed  that  particles  striking  the  piston  re- 
bound from  it  as  perfectly  elastic  bodies.  Let  U',  u^  denote 
the  values  which  U,  u^  assume  after  impact.  Then  we  have 

+  M  U  =  ra^'  +  M  U', 
+  MU*  =  m&i*  +  M  U'\ 
From  which  deduce 

-  4,Mmi  (M  -  mj  Uu^  -  IMmfuf  ,  . 


Now  add  together  the  equations  of  this  form  for  all  collisions 
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which  take  place  in  unit  of  time.  Let  Nj_  be  the  number  of 
such  collisions.  Also  let  0  =  ^MUZ  be  the  kinetic  energy  of  the 
piston,  and  0l  the  mean  kinetic  energy  of  the  particles  of  the 
set  A  which  strike  the  piston.  As  I  understand  Professor 
Thomson,  the  mean  is  taken  for  all  collisions,  not  for  all 
particles. 

It  is  assumed  that  the  velocity  U  of  the  piston  is  zero 
on  average  of  time,  because  the  piston  is  being  struck  on  the 
other  side  by  particles  of  the  B  set  and  the  pressure  is  supposed 
to  be  the  same  on  either  side.  Hence  we  may  neglect  in  the 
expression  (a)  the  term  which  contains  U  in  the  first  degree. 

The  expression  then  takes  the  form 

...............  (i). 


NOTE.  This  should,  I  think,  be  §^  instead  of  J^  because  as  MU2=  2<9t 
so  m>iUlz  =  %61,  &i  denoting  the  whole  mean  kinetic  energy  of  all  molecules 
which  collide  with  the  piston  or  with  other  portions  of  the  elastic  surface 
bounding  the  A  molecules,  i.e. 


But  this  does  not  affect  Professor  Thomson's  argument. 
Similarly  for  the  B  set  we  have 


and  since  90  =  —  N1d6l  —  Nf>6z  =  0  on  average 

t  2Mm2 

V    -W  +  ^          mfW-W* 

N1ml 


This  gives          20  -  Jft  =  — -r  J  (6,  -  ft), 
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Substituting  in  (1)  and  (2) 

,         2Mab 


2Mab 

~(*~  ll 


and  therefore  3^  =  302  =  0,  if  ft  =  02. 

It  follows  that  if  02  >  0t  the  effect  of  all  the  collisions  is 
to  increase  6l  at  the  expense  of  02,  and  vice  versa.  And  the 
condition  that  no  energy  shall  pass  from  A  to  B  or  from  B  to  A 
is  that  #2  =  #1,  or  the  mean  kinetic  energy  of  translation  for  all 
particles  colliding  with  the  piston  is  ffifi  aa-me  for  ,4.  as  for  B. 


126.  Now  let/(V)cfa£  be  the  number  per  unit  of  volume  of 
particles  of  the  A  set  whose  velocities  normal  to  the  piston  lie 
between  u  and  u  +  du.  Then  the  number  of  such  particles 

/•» 
which  strike  the  piston  per  unit  of  area  and  time  is  /    f(u)  udu. 

Jo 

The  sum  of  the  kinetic  energies  of  their  motion  normal  to  the 
,*> 

piston  is  Jraj  I    f(u)u3du.     The  mean  kinetic  energy  of  motion 
Jo 

normal  to  the  piston  for  all  A  particles  irrespective  of  their 

-00 

striking  the  piston  is  ^m^  I   f(u)  u*du.     Let  f(u)  be  so  chosen 

Jo 
that  the  ratio 


^•00  /*°^ 

I  f(u)u3du     I  f(u)u2du 
^o .  Jo 

/•OO  •  (-00  > 

I  f(u)udu        I  f(u)du 
Jo  Jo 


is  an  absolute  constant,  and  therefore  the  same  for  A  as  for  B. 
That  being  the  case,  when  6l  =  6z,  or  the  mean  kinetic  energy 
per  collision  is  the  same  for  A  as  for  B,  the  mean  kinetic  energy 
per  particle  is  also  the  same  for  A  as  for  B.  And  therefore  the 
condition  that  on  the  whole  there  shall  be  no  transfer  of  energy 
through  the  piston  from  A  to  B  or  vice  versa  is  that  the  mean 
kinetic  energy  per  particle  is  the  same  for  A  as  for  B. 

B.  9 
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.     If  we  assume  for  A  particles 


and  for  E  particles        /(tia)  =  C 
we  satisfy  the  condition. 

Possibly  there  may  be  other  forms  of  f(u)  which  satisfy  it. 
But  however  this  may  be,  some  assumption  with  regard  to  the 
relative  frequency  of  different  values  of  u  appears  to  be  necessary. 

Further,  f(u)  being  properly  chosen,  we  satisfy  at  the  same 

time  the  two  conditions  (1)  that  the  pressure  shall  be  the  same 

O,      on  either  side  of  the  piston,  the  number  of  particles  per  unit  of 

volume  being  the  same  for  A  as  for  B,  and  (2)  that  there  shall 

be  no  transfer  of  energy. 

The  law  of  distribution  of  velocities  thus  appears  to  be 
-i  of  the  essence  of  the  matter,  if  the  .molecular  motion  is  to 
/  represent  heat. 

The  Second  Law  of  Thermodynamics. 

128.  Our  gas  system  being  in  stationary  motion  according 
to  the  laws  above  investigated,  we  might  impart  to  it  a  certain 
quantity  of  heat  BE.     That  being  done,  it  would  assume  a  new 
form  of  stationary  motion,  in  which  the  mean  kinetic  energy  T, 
the  volume  co,  and  the  mean  potential  energy  %,  will  be  gener- 
ally different  from  what  they  were  in  the  original  motion.     If 
P  be  the  external  pressure  which  in  the  original  motion  was 
just  sufficient  to  prevent  expansion,  the  change  of  volume  dco 
can  be  effected  only  by  doing  an  amount  of  work  P9o>  against 
the  external  force  P.     A  part  of  the  energy  dE  supplied  to  the 
system  is  spent  in  doing  this  work.     The  remainder  is  spent  in 
increasing  either  (1)  the  mean  kinetic  energy  T,  or  (2)  the 
mean  potential  energy,  ^,  of  the  intermolecular  forces,  so  that 
dE  —  A9T  +  9%  +  P9&>,  X  being  the  number  of  molecules. 

129.  If  the  external  forces  P  also  have  a  potential,  there  is 
in  general  a  determinate  relation  between  dco  and  dT,  and  a 
determinate  relation  between  day  and  9^.     So  that  in  this  case 


X.]  THERMO  DYNAMICAL   RELATIONS.  LSI 

there  is  only  one  independent  variable,  and  no  question  can 

7)  J? 
arise  whether  or  not  -^   is  a  complete  differential.    Any  theorem 

proved  on  the  hypothesis  of  P  having  a  potential  is  not  the 
Second  law  of  Thermodynamics  as  usually  understood. 

But  the  external  force  P  is  in  Professor  Thomson's  language 
a  "controllable"  force.  It  is  possible  therefore  by  suitably 
varying  P  to  maintain  the  system  in  stationary  motion  with 
any  arbitrarily  assigned  values  of  T  and  o>.  But  dT  and 
8ft)  being  given,  3%  is  determinate.  There  are  then,  P  being 
controllable,  two  independent  variables.  It  is  convenient  to 
use  P  and  co  for  independent  variables. 

130.  If  in  rectangular  coordinates  the  abscissa  represents 
a),  and  the  ordinate  represents  P,  any  point  in  the  diagram 
represents  a  determinate  state  of  the  system,  in  which  it  is  in 
stationary  motion  with  the  values  of  P  and  ft)  corresponding  to 
the  point,  it  being  understood  that  the  controllable  force  P  is  to 
have  such  value  given  to  it  as  will  be  necessary  to  maintain 
that  motion  stationary.  So  also  any  curve  drawn  on  the  plane 
of  the  diagram  represents  a  series  of  states  through  which  the 
system  might  be  made  to  pass  successively,  P  being  made  to 
vary  as  required.  But  it  is  understood  that  the  change  of 
state,  that  is  the  passage  from  one  point  in  the  curve  to 
another,  is  effected  so  slowly  that  stationary  motion  is  always 
attained,  and  therefore  the  kinetic  energy  corresponding  to  the 
controllable  coordinate  co  is  for  our  purpose  to  be  taken  as  zero. 

If  the  curve  be  a  closed  curve,  the  system  having  passed 
through  all  the  states  denoted  by  the  curve,  has  at  the  end 
of  the  process  the  same  values  of  P  arid  ft),  and  is  therefore 
for  all  experimental  purposes  in  the  same  condition  as  it  was 
at  the  beginning.  Contemplating  the  system  itself  alone, 
nothing  that  we  can  observe  by  experiment  has  happened  to  it. 
Now  if  6  denote  the  absolute  temperature,  the  second  law  of 

/•  £\  TT» 

Thermodynamics  asserts  that  in  this  complete  cycle  I  —  =  0, 

r\T? 

or      -  is  a  complete  differential  of  some  function  of  0  and  co. 
u 

9—2 


132  KINETIC   THEORY    OF   GASES.  [CHAP. 

That  function  was  called  by  some  writers  the  Thermodynamic 
function,  by  others  the  Entropy  of  the  system. 

r)7? 

131.  As  I  have  said,  -^-  being  a  complete  differential,  no 
change  takes  place  in  the  system  itself  as  the  result  of  any 
complete  cycle  of  operations  of  the  kind  described.  But  wE^O, 

and  jPdft)  =j=  0-  Generally  a  quantity  of  heat  denoted  by  IdE, 
drawn  from  some  external  source,  has  been  converted  into 
mechanical  work  iPda)  outside  of  the  system,  and  a  further 

quantity  has  been  transferred  from  a  hotter  to  a  colder  external 
body,  as  shown  in  treatises  on  Thermodynamics.  But  if  our 
system  be  enclosed  in  a  non-conducting  envelope,  so  that  no 
heat  or  other  form  of  energy  is  allowed  to  pass  into  it  from 
without,  or  out  of  it  into  external  space,  it  may  by  suitably 
varying  the  controllable  P  be  made  to  expand  or  be  compressed, 
and  in  so  doing  it  does  external  work,  positive  in  case  of 
expansion,  negative  in  case  of  compression.  This  work  is  done 
at  the  expense  either  of  the  absolute  temperature  6  or  of  the 
potential  ^  of  the  system  itself.  The  line  traced  on  the 
diagram  in  such  a  process,  a  line  that  is  for  which  dE  =  0,  is 
called  an  adiabatic  line.  If  the  cycle  be  A  BCD  on  the  diagram, 
AB  denoting  expansion  at  constant  absolute  temperature  0,  BG 
adiabatic,  CD  compression  at  constant  absolute  temperature  & ', 
and  DA  adiabatic,  it  is  found  that  the  heat  which  must  be 
supplied  from  without  to  effect  the  expansion  AB  is  to  the 
heat  withdrawn  during  the  compression  CD,  both  measured  in 
mechanical  units,  as  6  :  6'.  If  temperature  is  represented  by 

O   ITT 

the  kinetic  energy,  we  ought  to  find  T  :  T'  : :  0  :  0',  or  ~  is 
a/complete  differential. 

/       132.     In    the   Philosophical   Magazine,   January,   1876,   a 

\/       Pro°f  of  this   proposition  was   given   by  the   present   writer. 

This    was    employed    afterwards    in    an    improved    form    by 

Dr   Watson  in  his   Kinetic   Theory   of  Gases,   and    has  been 

accepted  by  Professor  Bryan.     It  was  founded,  as  will  be  seen, 
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on   Boltzmann's   theorem,   that  the  chance   of  a   group  of  X 
molecules  having  their  coordinates  between  the  limits 


(in  which  position  their  potential  is  ^)  is  proportional  to 


This  is  established  only  on  the  assumption  of  condition  A. 
I  will  here  give  that  proof  of  the  second  law  in  order  that  we 
may  subsequently  see  what  modifications  are  required  in  it 
when  we  abandon  condition  A,  and  employ  the  more  general 
method  of  this  work.  I  here  use  h  instead  of  2h  of  Art.  44. 

We  have  dE  =  KdT  +  3 %  +  Pdco. 

Let  T=yj-=  the  mean  kinetic  energy  of  a  molecule. 

The;  '-» 


Aa/Fi        . 
and  T         3~  dT  +  ^dx  +  ^Pdv  ...............  (I). 

Now  substitute  for  Pd(o  from  Clausius'  equation,  Art.  24, 
f  Pa  =  24ra  (w2  +  v2  +  O  +  422  Er  , 

or  Pdco  =  &m  (tf  +  v2  +  O  —  +  J22lfr  —  . 

ft)  Ct) 

Now  if  there  be  a  general  cubical  expansion,  o>  becoming 
9o>,  every  line  r  in  the  system  becomes 


, 
or  dr  =  4  —  r. 

3    Q) 

Therefore  the  last  equation  becomes 

Pdco  =  £2ra  (u2  +  v2  +  w*)  —  + 
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Again,  for  each  configuration 

22Rdr=-dx  =  -CT^dco. 

dco 
Hence  taking  mean  values 


Here  -^  is  the  mean  for  all  configurations  of  ^—  ,  that  is 
the  mean  of  a  variation,  while  9%  is  the  variation  of  a  mean. 

The  distinction  between  9y  and  ~  dco  is  of  the  essence  of  the 

dco 

matter. 
Again, 

iSmO"^*1  +  «f)-  =  T-. 
7  co       h  co 

if  X  be  the  number  of  molecules  of  the  system. 

We  have  then 

D^        X  dco      dy^ 

PdcO  =  y-  ---  -£  d(0, 

fi  co       dco 

'2fi  dco      2h(fa^ 

and  -s-  Pdco  =  f  X  -  --  -  -  y*  9o>. 

3  w        3  do> 

And  therefore 


133.     So  far  our  equation  is  derived  from  physical  con- 
siderations. 


Now  consider  the  auxiliary  function 
u  =  log  Ml 


1  hen  du  =  -77  oh  +  -=— 

dii  dco 

du^L  1 

and  -—  dh  = 

dh 
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because,  since  the  chance  of  the  configuration  o^  .  .  .  z^  is 

e-**^,...  dzK, 
evidently 


hx  h 


Again,  -^dco  =  -  -^r-  -  I    I  e~ 

e-**d*...d*JJ 


the  last  term  being  necessary  because  every  element  of  length 
dxl  ...    is   increased   in    the  ratio   1 
cubical  expansion  above  mentioned. 


dxl  ...    is   increased   in    the  ratio   1  :  1  +  1  -  -   by  the   general 


mi  '  d>U  ^  ,    C^V  X 

1  hat  is  -j —  9ft)  =  —  h  -j^  9ft)  -f-  —  9&>, 

dco  dco          co 

because  9  (dxl . . .  dz^)  =  \  —  , 

ft> 

^  — ^7       j  dy      \9ft) 

and  du  =  —  vdh  —  h-^-\ , 

dco        co 

2\  —  —  Zh^d    —  29    4- 2~dh  (III} 

Substitute  this  in  (II),  and  we  have 
dE 


which  is  a  complete  differential. 

134.  Let  us  now  consider  what  modification  has  to  be  made 
in  the  above  proof  when  Q,  instead  of  being  the  sum  of  squares 
of  the  velocities,  has  the  value  given  to  it  in  this  treatise. 

With  the  form  of  Q  containing  the  b  coefficients,  the 
relative  frequency  of  the  several  configurations  is  different 
from  what  it  was  when  Q  was  a  sum  of  squares  only,  but  the 
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values  of  %  and  -~  are  not  altered  for  any  given  configuration. 

Hence  y  and  -^  have  different  values  from  those  which  they 
do) 

had  in  the  former  case,  v  and  ^  are  the  same  as  in  the  former 

da) 

3         3 

case. 


Nevertheless  the  physical  equation 

dE  =  \dT+dx  +  Pdo>    ..................  (I) 

remains  unchanged  in  form  with  the  altered  value  of  %. 
Further  when  we  substitute  Clausius'  expression  (Art.  24)  for 
Pdco,  the  equation 


,  - 

ax  +  fx--  3      80,  ......  (ii) 

remains  unchanged  in  form   with  the  altered  values  of  h,  %, 

and§ 
dco 

Now  let  us  consider  the  auxiliary  function  u.     We  might 
in  the  proof  above  given  have  put  it  in  the  form 


u  =  log 
In  our  present  problem  it  would  then  take  the  form 


U 


=  log  JJJ  €-*'*  (^  VS6-*«  dx,  .  .  . 
A 

=  log  ff  e-^f- 
JJ         \w 


where  ds  stands  for  dx^...dz^  and  do-  for  d^  ...  efyx.     Then 
h'  and  o>  are  the  two  independent  variables.     Let 

\ 


Then 


X.J 
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,T  du  ~j, 

Now    j-,  dh 


f 

ff  <r 


du  0 
Again   -=-  dt» 


\ 
-        \ 
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,dsd<r  ......  (1). 


the  last  term  referring  to  the  variations  of  all  the  elements  of 
length  (&e,  dy,  ck,  consequent  on  the  expansion  denoted  by  9o>. 

That  is  --  9o> 


'x  Fdsdo- 
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Adding  together  (1)  and  (2)  we  find 


CO 

1 


c  c          sx 

But  now  1  1  Fda  =  —  is  constant  whatever  be  the  values 

of  h  and  co  or  h'  and  co.    It  is  therefore  not  affected  by  the 
variations  of  h'  and  co. 

These  variations  affect  neither  the  quantities  «i...7x,  nor 
the  limits  of  the  do-  integration. 

We  have  then 


a  result  of  precisely  the  same  form  as  in  the  former  case. 
135.     In  any  case  in  which 


the  proof  of  the  second  law  for  the  two  independent  variables  T 
and  to  can  be  presented  in  a  simple  form. 

For  we  have  in  that  case  from  II 


a> 
Now  let  i  —  --J  ,  that  is  i  is  the  time  in  which  a  particle 

moving  with   the   velocity  of  mean  square   would  describe   a 
distance  proportional  to  the  linear  dimensions  of  the  system. 

Then  |8  log  co  =  29  log  i  +  d  log  T, 

and  therefore 


2X3  log  (iT). 


X.]  THERMODYNAMICAL   RELATIONS.  139 

We  might  define  \2iT  as  the  Action  of  the  system  par  excel- 
lence. With  that  definition  the  Action  is  the  same  for  all 
states  of  the  system  defined  by  dE=  0,  that  is  for  all  states  of 
the  system  on  the  same  adiabatic  line.  And  we  may  say  that 

//}  P          \ 
the  thermodynamic  function  <£,  f-y  =  d<j>}  either  is,  or  is  propor- 

tional to,  the  Action. 


136.  But  d%  =  -r^  3&>  cannot  exist  for  any  system  such  as 

we  have  treated  as  a  gas,  except  for  a  so-called  perfect  gas  in 
which  %  =  0.  If  there  be  a  potential  %  at  all,  the  relative 
frequency  of  different  configurations  of  the  molecules  must 
generally  vary  when  T  varies.  And  therefore  although  for  any 
configuration  the  change  of  %  as  G>  varies  is  independent  of  T, 
yet  the  mean  of  all  values  of  %  will  not  be  independent  of 
T,  because  the  variation  of  T  causes  some  configurations,  i.e. 
some  values,  to  be  represented  with  greater,  and  others  with 
less  frequency.  We  may  for  instance  conceive  a  case  in  which 
the  molecules  placed  at  equal  distances  from  each  other  would 
exert  no  sensible  forces  on  each  other  by  reason  of  their 
distances,  and  so  %  =  0.  If  we  set  them  in  motion  they  will 
in  some  cases  approach  so  near  to  each  other  that  ^  =(=  0- 

137.  Professor  J.  J.  Thomson  in  the  work  above  referred 
to  gives  a  proof,  not,  as  he  says,  of  the  second  law  of  Thermo- 
dynamics, but  of  a  proposition  analogous  thereto.     He  does  not 
confine  himself  to  the  two  variables  T  and  co,  but  employs  T 
with   any   number   of    controllable   coordinates   </>.     Then   he 
assumes  that  the  mean  potential  energy  %,  or  in  his  notation 
F,  is  completely  fixed  by  the  controllable  coordinates,  that  is, 
is  not  a  function  of  T.     That  restriction  prevents  Professor 
Thomson's  proof  from  being  applicable  to  gases,  as  a  proof  that 

r)7T 

-^  is  a  function  of  temperature  and  volume,  except  in  the 

limiting  case  of  a  perfect  gas  in  which  V  is  non-existent.  Nor 
does  Professor  Thomson  propose  to  apply  it  to  such  cases. 
It  may  well  be  that  in  the  more  general  class  of  problems 
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treated  of  in  his  work  the  condition  that  V  shall  be  a  function 
of  the  controllable  coordinates  is  not  a  serious  restriction. 

138.  The  following  is  another  form  of  solution  of  this 
problem.  Let  there  be,  instead  of  the  one  controllable  coordi- 
nate G),  any  number  of  them  as  ^  ...  con.  We  will  now  use  to 
as  a  type  controllable  coordinate.  It  is  assumed  that  the  rate 
of  change  of  any  o>  with  the  time  is  so  slow  that  we  may  leave 
the  velocities  Wj  . . .  wn  out  of  account.  Let  Vl . . .  Vn  be  the 
generalised  components  of  momentum  corresponding  to  a)1 . . .  o>n. 
Then  it  is  a  characteristic  of  the  stationary  motion  that  every 

-j     is  zero.     And  therefore  if  p1...pn  be  the  components  of 
cut 

external  force  required  to  maintain  the  motion  stationary 

dT     dXf 

p  =  -= r*  for  each  o>, 

da)     dco 

by  Lagrange's  equations.     See  Art.  24. 

Generally  let  F  be  the  function  denoting  frequency,  that  is 
F  is  a  function  of  the  coordinates  xl ...  zn,  and  the  velocities 
«i"-7n>  °t  the  molecules  of  the  system  such  that  the  time 
during  which  on  average  of  any  very  long  time,  they  lie  between 
the  limits 


zn...zn+dzn, 
a,  ... 


is  Fdxl  ...  d2nda.l  ...  dyn,  or  as  we  may  write  for  brevity  Fdsda. 
Then  the  mean  value  of  any  function  as  %  is 

^=  [[  xFdsda  /  (JFdsda, 


And 
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In  the  same  case  we  have,  if  r  represent  the  kinetic  energy 
in  any  particular  state,  T  the  mean  value  of  T, 


rdFdsda-  + 
=  jjrctFdsdo-  + 

139.  Let  us  at  this  point  assume  the  function  F  to  be  a 
product  of  two  functions,  F  =ff'y  where/  is  a  function  of  the 
coordinates  only.  In  that  case  we  have  also 

X 


%-&*.-  j  tffib/Jjto. 
And  now 


=  2dT  +  j  xdfds  -  jj  rdFdsda, 

because  -=—  9  &>  =  1  1  F  -^-  dcodsd 

do)  jj      aw 


therefore 

^  =  23  log 

O  Tfl 

And  now  we  make  -     a  complete  differential  by  assuming 


where  <^>  and  -^  are  any  functional  symbols. 
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A  solution  of  these  equations  is 

F=  Cfe-fcfr+r),    f=C'e-h*t 
where  G,  C'  are  the  known  constants. 

This  solution  is  the  one  we  should  take  under  the  assumption 
of  condition  A.     Another  solution  is 


f=  G'e~h™,     F  =  C<T1 
if  Q  =  n  T.     This  is  the  solution  corresponding  to  Chap.  VIII. 

The  subject  of  the  Second  Law  of  Thermodynamics  has 
been  very  fully  investigated  by  Professor  G.  H.  Bryan  in  his 
Report  on  the  present  state  of  our  knowledge  of  Thermo- 
dynamics, Part  II.,  British  Association,  1894 


APPENDIX. 


MATHEMATICAL   PROPOSITIONS. 


IN  the  preceding  chapters  I  have  assumed  without  proof 
certain  propositions  relating  chiefly  to  determinants.  They  are 
for  the  most  part  very  simple,  but  as  some  readers  may  think 
they  require  proof,  I  will  devote  this  chapter  to  proving  them. 
I  shall  require  to  make  use  of  the  following  elementary 
property  of  determinants.  If  each  of  the  constituents  in  one 
row  is  equal  to  the  corresponding  constituent  in  a  second  row 
multiplied  by  a  constant  factor,  the  determinant  is  zero.  The 
two  rows  for  which  this  is  the  case  are  said  to  be  similar. 
Whatever  general  proposition  is  proved  of  rows  is  of  course 
equally  true  of  columns. 

(a)     Let  us  take  for  a  type  determinant  of  n2  constituents 


...  C 


in  which  Gpq  is  the  constituent  (whatever  its  value)  in  the  pih 
row  and  qih  column,  and  (712  is  not  necessarily  equal  to  (721,  &c. 
This  can  be  made  to  represent  any  given  determinant  by- 
assigning  proper  values  to  the  constituents  C. 

Let  Dpq  denote  the  minor  of  D  formed  by  omitting  the  pth 
row  and  qih  column,  that  is  by  striking  out  Cpq  and  the  row 
and  column  in  which  it  is  ;  Dqp  the  minor  formed  by  omitting 
Cqp.  In  the  determinant  or  any  of  its  minors  the  diagonal 
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from  the  first  constituent  in  the  top  row  to  the  last  in  the 
bottom  row  shall  be  called  the  axis,  and  its  constituents  axial 
constituents. 

(b)     The  sign  of  Dpq  shall  be  such  that 


=  &c. 

In  order  that  this  may  be  the  case,  we  must  in  forming 
any  minor  Dpq  apply  the  following  rule  of  signs,  namely  to  the 
product  of  all  the  axial  constituents  of  that  minor  the  positive 
or  negative  sign  is  to  be  attributed  according  as  p  +  q  is  even 
or  odd,  and  we  deduce  therefrom  the  sign  of  the  product  of  any 
other  set  of  constituents  of  the  minor  according  to  the  usual 
rule  in  expanding  determinants.  In  like  manner  Dprqs  shall 
denote  the  minor  formed  by  omitting  Cpq  and  Crs  and  the  rows 
and  columns  in  which  they  are.  It  follows  from  the  above 
convention  concerning  signs  that 

p=ng— n 

nD=  2    2  GpqDpq. 

p=l  q=l 

Then  evidently  -^—  =  Dpq,  whatever  values  p  and  q  may 
have. 

(c)     Now  consider  the  new  determinant 

C   —  6  x       C    6 


•  •  •  ^nn      vnXn. 

If  all  the  (7s  in  any  two  rows,  as  the  pih  and  qih,  were  zero, 
jy  would  be  zero,  because  the  pih  and  qih  rows  would  then  be 
similar.  It  follows  that  every  term  in  D  contains  at  least 
n  —  1  C  factors,  and  therefore  that  D'  contains  no  products  of 
any  two  or  more  of  the  new  partial  constituents  6x.  For  such 
products  necessarily  enter  in  pairs  equal  and  of  opposite  sign, 
as  for  instance 

Therefore 
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(d)  Now  form  a  new  determinant  of  (71+  I)2  constituents 
by  adding  to  D  an  (n  +  l)th  row  and  (n  +  l)th  column.  The 
(n  +  l)th  row  shall  be 


and  the  (n+  l)th  column  shall  be  symmetrical  with  it,  contain- 
ing ft(M+1)  in  the  first  row,  £a<n+i)  in  the  second,  and  so  on. 
Also  let  f3pq  =  pqp,  for  all  values  of  p  and  q. 

Let  A  be  that  new  determinant,  namely 

11  >    ^12  •  •  •  ^in  >    Pi  (n+i) 

21,     #22...    #271 1 


Next  in  D'  write 


_        (n+i)  P2  (n+i) 


& 


2an+1  ' 

/J  £»2  (W+l)  /) 

"2^2  =  "o >      "p57?  = 

Then  we  have 

Q 

D'  =  D-22' 


« 

' 


g  (n+i) 


2an+1 


and  by  the  ordinary  expansion 

A  =  2an+1D- 

Therefore  D'  = 


(/)     In  this  last  deduction  we  have  attributed  no  particular 
values  to  the  C  constituents. 

It  would  therefore  hold  for  all  the  minors  of  D'  not  contain- 
ing the  suffix  n  -\- 1  as  well  as  for  D'  itself,  namety 

1 


An, 


*•'-* 


B+1 


and  therefore 


IK  "    A  ' 
^-A» 

rr  ~  A  ' 


10 
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(g)  Now  consider  the  function  e~h®  in  which  h  is  any  con- 
stant, and  Q  is  a  quadratic  function  of  n  variables  uly  u2  ...un, 
namely 

Let  us  effect  the  integration 

•~h®dun  between  limits  +  oo  . 
The  terms  in  Q  which  contain  un  are 

which  may  be  put  in  the  form 

&c.\2     (blnu!  +  b^u2  +  &c.)2 


( un  \/an  +  — 

\  2  van 

Let  us  for  a  moment  write 


&c. 


=  z. 


Then 


f   6-*Qdlln 

J  —oo 


~ftQ'  multiplied  by  a  numerical  factor, 


and 


Comparing  these  results  with  (e),  (/)  we  see  that  if 
2«i,  bl2,  bK...bm, 


D  containing  (n+  I)2  constituents 
and      D'  = 


l»- 


containing  %2  constituents ; 
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then  W*L       AL'_A, 

jy     D  *    iy     D' 

That  is,  these  ratios  are  not  altered  by  the  change  in  the 
coefficients  o^,  612,  &c.,  consequent  on  the  integration  accord- 
ing to  un.  Evidently  the  same  law  holds  for  every  subsequent 
integration.  Let  us  then  perform  the  integrations 


and  let  the  result  be 

...dun 


If  e-tUrt+BiM+AM  dujfa  is  proportional  to  the  chance 
that  two  variables  ^  and  u2  shall  lie  between  the  limits 
^...u^  +  dUifor  the  one,  and  u*...u*  +  dui  for  the  other,  we 
easily  find  the  mean  of  %2,  thus 


the  limits  of  integration  being  for  each  u  ±  oo 


if  d  be  the  determinant 


But  -^  =  -jj-  by  the  last  article. 
Therefore  u?  —  ^  • 

rljj 

T\ 

Similarly 


and  so  on. 

10—2 
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(h)     Under  the  same  circumstances  to  find  the  mean  value 

of  -=-£ ,  or  2^!  +  &12^2  +  blsu3  +  &c.,  M!  being  given  (see  Art.  97). 

That  is  to  effect  the  integration 

&c.). 


HI- 


Let 


2  + 


The  coefficient  of  u^  in  P  is  twice  the  coefficient  of  u^  in  Q. 
The  coefficient  of  every  other  u,  as  u2,  in  P  is  equal  to  the 
coefficient  of  u^  in  Q, 

After  integrating  according  to  un  the  coefficient  of  u^  in  P 
is  2ax  -  ^  while  that  of  u*  in  Q  is  a,  -  ^ . 


The  coefficient  of  uz  in  P  is  612  —  ™  m,  and  the  coefficient 
of  ^^2  in  Q  is  the  same. 


The  same  relation  between  the  coefficients  in  P  and  those 
in  Q  which  existed  before  integration  according  to  un  exists 
also  after  that  integration.  The  same  rule  holds  for  every 
subsequent  integration  to  u2  inclusive.  But  when  all  the 
integrations  have  been  effected  to  w2  inclusive,  Q  has  become 


Therefore  P  has  become         u    which  is  the   mean 


value  sought. 

(i)     Let  us  now  consider  the  determinant 


D= 


-/A  -/A... 

-/A  (!+/,'),      /,/„       /,/„... 

-/A    /,A     (!+/,').  A/.,- 

-/A      /./a,         A/s,       (1 


treated  of  in   Chapter  v,  Art.  63.     Here  5  is  written  for 
of  that  Article. 
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Comparing  it  with  the  type  determinant,  we  find 

*»-. 

&C. 

Let  it  be  required  to  find  the  values  of  D  and  its  first 
minors. 

(j)  We  observe  that  in  D  every  row  except  the  top  row  con- 
tains one  constituent  of  the  form  1  +/2,  e.g.  the  (p  +  l)th  row 
contains  1  +fp.  If  it  were  //  instead  of  1  +//,  the  (p  +  l)th 
row  would  be  similar  to  the  top  row  and  therefore  D  would  be 
zero. 

If  we  strike  out  any  column  other  than  the  first,  say  the  qih 
column,  we  deprive  the  qth  row  of  its  constituent  1  +/9-i2. 
Hence  the  ^th  row  is  now  similar  to  the  top  row. 

And  if,  to  form  a  minor,  we  strike  out  any  other  row,  neither 
the  first  nor  the  qih,  the  determinant,  i.e.  the  minor  so  formed, 
is  zero.  We  see  then  that  if  p  >  1,  and  q  >  1,  and  p  =(=  q,  Dpq  =  0. 
But  if  we  expand  D,  the  constituent  Cm  appears  only  in  the 
form  DpqCpq,  and  Dpg  =  0.  Therefore  D  when  expanded  does 
riot  contain  any  anaxial  constituent  except  those  in  the  first 
row  or  first  column. 

From  (i)  and  (j)  we  see  that  /)=  $2.  For  the  same  reason 
D^=S2,  and  so  on  for  every  coaxial  minor  except  Du,  but  as  we 
shall  see  later  Du  =  2/2  +  1. 

(k)  We  have  next  to  find  D12,  D13,  &c.  Striking  out  from 
D  the  top  row  and  second  column  to  form  D12,  we  find  that  in 
D12  every  row  except  the  top  row  has  one  constituent  of  the 
form  1  +/2.  Therefore  by  the  same  reasoning  as  in  (j)  and  by 
the  rule  as  to  the  sign  of  the  product  of  the  axial  constituents, 
A2  =  Ai  =  +/i&  When  we  form  D13,  the  product  of  the  axial 
constituents  is  by  our  rule  to  be  positive.  But  in  order  to 
reduce  this  minor  to  one  in  which  the  top  row  does  not 
contain  a  constituent  of  the  form  1  +/2,  we  must  interchange 
the  two  top  rows,  and  thereby  change  the  sign,  of  the  minor. 
Therefore  Dn  =  +/2$  and  so  on.  In  forming  any  minor  D^  we 
have  if  p  is  odd,  an  odd  number  of  changes  of  sign. 
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(I)     Next  let  us  deal  with  Du. 


Let 


l+/i2,    /i/a,    A/3,... 


/I/.,       /I/., 


or 


All  products  of  the  form  AV?  evidently  disappear  in 
expansion,  and  therefore  all  terms  involving  y*2  except  in  the 
first  degree  disappear,  and  therefore 


=  S/»  if/,,-1. 

(m)  Again,  if  we  strike  out  any  column  except  the  first,  say 
the  second,  we  thereby  deprive  the  second  row  of  its  constituent 
1  +/22.  And  therefore  if  we  now  strike  out  the  top  row  we 
shall  have 


,'),      /s/4,  ... 


=  —  /!/2  by  the  rule  as  to  the  sign  of  the  axial  constituents. 
Similarly  A13  =  -/i/3  . 

Also  A.^22-,2, 


(^)  In  the  type  determinant  let  now  every  C  become 
C7+9C7,  and  let  all  the  30's  be  small,  so  that  we  may  neglect 
products  of  them.  That  being  the  case 


2     Z, 

=  l  q=l 


'PI' 
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Let  us  apply  this  to  the  determinant  of  (*)  in  which 

cu  =  &,      cl2  =  -/A      cn  =  -/A  &c., 
?,  &c- 


Here  as  we  have  seen  -^    =  0,  unless  either  p  =  1,  or  q  =  1, 


Therefore 

9/>  =  DudCu  +£&€„  +f£dC13  +  &c. 


Now  as  in  Art.  83  let  us  make 

dcu  =  o,         dcl2  =  ±blns,       dC13  =  i&^s,  &c. 

8021  =  J6ln>Sf,       90^  =  16^,  &c., 
aOa  =  -  6ln/1}    9(733  =  -  6^,  &c. 

Then  as  products  of  the  6's  are  to  be  neglected,  we  have 

-  &c. 


=  0. 

(o)     Let  us  treat  Du,  or  A,  in  the  same  way.     Here 

Cu=l+/i2,      012  =/,/,,   &C. 

And  8Ai  =  9A, 


And  8(7U  =  -  6ln/!,     aO^  =  -  62n/2,  &c., 

9012  =  i(612  -  6m/2  -  6^/0,  &c. 
Therefore,  since  each  anaxial  term/j/2,  &c.  appears  twice, 

a  A  =  -  (2/*  -/^  6^  -  (I/*  -/22)  62n/2  -  &c. 


+  &C.) 
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That  is  9A  =  —  612/i/2  —  613/1/3  -  &c. 


In  this  result  every  bpgfpfq  is  to  be  counted  once  only,  that 
is,  we  do  not  distinguish  bpq  from  bqp.  If  bpqfpfq  and  bqpfqfp 
are  counted  separately,  we  must  write 


and  this  form  is  the  one  required  in  Art.  83. 

(p)    To  prove  that  for  certain  values  of  b,  %fp%bpgfq  —  ^fp^b 
see  (33)  of  Art.  83.     Let /(Art.  61)  =  — ^  when  r  >  a. 

Taking  0  for  origin,  let  p  be  any  point  Op  =p.     About  p 
describe  a  sphere  of  radius  r  less  than  p,  and  let  q  be  a  point 

on  its  surface.     If  n  +  2  =  1,  or  /=  - ,  the  mean  value  of  fq  for 

all  points  q  on  the  sphere  is  ^j~  >  °r/>-     And  equation  (33)  is 
accurately  satisfied,  so  long  as  r  <  p. 

If  w  +  2  >  1 ,  we  have  in  usual  notation 

sin  OdO 

o  ~^+* 


- 

.  o 


-  r8) 


n 


r 

=/;  +  terms  containing  fp  multiplied  by  —  2&c. 


And  if  now  b  =  —  f  -  J  ,  where  fju  >  3, 


fp 

c 
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s* 

because,  -  being  negligible,  all  the  terms  in  the  integral  which 
follow  the  term  fp2  £6  are  negligible  compared  with/p2E6. 
(q)     To  find  the  condition  that  the  quadratic  function  (q) 

shall    be    always    positive  whatever   values   be   attributed   to 
ult  u2...  un. 

Evidently  since  u2 ...  un  may  be  all  zero,  at  must  be  positive. 
Similarly  «2,  a3,  &c.,  must  be  positive. 

Also,  given  un,  Q  is  minimum  when 

dQ  _. 
dun 

And,  when  minimum,  Q  must  be  positive.     Substitute  in  Q 
for  un  the  value  found  for  it  from 

'^2 +  &C.  =0, 

that  is  un  =  -  ™^-'^«+"-  t 

Q  then  becomes 

-&c., 


where  a^^-,     a/  =  aa-         ,  &c., 

4anJ  4an' 

and  every  a7  must  be  positive.     Also  612'  =  612  —  ^LJ?  &C- 

^ttjj, 

Proceeding  in  the  same   way  Q  is  ultimately  reduced  to 

Q  =  ^yr-  ^j2,  and  as  this  must  be  positive,  yr-  must  be  positive. 
Z//n  *TI 

Similarly  -=-  must  be  positive,  and  so  on.     Therefore  D  and  all 

^22 

its  first  coaxial  minors  have  the  same  sign.     But  if  we  only 
reduce  Q  to  a  quadratic  function  of  ut  and  u.2 
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From  which  and  the  corresponding  equations  we  see  that 
all  the  second  coaxial  minors  have  the  same  sign  as  the  first 
coaxial  minors,  and  therefore  the  same  sign  as  D. 

Similarly  we  can  prove  that  the  third  coaxial  minors  must 
have  the  same  sign  as  the  second  and  so  on,  so  that  D  and  all 
its  coaxial  minors  must  have  the  same  sign.  And  this  sign 
must  be  positive  because  al}  a2,  &c.,  are  coaxial  minors.  This 
is  the  required  condition. 

(r)  Referring  to  Chapter  v,  in  which  it  is  proved  that  with 
(r)  spheres  of  finite  diameter  c, 

dMxz  1  fd£ 

dt  2h  \dz     dx 

where  Mxz  or  M  =  \\  Idxdydzay  '  ,/  +  T  )  > 

to  prove  the  same  thing  for 


Consider  the  space  between  two  infinite  planes  x  ==  0  and 
x  =  dx,  between  which  ~  is  supposed  constant  and  negative. 

Let  u  =  a.  —  %.  Then  it  is  proved,  as  in  Art.  64,  that  the  number 
per  unit  of  volume  of  molecules  for  which  a  —  f  lies  between  u 
and  u  4-  du,  when  x  =  0,  exceeds  the  corresponding  number  when 
x  =  dx  by  the  quantity 

C<rv*Zh'u^  dxdu. 

dx 

Therefore  the  number  which  enter  the  space  between  the  planes 
exceeds  the  number  which  leave  it  per  unit  of  time  by 


Hence  as  in  Art.  66,  it  follows  that,  for  each  element  of 
volume  dxdydz, 
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d  (ua)  =  -  ( r    Ce-^e-^Zh'iWduda  ^ 
at  JJ  -a,  CUB 


dx' 

Now  u  and  a  are  not  independent,  as  u  and  7  were  in 
Art.  66.     Hence 

'•)L  C)L       O  O 

~  =  —  »*  =  _  A  =  _^ 

OZ,/  u  OL'    yl.ta         .I/.'.T/    > 


and 


Hence  dfdu*_dtd*_     _3 

^  ^  "cte'cft  "      2/i  ...............  (  >J 

which  corresponds  to 


of  Art.  66. 

It  follows,  the  motion  being  stationary,  that 


= 

cfc  dx     2h  \dx    ' 
or  by  the  process  of  Art.  68 


_= 

r  dr  ^  2h  \ 

When  the  material  points  are  replaced  by  spheres  of  finite 

diameter  c,  -7-  acquires  a  new  term  due  to  collisions,  which  we 
cLt 

will  denote  by  ^-  ,  and  we  have 
vt 


dx'dt      dx    dt    ' 

Vx>     Vyt     Vz,    being    the    components   of    the    half    relative 
velocity  of  two  colliding  spheres. 
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r)F2 

Also  ^r  =  27rc2pF(\'2  -  X2)  V\ 

dt 

where  \V=  VX)     pV=Vy,     vV=Vz, 

and  using  the  same  notation  as  in  Arts.  70  —  71,  we  find 


sin  0  cos  0 


I  F3X  f  2d0  f 
dt        TT  dx"l          J  o      J  o 

x  (-  X  cos  0  +  Vl  -  X2  sin  6  cos  </>')  |X2  (4  cos4  6  -  4  cos2  0) 


+  (  1  -  A,2)  sin2  20  cos2  (/>'  -  2\  Vl  -  X2  cos  2(9  sin  20  cos  0xj 


4 

whence  introducing  the  factor  Ce~2hV2V2dV  where  C  =  -j- 

and  integrating  from  V  =  0  to  V  =  oo  , 

9  7*"      8  1  df 

we  get  -5—  =  ^  7rc3/o  ^-  -=2  , 

8^        15       ^  2hdx 


We  have  next  to  calculate  the  term 


due  to  collisions.     That  gives  by  the    same    method  as  that 
employed  in  Art.  76 


[2 

Jo 


>   sn     cos 

dt  dx         TT  J        Jo 


x  (-  X  cos  d  +  Vl  -  X2  sin  d  cos  </>') 
x  (-  2X  cos2  6  +  Vl-X2  sin  20  cos  </>')  i  ^(r  =  a). 
Whence  we  obtain  for  a  pair  of  colliding  spheres 


45 
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and  for  each  sphere  separately 
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Adding  this  to  (3)  we  find  for  the  change  of  Mxx  with  the 
time  due  to  collisions 

so 


Using  Mxx  for  Mxz  or  M  in  Art.  79,  we  find 

-2dd£_3_(d];]* 
dt  dx      2h  \dx) 


The  first  line  is  zero.     The  second  can  be  made  zero  by 
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